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Final Management Report 

ADP: ‘Synoptic Study of Quasar Continua* 

P.I.: Martin Elvis 

This report summarizes the work undertaken on the above ADP grant. The 
work is now complete, with the submission of the major product, the Atlas of 
Quasar Energy Distributions’ (Appendix A) being finalized for submission to the 
Astrophysical Journal Supplements. Other papers written as a result of this grant 
are included in Appendix B. The overall themes developed under this grant will b» 
pursued in a Long Term Space Astrophysics Program grant, which will build on 
the data and methodology established under the ADP. The aim of the Long Term 
program will be to extend our knowledge of the quasar continuum shapes to high 

redshifts. 

Once the QED (Quasar Energy Distributions) data base became sufficiently 
complete we began our ‘synoptic’ ( i.e . multi- wavelength) analysis. This bore fruit 
immediately in the study of the ‘accretion disk’ Big Bump. Our first detailed study 
of the strength of the opt/uv bump strengths and shapes was undertaken. Four 
different measures of bump strength were made. These all correlate, albeit with 
scatter. All the measures clearly showed a range of bump strengths, most of which 
can be reproduced by reddening. A simple change in normalization of a ‘bump’ 
component also cannot reproduce the observed range of shapes. The resulting 
paper on 'Weak Bump Quasars’ was published in Ap. J. Letters. Follow-up work 
looking at the whole range of ‘Big Bump’ strengths, and shapes followed. 

Accurate Galactic column densities derived from Green Bank 140ft data were 
published in the Astronomical Journal. These column densities are some 10 times 
more accurate than those derived from the standard ‘Bell Labs’ survey of Stark et 
al. alone, and allow valuable constraints in fitting soft excesses to x-ray spectra. 
The work also has value for soft x-ray and EUV surveys (eg ROSAT, EUVE). 
Further work on this subject is planned due to the collapse of the 300ft, which 
would otherwise have undertaken an all-sky survey to this level of accuracy . 

A paper on the variability of the soft excess in PG1211 + 143 resulted in a 
paper to the Ap. J. It is clear that the soft excess in PG121 1 + 143 is persistent over 6 
years (from an early, serendipitous, Einstein observation through the last EXOSAT 
observation), yet varies in amplitude on a short (weeks) timescale, demonstrating 
that it comes from a compact source region. If the variation is intrinsic, ie not due 
to changing absorption, then the emission cannot be from optically thin plasma. 


A visit to Meudon to complete work on the ‘soft-excesses’ paper with 
Dr. Masnou was undertaken. The visit, by B. Wilkes, resulted in considerable 
progress on the ‘soft-excesses’ paper with J.-L. Masnou. Dr. Wilkes and Dr. 
Masnou characterized the ‘ultra-soft’ excesses in a more definitive fashion, modifying 
the conventional fitting procedure to make the strength of the excess the fitted 
parameter not the normal ‘power law index’. The intercalibration of the Einstein 
MPC/IPC was investigated (with J-L. Masnou, Meudon, FRANCE) in order to 
confirm the new X-ray spectral results derived from this instrument combination. 
A paper on the soft X-ray excesses is in press with Astronomy and Astrophysics. 

The EXOSAT ME and LE spectral data on 12 quasars (the complete set of 
good quality EXOSAT quasar data) were analyzed by A. Comastri (visiting from 
Univ. Bologna, ITALY). Results show interesting differences from the IPC results 
and use of the Boron filter suggests in a few cases that the high energy spectral 
slopes continue clown to ~1 ke\ . The 3LX and A1P filter data show strong excesses 
in most cases. A paper is in press in the Ap. J. 

Our sub-millimeter observations from the JCMT were used in conjunction 
with the IRAS and other, ground-based, infrared data to constrain synchrotron and 
dust models for the infrared continuum in radio-quiet active galaxies and quasars. 
A paper has been published in MNRAS. 

X-ray data for the Piccinotti sample were included in the data base. Optical 
and ultraviolet data, from our collaborators C. Boisson and M.J. Ward, were 
included for the Piccinotti sample. The data base software was enhanced to include 
basic model fitting functions (polynomial fitting to the data, integration under a 
curve, power-law and black body fitting) to allow a more quantitative investigation 
of the strength of the Big Bump, in particular to define the ‘weak bump’ quasars 
more carefully. 


The ‘TIGER’ software was updated to use the new release 2.0 of SMONGO, 
which allows the use of color plotting. This makes our crowded data plots far easier 
to examine by distinguishing, e.g. different epoch observations from one another, or 
data from model fits. The SUN workstation color upgrade arrived and was installed. 
It has proved valuable in visualizing the multi-epoch data and in comparing data 
with models. 

A Harvard astronomy department graduate student, Olga Kuhn, elected to 
work on the range of ‘Big Bump’ strengths and shapes for her research exam. 
She is concenti ating on modeling the optical/ultraviolet ‘Big Bump’, in particular 
examining the iange of accretion disk properties needed to explain the variety of 
shapes that are seen. She modified her accretion disk code to incorporate the 



effects considered by B. Czerny which will allow more sophisticated modeling of the 
optical/ultraviolet ‘Big Bump’ color-color plots. She then made ‘mixing curves’ of 
these models with various fractions of power-law slopes with a number of spectral 
indices. These were compared with the IR/V/UV color-color plot for the data (as in 
the Weak Bump Quasars paper), and with other color-color plots, e.g. IR/UV/SX. 
0. Kuhn attended the SAAS-FEE School (Switzerland) on ‘Active Galactic Nuclei’, 
a week-long course of lectures given by leaders in the field designed for graduate 
students. Her research exam is due soon and her write-up for that exam will be 
recast as a paper for the Ap. J. 

The accretion disk modeling code being developed by O. Kuhn and J. 
McDowell, based on the original code of B. Czerny and was used to generate the 
model locus described above. The code has shown that the effects of General 
Relativity and of electron Scattering do not change the shape of this locus, but 
do strongly affect the derived physical parameters (mass and accretion rate). This 
implies a degeneracy in these parameters that the present data cannot lift. The 
effects of inclination and of using a Kerr metric for the black hole are now being 
explored. 

Two summer students, Elisha Polomski and Sally Oey, were employed. Their 
main task was the reduction of the optical continuum and line data. They made 
major progress in reducing the optical data that is part of the data base. Many of 
the optical spectra are now ready to be included in the TIGER data base. Several 
more optical observing runs remain to be reduced however. 

The peculiar object PG 1407+265 has an uninterpretable spectrum. We no 
longer believe that we know its redshift; we are not even certain that it is a quasar. 
A paper is in preparation describing this unusual object. 

A visit to the IUE RDAF at GSFC by J. McDowell resulted in the complete 
reduction of some 60 IUE spectra of quasars in our sample. This completed our data 
base. It also added information on time variability in the ultraviolet which is vital 
to our confidence in our non-simultaneous multi-wavelength energy distributions. 
The sample of 37 quasars with ‘complete’ continuum coverage has been analyzed 
to examine the shapes of the optical/ultraviolet ‘Big Bump . A new color-coloi 
diagram was constructed using octave wide near-IR, optical and UV bands. This 
diagram allows the distribution of the slope/strength and the curvature of the Big 
Bump to be displayed for the whole sample in a model independent fashion. The 
locus of points described by accretion disk models can be plotted on this diagram, as 
can the effects of dust reddening, the inclusion of a galaxian (starlight) component, 
and of a non-thermal power-law component. The sample objects are restricted to 
one part of this diagram. The accretion disk locus forms one boundary to this 


region, suggesting that it has some physical significance. However the mixing of 
disk models with reddening, galaxy starlight or power-law components does not 
naturally reproduce the distribution of data points. This promises to be a valuable 
technique for understanding this major component of the quasar continuum. 

The complete sample of Quasar Energy Distributions has been assembled 
with care taken over variability and inter-wavelength calibration. Host galaxy 
contributions to the quasar continuum light have been estimated from the literature 
and subtracted. The set of QEDs are being made available in digital form via 
anonymous ftp from SAO. 

As a result of this study of quasar energy distributions we find that: 

1. The bolometric luminosity of the quasars is a factor of 12.7 ± 4.2 times the 
monochromatic visual luminosity, using a simple interpolation in the EUV; 
the more accurately defined ‘UVOIR luminosity’ is a factor of S.5 ± 2.9 times 
the monochromatic visual luminosity. 

2. The strength of the 'ultraviolet bump’ varies by a factor of 10, but typical ob- 
jects have an ultraviolet luminosity three times their near infrared luminosity. 

3. "The submillimeter break” often occurs at wavelengths of 60pm or shorter. 
The IR emission can be weak relative to the optical. 

4. Objects with weak ultraviolet bumps are luminous in the far infrared. 

5. We confirm the near constancy in wavelength of, and ubiquity of, the “1 pm 
inflection” first noticed by Neugebauer et al (1979) and studied by Elvis et al 
(19S6) and Neugebauer et al (19S7). 

6. The shape of the ultraviolet bump varies significantly from object to object, 
rising steeply in the optical for some objects, but remaining flat there for 
others. 

1 ■ The ultraviolet variability of these objects (chosen to exclude strongly opti- 
cally variable quasars) is typically less than a factor of 2 on a timescale of 
years, and many vary by less than twenty per cent. 

8. While we confirm the correlations found by Kriss (1988), we suspect that they 
may be induced by host galaxy contamination at low luminosities. 

9. The ‘infrared to x-ray connection’ discussed by Carleton et al (1987) is not 
present in this sample. 

The publication of this Atlas completes the work under the ADP Grant on 
the 'Synoptic Study of the Quasar Continuum’. 
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Abstract 


We present an atlas of the spectral energy distributions (SEDs) of 33 normal, non- 
blazar, quasars over the whole available range (radio to 10 keV X-rays) ol the 
electromagnetic spectrum. XX radio-quiet and XX radio-loud quasais aie included. 

We discuss the effects of the limited dynamic range of flux density in the 
different parts of the spectrum on the range of SED shapes that aie seen. 

We present a mean energy distribution for radio-loud and radio-quiet quasars, 
along with an analytic representation appropriate for theoretical modelling. 
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Abstract 


We present an atlas of the spectral energy distributions (SEDs) of 33 normal, non- 
blazar, quasars over the whole available range (radio to 10 keV X-rays) of the 
electromagnetic spectrum. XX radio-quiet and XX radio-loud quasars are included. 

We discuss the effects of the limited dynamic range of flux density in the 
different parts of the spectrum on the range of SED shapes that are seen. 

We present a mean energy distribution for radio-loud and radio-quiet quasais, 
along with an analytic representation appropriate for theoretical modelling. 


1 Introduction 


One of the reasons that main sequence stars are much better understood than quasars 
is that they radiate (almost) black body spectra at between 4000 and 30,000 K, 
so that their black body peak moves conveniently through the optical band. The 
resulting strong color changes allowed the early recognition of the main sequence in 
the Hertzprung-Russell diagram. By showing that most starts lay in a restricted band 
of color and luminosity this diagram provided a crucial input to theories of stellar 
structure. The current lack of understanding of quasars may similarly be due to 
the distribution of their continuum light. Because the quasar phenomenon covers an 
extremely broad range of wavelengths it is hard to see continuum features analogous 
to the black body peak in normal stars. Quasars seem in fact to emit almost constant 
power per decade of frequency from lOO/im to at least lOOkeV (see e.g. figure 1 of 
Carleton et ah, 1987). While this equipartition is surprising (and may be to some 
extent an observational artifact, see §6) it contains too little information to constrain 
theoretical ideas. Overcoming this problem requires the assembly of Spectral Energy 
Distributions (SEDs) for sizable samples of quasars over the whole accessible range 
of the electromagnetic spectrum from far-infrared to hard X-rays. 

In this paper we present SEDs for a sample of 33 quasars. We concentrate 
on ‘normal’, (he. non-Blazar)quasars. Our primary selection criteria were: existing 
Einstein observations at good signal-to-noise ( to ensure good x-ray spectra); an 
optical magnitude bright enough to make an IUE spectrum obtainable. This 
combination leads to some IRAS upper limits. We consider the effects of these 
selection criteria in some detail in section 4. 

Several previous studies of quasar SEDs have been published. Each emphasized 
a particular region, and none considered the X-rays in detail. Infrared to optical 
for the PG quasars were presented in Sanders et al (1989) with an emphasis on 
explaining the infrared ; near-infrared to ultraviolet SEDs for IUE observed quasars 
were presented by Sun and Malkan (1989). Our study is different in that it includes 
X-ray data for all the objects, divided into three energy bands in many cases, and 
takes care to include both IUE and IRAS data. The sample is fairly evenly divided 
between radio-quiet and radio-loud objects. 


2 Continuum Features 

In Fig. 1 the radio to x-ray rest frame energy distributions of two typical quasars, 
3C273 and Mkn586, are shown. All energy distributions are plotted as log i/f„ vs. log 
such plots give the best indication ol the frequency ranges where most energy is 
released (since j/f„ is the flux per logarithmic frequency interval). The x-ray data are 



indicated with a ’bow-tie’ symbol representing a power law fit with both the best fit 
slope and the 90 per cent confidence limit slopes drawn in. No data are available in 
the extreme ultraviolet 'gap' beyond the Lyman limit where our galaxy is opaque. 

In the 1-100 gm infrared both quasars are almost flat (Ward et al. 1987, 
Neugebauer et al. 1987, hereafter N87). A single, nearly horizontal, power-law fits 
the IR points reasonably well and intersects the x-ray point at about 1 keV. We will 
call deviations from this power-law ‘continuum features’. There are four prominent 
features of this kind. 

The power output always drops in the millimeter band ( the ‘ mm-break , Fig. 
1) but the size of the drop varies dramatically from object to object. Quasars in which 
the drop is only 2 decades are called ” radio-loud” (e.g. 3C273). The great majority 
of quasars have a much stronger mm-break of 5 or even 6 decades (Condon et al 1981, 
Kellermann et al. 1989) and are called ‘'radio-quiet” (e.g. Mkn586). This distinction 
between radio-loud and radio-quiet quasars is the oldest in the quasar literature and 
goes back to the ‘blue interlopers’ found in early radio source identification work 
(Sandage 1965). Radio-quiet objects are much the more common, by about a factor 
10. The mm-break is the strongest feature in normal quasar continua. 

The optical-ultraviolet continuum rises above the infrared and forms a “UV 
bump” (Shields 1978, Malkan and Sargent 1982, N87). Variability studies show that 
this is a separate component lroin the infrared (Cutri et al. 1985) since it varies much 
more strongly. This big bump is most often interpreted in terms of thermal emission 
from an accretion disk ( e.g. Malkan 1983, Czerny and Elvis 1987). The beginning 
of the bump is marked by an inflection between 1 and 1.5 //m in the rest frame, this 
“1 n m inflection” is the only continuum feature whose wavelength is well defined. 

X-ray spectra ol many radio-loud quasars and high luminosity Seyfert 1 
galaxies rise to high frequencies in log j/f„ vs. log u energy distributions (Mushotzky 
1984, Bezier et al. 1984. Turner and Pounds 1989, Williams et al 1991). They cannot 
then be an extension of a flat, or slightly falling, infrared power-law, as has been 
suggested for the radio-quiet objects (Carleton et al 1987). A new emission component 
must be emerging in the x-rays in these objects. This ‘hard x-ray component at 1 
keV is the third clear feature of normal quasar continua. It is not always seen. 

The most recently identified continuum feature is the ’XUV excess’ (e.g. 
Arnaud et al 1985, Wilkes and Elvis 1987, Turner and Pounds 1989, Masnou et 
al 1991). The 1 keV spectrum of quasars is often fitted with a power law spectrum 
since the available spectral resolution ’.s very low; studies using the Einstein IPC and 
EXOSAT showed that excess flux above this power law was often present in the ultra- 
soft (0.3 keV and below) region of the spectrum. It is possible that the XUV excess 
is the same physical component as the ultraviolet bump since the energy distribution 
rises towards the EUV from both sides of the unobserved spectial legion, it is piesent 
in about half of objects (both Seyferts and quasars) studied. 


3 The IPC UV sample 


We have carried out an extensive program of quasar continuum observations over 
the past five years, the results of which we present here. Since we believe it to be 
important to include x-ray observations, we restricted ourselves to objects that had 
been observed with the Imaging Proportional Counter (IPC, Gorenstein, Harnden 
and Fabricant 1979) instrument on the Einstein Observatory (HEAO 2, Giacconi et 
al 1979) satellite. We note that this introduces a bias towards objects which are low 
redshift, moderate luminosity and have strong x-ray to optical flux ratios; the sample 
is heterogenous and not flux-limited. The objects were selected mainly from the PG 
and 3C and Parkes catalogs. No single completeness criterion could be used however, 
because of the quasi-random way in which the original Einstein observations were 
made. 

We have selected a subsample of quasars which have sufficient counts in the IPC 
to give a reasonably constrained power law spectral fit, and which are optically bright 
enough (V< 17) to be observable with IUE: the TPC UV sample’. Approximately 
half the sample is radio-loud and half radio-quiet. The x-ray and ultraviolet data 
have been combined with IRAS space-based infrared data and ground-based optical, 
infrared and radio observations to construct the complete energy distributions. 

The 33 quasars in the Einstein UV' sample are listed in Table 1(a). The 
sample has substantial overlap with the work of Elvis et al (1986) and Wilkes and 
Elvis (1987). The table gives the common name of the quasar as well as its catalog 
number in the Einstein Observatory Source Catalog (Harris et al 1991), and the 
name of the associated host galaxy where appropriate. Both B1950 and J2000 
coordinates are provided, as well as the redshift and typical V magnitude. The 
estimate of the foreground Galactic hydrogen column density (in units of 10 2 ° cm -2 ) 
is listed together with a reference. Most of the column estimates are accurate 
values from Elvis, Lockman and Wilkes (19S9). Finally, each object is given a 
classification. The classifications are as follows: radio-quiet (RQ), broad absorption 
line (BAL, a subset of RQ), and radio-loud (RL). Radio-loud quasars are further 
subdivided into superluminal (SL, a subset of RL), flat spectrum compact (FSC), 
steep spectrum compact (SSC), and Fanaroff-Riley class 2 steep spectrum doubles 
(FR2). Further, following the convention of Veron-Cetty and Veron (1987) and 
Schmidt and Green (1983), radio-quiet objects with absolute visual magnitude fainter 
than —23.0 calculated according to their prescription are designated as Seyfert 1 
(Syl); two objects in the UV sample satisfy this criterion. Table 1(b) lists the 
corresponding properties of 37 other quasars for which ultraviolet, optical and infrared 
data is presented in this paper. These other quasars were observed by Einstein or 
EXOSAT but are not included in the main sample either because of poor quality 
x-ray data, strong optical variability, low luminosity or low optical flux. 

Further properties of the quasars in the UV sample are listed in Table 2, and 


illustrated as histograms in Fig. 2, with the shaded regions corresponding to the 
radio-loud objects. The numerical value of the radio-loudness, Rl, is defined to be 
the logarithm of the ratio of observed core 5 GHz flux to flux in the optical B band, 
and an object is considered to be radio-loud if Rl = log(/5G//z//s) > 1- Note that 
this definition is based on observed frame fluxes, but corrections to the rest frame 
do not make an important difference. In each case an attempt has been made to 
estimate a core flux on arcsecond scales, since we are studying the properties of the 
compact central source at other wavelengths and wish to neglect the extended radio 
source, if any. We were unable to find maps of the two most southern sources, and 
the VLA fails to resolve the core of 3C48 from the compact steep spectrum source 
in which it is embedded (Spencer et al 1989). For sources where information on the 
fluxes of individual radio components was available, we have calculated the Browne 
parameter Rg which gives the ratio of core to extended emission at 5 GHz, and is 
believed to be correlated with source orientation (Orr and Browne 1982). 

Monochromatic luminosities (log(i/L„/erg/s)) in the visual (Ly at 5400A in the 
rest frame) and x-ray {L \ at 2 keV) are tabulated, as are the bolometric luminosities 
derived later in this paper. The absolute visual magnitude is also given, derived 
using the formulae of Veron-Cetty and Veron (19S7) but using our chosen cosmological 
parameters. q x is the x-ray spectral index derived from the IPC results; the ultraviolet 
to x-ray two-point spectral slope a 01 - is defined from 2500A to 2 keV in the rest frame. 
These spectral indices are conventional energy indices, a = — dlog .G/dlog v. 


4 Data 

4.1 Overview 

To assemble the SEDs for the 33 quasars required observations on 12 different 
telescopes, using 16 different instruments, in locations from ground-based to space. 
Because of this diversity of observing techniques it is necessary to describe each 
data set and the corrections made to it carefully. Table Nof their characteristics 
(wavelength range, A A/ A, typical signal to noise ratio). Also noted is the faintest 
detectable flux for each detector, and the brightest observed flux of an AGN in the 
band. These will be used in section §XX when discussing selection effects. Figure 
XX shows the bands covered graphically, and notes the type of detector, observing 
site, conventional names for the bands, and regions for which interstellar opacity is 
important. In the following sections each band is discussed in turn. 


4.2 X-ray 


The ’x-ray band’ covers two full decades of the spectrum. Even though most x- 
ray spectra have low resolution it is still useful to divide the band into three parts: 
‘hard x-ray’, covering 2-10 keV; ‘soft x-ray’, covering 0.3-2 keV; and ‘ultra-soft x-ray’ 
covering 0. 1-0.3 keV. All the 33objects in our sample were observed in the soft x-ray 
band in 1979 or 1980 with the IPC instrument on Einstein. The limited spectral 
resolution of this instrument allows us to characterize the 0.1 — 3.5 keV spectrum as 
a single power law modified by foreground absorption. For most of the objects in the 
sample, the data and spectral fits are presented in Wilkes and Elvis (1987). Although 
the foreground Galactic absorption is now well determined for our objects, we choose 
to retain the earlier fits in which the Galactic absorption is allowed to vary freely, 
since the presence of an ultra-soft excess in the incident spectrum is almost equivalent 
to a reduction in the absorption, and thus we get a better estimate of the slope in 
the soft band (beyond the excess) than if the absorption is fixed at its true value. 

Ultra-soft excesses were searched for in thirteen of the highest signal-to-noise 
observations. The data were re-analysed by Masnou et al (1991) using a two power 
law model and including data out to 10 keV from the MPC instrument, revealing 
the presence of an ultra-soft x-ray excess component in eight of the objects. This is 
indicated as a separate data point at 0.2 keV on our energy distribution plots (Fig. 
4). In these two power law fits, the foreground absorption was fixed at the known 
(Elvis, Lockman and Wilkes 1989) value. The object PGl211-fl43 was studied in 
a separate paper by Elvis et al (1991). These results are summarized in Table 3 
together with several previously unpublished spectral fits. The analysis for the new 
fits was identical to that of Wilkes and Elvis (1987). Contour plots for the new fits 
are presented in Fig. XX. 

Hard x-ray data are available for about half the sample from EXOSAT ME 
(Turner and Pounds 1989, Comastri et al 1991), and Ginga observations (Williams 
et al 1991) and the corresponding power law fits are listed in Table 3. The effects of 
foreground absorption are relatively unimportant in this band. 


4.3 Ultraviolet 

All of the 33 UV sample objects have been observed with both the long and short 
wavelength cameras on the International Ultraviolet Explorer between 1978 and 
1989. A total of 19 exposures were made explicitly for this program, many were 
long wavelength observations needed to fill in pre-existing short wavelength data. 
In addition, 108 further spectra were extracted from the Regional Data Analysis 
Facility archive. Both sets of data were analysed uniformly using the GEX Gaussian 
extraction algorithm (Urry and Reichert 1988) which is the most effective for our faint 
targets. Bad data points (reseau marks, cosmic ray hits and microphonic noise in the 



LWR) were removed interactively. We then averaged the ultraviolet flux within each 
of a set of wavelength bands chosen to avoid strong emission lines, converting the 
somewhat noisy spectra to a small set of relatively well determined continuum flux 
estimates. The continuum wavelength bands are defined to be fixedin the rest frame 
of the object (Table 4) and are 50A wide shortward of 1900A, 100A wide shortward 
of 5000A and 200A wide beyond that wavelength. Bands which overlapped a region 
of avoidance 8000 km/s to either side of strong broad lines or 2000 km/s to either 
side of strong narrow lines were omitted. We avoid the lines LaAl215, OIV A1402, 
CIV A1549, C7//JA1909, Mgl 1X2798, ATVA3426, [Oil] A3727, [0111] A5007/5009, 
JVe//A3869/3968, 01 A6300, and the Balmer lines Ha to HS. We note that weak 
lines and the blended [Fell] lines are not avoided and so these are included in our 
‘continuum’ fluxes. 

Table 5 lists the IUE exposures and the observed frequencies and fluxes corre- 
sponding to each continuum band. Note that the fluxes are given as log(j// t/ /Jy Hz) 
in the observed frame, and that they are not corrected for foreground extinction, 
(ljy Hz = 10~ 23 erg cm -2 s' 1 )- The first row on each page gives the rest frame 
wavelengths of each band. Then there is a header line for each object which gives the 
name of the object and the logarithm of the corresponding observed frequencies in Hz. 
This is followed bv two lines for each observation giving the date of the observation 
(in the first column) followed by the logarithmic fluxes (log(t/.F(i/)/l Jy Hz)) on the 
first line and the errors in the logarithmic fluxes on the second line. The etrors are 
1 a internal statistical errors derived from the scatter in the individual data points 
prior to binning. 

Table 5(a) contains data for the UV sample, while Table 5(b) contains data 
for other objects observed in the program. The optical spectrophotometry for the 
probable high redshift object PG 1407+265 (see discussion below) is also included 
in Table 5(a) for convenience, since the rest wavelengths for the observed optical 
spectrum lie in the ultraviolet. 


4.4 Optical 

Spectrophotometric observations for 14 objects were obtained with the Red (~ 
4500 - 7500 A) and Blue (~ 3200 - 6400A) Spectrographs or the FOGS (Faint Object 
Grism Spectrograph, ~ 4500 - 7500A) on the MMT. Blue spectrograph observations 
were made through a 5” circular aperture at air mass below 1.4 to minimize light lost 
due to atmospheric dispersion. Objects were then reobserved at highter resolution 
with a 1x3” aperture. A nearby standard star was observed immediately before or 
after the quasar observation. For both red spectrographs the large aperture used 
was 10x20” and the small one was a 1” long slit. In all cases the large aperture 
observations were used to flux calibrate the accompanying, higher S/N small aperture 
observations in order to obtain ~ 5 — 10.4 spectral resolution. The data were reduced 



in the standard manner, using IRAF. To ensure the photometric accuracy of these 
spectra, BVRI CCD photometric data were obtained on the FLWO 24-inch telescope 
within one week of the MMT observations. Table 7 gives the BVRI photometry 
values, estimated within a 14 arcsecond beam; Table 8 lists the continuum fluxes 
observed with FOGS and MMT spectrograph. As with the IUE data, emission lines 
have been avoided by averaging the logarithmic fluxes in the line-free, rest frame 
continuum bands listed in Table 4. 

Optical spectrophotometry for the 18 PG objects in our program were pre- 
sented by Neugebauer et al (1987, N87). These data were already corrected for 
Galactic reddening; this correction is removed here using the same law used by those 
authors for consistency with our database (Neugebauer, G., private communication). 
This allows us, and others, to apply a uniform Galactic dereddening to all of our 
optical, ultraviolet and soft x-ray data. Since the values of E(B — V) used in N87 
were not given in that paper, we tabulate them here (Table 6) for all the objects in 
that paper together with the values obtained using the prescription described below 
(Section 6.2). 

Optical data have also been included in the figures and analysis from Neuge- 
bauer et al (1979), Treves et al (1988), Sitko et al (1982), Condon et al (1981), Adam 
(1978,1985), and McAlary et al (1983). 


4.5 Near Infrared (1-3.5/um) 

Table 9 records measurements of JHKL (1.2 — 3.5 //m) photometry obtained at the 
MMT and the IRTF. The table gives the magnitudes and errors on the Johnson scale, 
the date of observation and the telescope and aperture used. 

We also included infrared photometry from Glass (1982), Rieke (1978), Condon 
et al (1981), Hyland and Allen (19S2), Sitko et al (1982), Rudy, LeVan and Rodriguez 
(1982), Ward et al (1987), and Neugebauer et al (1987). 

4.6 Far Infrared (> 10^m) 

Table 10 gives measurements at N and Q (10, 20 fim) made at IRTF and the United 
Kingdom Infrared Telescope (UKIRT). At the IRTF, the CT1 bolometer was used 
with a 6” beam and an east-west chopper throw of 30”. At UKIRT the UKT8 system 
was used with an 8” beam and a 20” east-west chopper throw. All magnitudes were 
derived from comparisons with standard stars (Tokunaga 1984, Elias et al 1982). 
Magnitudes are listed on the instrumental system, i.e. without color corrections. 

The Infrared Astronomical Satellite (IRAS, Neugebauer et al 1984) surveyed 
the sky in 1983 at far infrared wavelengths (12— 100//m). We have determined fluxes 



or upper limits at the positions of each of our sources. Where pointed Additional 
Observations (AO) were made, these were used; otherwise ‘lineadd’ (LA) estimates 
were made from the survey scans. This procedure allows a better estimate of the 
local background and errors in a particular measurement. The coadded survey maps 
were used to check for the presence of contaminating cirrus. The results are listed in 
Table 11; they are consistent (This is currently a lie!) with the results of Neugebauer 
et al (1986) and Sanders et al (1989) where we have objects in common (for the UV 
sample, 6 objects are previously unpublished). In some cases the 100/im upper limits 
are rather weak due to the presence of cirrus in the region of the source. Upper limits 
(3a) are listed for each source when no detection was made in any of the IRAS bands. 

4.7 Radio and Millimetre 

We have gathered core radio fluxes at 5 GHz for the UV sample from the literature 
(Table 12(a)). The PG sample study with the VLA by Kellerman et al (1989) was 
given preference over other references in the calculation of radio-loudness. By core we 
mean the flat spectrum compact component which appears to be physically distinct 
from the steep spectrum diffuse emission. Since this steep spectrum emission can 
itself be relatively compact in angular size (e.g. 3C4S), we do not use a fixed angular 
size to define the core, although in practice because detailed spectral information is 
not usually available we often use flux unresolved within a 1 arcsecond beam as our 
criterion. In the absence of spectral data this is an upper limit on the flat spectrum 
core flux. Beamsizes are given in the table. In Table 12(a) we also list measurements 
at other frequencies; and in Table 12(b) we list estimates of the flux of any extended 
radio source associated with the object. 

Millimetre wave data or upper limits are available for about half the sample 
(Table 13). Only III Zw 2 and 3C 273 are strong millimetre emitters, but four other 
sources have weak detections. 


5 Corrections 

5.1 Variability and averaging 

One limitation on our dataset is that the observations are typically not simultaneous, 
although the optical and IR data were generally obtained within about one month. 
For many of the objects we have observations at at two epochs (rarely more) in a given 
waveband, so we can make a crude estimate of the degree of variablilty. In table 14 
we list the observed range of variablity, F rnax jF m i n> and the associated timescale, at 
rest wavelengths in t lie near infrared, optical, mid ultraviolet and far ultraviolet. Fig. 



XXX. It can be seen that in the optical and infrared variability is not a serious problem 
for these ‘normal’ quasars, but that in the ultraviolet the variability is significant on 
timescales of a few years, although typically (13 out of 18 cases) it is less than a factor 
of two. These results indicate somewhat less variability found by Kinney et al (1991). 
We have to address this problem! 

To generate a single mean energy distribution for each quasar, we have taken 
an average (in logr'F(i')) of all the data in each frequency bin. However, for the 
IUE data, we have been selective in the exposures we chose to include in the average. 
Specifically, where simultaneous data from the long (LWP/LWR) and short (SWP) 
wavelength cameras were available, we have included them and excluded ‘orphan’ 
LWP/LWR or SWP exposures in order to avoid spurious steps in the data due to 
variability. Objects affected are Q0007+106, Ql 100+772, Ql 146-037, Q1202+281, 
and Q1613+65S. Further, when one exposure had significantly lower signal-to-noise 
than the rest available for a given object, it was omitted (Q1545+210, Q1613+658, 
and Ql 721 +343). 

5.2 Magnitude scales 

To include the optical and near infrared photometric data in the energy distributions, 
we have adopted for each band an effective wavelength and an absolute zero point 
appropriate for a flat energy distribution. We used the Hayes (1975) calibration of 
Vega, extended to other wavelengths by matching to a Kurucz (1979) theoretical 
model (9400K, log</ = 3.95). We convolved the Vega energy distribution with filter 
shapes from Tokunaga (1986) and Johnson (1965) to obtain absolute calibrations, 
assuming that Vega has a magnitude of +0.03 in all bands. However, we did not 
include the effects of atmospheric absorption in the calculation. The resulting absolute 
calibrations are listed in table 15. The color corrections are important for the optical 
bands but negligible in the infrared. 

5.3 Extinction corrections 

Corrections for foreground (Galactic) extinction are important in the optical, UV 
and soft x-ray. The x-ray spectral fits of table 3 already include correction for the 
line of sight absorption column as discussed above. A single extinction correction is 
applied to all the rest of the data, using an extinction law based on that of Savage 
and Mathis (1979) in the visible and ultraviolet, and Rieke and Lebofsky (1985) 
in the infrared beyond 3pm, Table 16. The magnitude of the correction has been 
estimated from the Galactic neutral hydrogen column by assuming a fixed conversion 
of N(HI)/E(B-V) — 5.0 x 10 21 cm 2 mag -1 (Burstein and Heiles 1978). The Galactic 
HI column (Table 1) has been accurately measured with a narrow beam and good stray 



radiation corrections in all but a few cases using the 140 ft Green Bank radio telescope 
(Elvis, Lockman and Wilkes 1989). The value for 3C 2<3 was taken from Dickey, 
Salpeter and Terzian (1978). In the three remaining cases it has been estimated from 
Heiles and Cleary (1979) or Stark et al (1984) 

5.4 Cosmological model 

We have adopted a standard Friedmann-Robertson- Walker cosmological model with 
D 0 = 2q 0 = 1 and H 0 = 50km s _1 Mpc _1 . After Galactic reddening corrections were 
applied, the data were blueshifted to the rest frame. Since in the rest frame we are 
working with the complete energy distributions, no k-corrections and no assumptions 
about the intrinsic spectrum are required. 


5.5 Host Galaxies 

Although the overall energy output of our sample objects is dominated by the active 
nucleus, in the near infrared and optical the host galaxy can make a significant 
contribution. In a few cases (Q0049+171, Q1426+015, Q1501 + 106, Q1613+658) 
where a strong host galaxy is present in optical images and large aperture photometric 
data shows a significant excess over the small aperture spectrophotometry, we omit 
the photometric data from the final energy distributions. 

Table 17 gives absolute magnitudes of the host galaxies and central point 
sources as found in the literature for 17 of the UV sample quasars. In cases where 
the change in the near infrared energy distribution is more than 5 per cent, we have 
subtracted a host galaxy template normalized to the quoted host galaxy magnitude. 
Our host galaxy template is based on the Sab galaxy model of Pence (1976). Fig. XX 
shows the host galaxy subtracted energy distributions. 


5.6 PG 1407+265 

The redshift of PG 1407+265 (= 2E 3196) is uncertain; the object has very weak 
emission lines, the only certain feature coming at 5500A. An identification of this as 
Mg II and the weak presence of GUI] led Schmidt and Green (1983) to propose the 
redshift as z=0.944. The absence ol a clear Lyman alpha line in our IUE spect a 
caused us to consider lower redshift identifications, including blueshifts, but no choice 
of redshift allows normal quasar line ratios. We have adopted the high redshift given 
by Schmidt and Green because of the position of a prominent continuum feature, 
namely the inflection at the beginning of the ultraviolet bump. The rest wavelength 
of this feature lies between 1.0 and 1.5 microns for all our other objects, and in this 
object occurs at an observed wavelength between 2.0 and 2.4 microns, which lends 



support to the high redshift estimate. A weak broad line is marginally detected in 
the LWP spectrum. If this line is Lyman alpha, it peaks at approximately z=0.99, a 
velocity shift of 15000 km/s with respect to the Mg II line. 

We continue to adopt the Schmidt and Green value despite the absence of 
definitely observed hydrogen lines. We note that the broad width of the line, the 
overall continuum shape, the lack of strong optical variability and the x-ray to optical 
flux ratio indicate that the object is indeed a quasar rather than, for instance, a star, a 
BL Lac object, or some more exotic object. The near infrared to ultraviolet spectrum 
of PG1407+265 is presented in Fig. 3. We believe that the sharp drop in flux 
level between the LWP and SWP spectra at 2000 ± 100A is likely due to ultraviolet 
variability rather than a Lyman absorption system (which would be at z ~ 1.1 , 
inconsistent with the adopted redshift). Note the good agreement in flux level at the 
overlap between the two optical spectra and the smooth continuation of the spectrum 
made by the infrared photometry, implying a lack of variability. 


6 Properties of the energy distributions 


6.1 Luminosities: Bolometric and individual bands 

In Fig. 4 we present two energy distributions for each object, an overall view and a 
closeup of the infrared to ultraviolet region. The overall view covers a fixed flux range 
of ten decades and illustrates the radio-loudness and x-ray properties of each quasar. 
The closeup view covers two decades in flux and illustrates the l//m inflection, and 
the strength of the blue bump. 

To characterize the large scale distribution of the energy output of the quasars, 
we calculate integral luminosites in a set of broad bands. The integrals are calculated 
by running a simple linear interpolation through the data points in log vL u space, i.e. 
connecting the individual points with a power law. The errors indicated below are 
estimated by performing the integrals using the one sigma high and one sigma low 
flux values instead of the nominal values. For upper limits we interpolate between 
detections on either side. The lower of the interpolated value and the upper limit 
is used as the nominal flux estimate, but the errors are estimated using zero as the 
lower error bar and the upper limit as the upper error bar. 

1. Bolometric The bolometric luminosity is typically well defined except for two 
regions of the spectrum: the mostly unobservable EUV region and the as yet 
unobserved hard x-ray and gamma-ray region. Indications are strong that the 



0.1-10 mm gap is energetically negligible. As a first crude estimate of the EUV 
luminosity, we simply make a linear interpolation between the ends of the IUE 
and Einstein ranges. A reasonable upper limit to the EUV luminosity can 
be made by finding the maximum blackbody curve which does not exceed the 
observed data; however this limit is not strong, as the luminosity implied is 
typically 10 to 100 times the luminosity observed in the rest of the spectrum 
(Fig. 5.) 

We perforce neglect the unknown luminosity above 10 keV (10 18-4 Hz). 

2 . UVOIR The three decades between lOO^rn and 0.1/im (the ultraviolet/ optical/- 
infrared or ‘UVOIR’ region) are relatively well sampled, so the corresponding 
UVOIR luminosity can be much more accurately calculated. This luminosity 
accounts for most of that which is directly observed, and so is a useful number 
for quoting as ‘the luminosity’ of a quasar. 

3. Decades We also tabulate the luminosity in individual decades across the 
electromagnetic spectrum. Outside the range 1 — 0.1/zm, these are often 
estimates from only one or two points and the errors are correspondingly large 
(typically 25 percent in the far IR). 

4. Octaves In order to describe the shape of the ultraviolet bump component, we 
define a set of narrower octave wide bands. The four bands we call IR (1 — 2/im), 
VIS (4000 - 8000 / 1 ), NUV (2000 - 4000 A ) and UV (1000 - 2000A). The IR 
band gives the luminosity longward of the 1 pm inflection, where there should 
be little contribution from the bump. The UV band measures the luminosity 
in the bluest observed part of the bump, while the VIS band samples the early 
part of the bump’s rise. The NUV band covers the near ultraviolet region 
of the spectrum dominated by the ‘small bump’ of blended Fe II and Balmer 
continuum emission (Wills, Netzer and Wills 1985). 

The calculated integral luminosities are tabulated in Tables 18 to 20. 

6.2 The Mean Energy Distribution and its Dispersion 

We use the interpolations obtained in the previous sections to obtain mean energy 
distributions for the radio-loud and radio-quiet quasars in the sample. Because we 
believe the 1 p.m inflection to be an important feature separating physically different 
continuum components, we normalize each energy distribution at that wavelength by 
dividing by T Bose . Fig. 8(a) shows the two mean normalized energy distributions. 
The difference in x-ray slopes between the radio-loud and radio-quiet quasars is 
immediately apparent, but the most striking feature is the close agreement in the 
UVOIR region, as also found by Sanders et al (1989). This is true despite the wide 


dispersion in shapes of the individual objects. Because of this agreement we feel it is 
useful to provide an analytic approximation to our mean quasar energy distribution. 
Fig. 8(b) shows a closeup of the overall mean distribution in the UVOIR region 
(solid line), together with the extreme envelope found in the sample (dotted lines). 
Note that the ‘small bump’ (A ~ 3500A) is present in the data. Table 22 contains 
approximate piecewise power law fits to the curves of Figure 8, together with the 
measured dispersion for each segment. 

We have also derived the mean bolometric corrections relative both to our 
near-infrared baseline normalization and the more usual normalization at V. The 
mean value of Ly/Lsase = 1.4 can be used to convert between the two normalizations. 
We derive an Lg^/Ly and Luvoir/ L v of 12.7 ± 4.2 and 8.5 ± 2.9 respectively. The 
mean ultraviolet to infrared color is 3.1 ± 1.3. The details of these and further mean 
colors are given in Table 23. 

Check Sanders et al 10 12 Lsw?x separation. Redo dispersion rel to UVOIR. Check 
dispersion is same for RL and RQ. How many objects not in common? 


6.3 The Spectral Window Function 

While an average spectrum gives a reasonable representation of the quasar continuum, 
the observed ranges of continuum shapes is quite broad. The observed diversity of the 
quasars in our sample is a factor of 10 in the UV and at 100/xm. This is as extreme as 
could be observed as a result of what we shall call the ‘Spectral Window Function’. 

Bounding the spectra in figure NN are two solid lines. The lower line is the 
limiting sensitivity of the telescopes/instruments used in each frequency range. The 
upper line is the upper envelope of the brightest objects observed at each frequency. 
This also is a limit to the range of observable continuum shapes. In effect this is a limit 
imposed by the space density of AGN; were they more common there would likely be 
one closer to us and brighter. Only AGN lying between these two lines are observable. 
The region between the two curves we call the ‘Spectral Window Function’, by analogy 
with the window function familiar in timing studies. Where this function is narrow all 
AGN will necessarily have spectra that are similar. Wherever the window function 
dominates we are unlikely to be seeing the whole range of spectral shapes in the 
quasar population. A detailed statistical treatment for the effects of the Window 
Function needs to be developed. 

The window function is particularly narrow in the far-infrared and in the X- 
rays. For the far-infrared the IRAS 12, 25, 60 and 100/xm band sensitivities determine 
the width of the function to be only a factor of ~10. Table NNN lists the width of 
the window function and the extreme of the dispersion of quasar shapes (normalized 
to Luvoir) in our sample. 


The similarity of the far-IR spectra of detected quasars seems likely to be a 
reflection of the narrow accessible range. IRAS upper limits for the ‘Weak Infrared 1 
quasars encourage this belief, since some seem to be IR-quiet ( e.g . PG0026+129, see 
below). One third of the PG quasar sample was not detected with IRAS. Sanders et al. 
argue reasonably that these are mostly the same as the detected quasars. However 
some IR-quiet quasars may exist, as discussed below. 


7 Discussion 

7.1 The 1/im Inflection 

The inflection around 1/im noted by Neugebauer et al (1979) is present in 28 of our 
objects, while in the remaining 5 it is difficult to rule out its presence because of 
the lack of IR and red observations. Note that in the ultraviolet weak objects, an 
inflection is still seen. 

The lowest point, which we will call the baseline monochromatic luminosity (cf. 
Carleton et al 19S7), is almost always coincident with the photometric H, or sometimes 

J, band measurement. While the existence of the inflection is independent of the H 
and I\ measurements, as in over hall the objects the mid IR and the optical have slopes 
of opposite sign, the precise wavelength and flux of the minimum are dependent on 
the H and K measurements and therefore also on the adopted calibration at H and 

K. We have checked our results using a number of other groups 1 calibrations and 
find that the existence of the inflection remains clear and is not an artefact of the 
conversion to absolute fluxes. 

The location of the inflection is relatively constant, as noted by Elvis et al 
(1986). The inflection wavelength always occurs between 0.8 and 1.5 microns, and 
usually between 1.0 and 1.5 microns, in agreement with Neugebauer et al (1987). The 
inflection wavelengths and corresponding baseline monochromatic luminosities L^ase 
using our magnitude scale zero points are listed in Table 21. 


7.2 Strength and Shape of the Infrared Continuum 

The infrared continuum consists in general of a broad bump, with the flux rising 
from 1/im towards longer wavelengths, peaking and falling by many decades, with 
very little emission in the submillimeter. Submillimeter observations of our objects 
are lacking; only four objects have IRAS and IRAM observations which constrain the 
slope of the mm-break to d\ogS„/d\og v > 2. 

Within this overall shape there is substantial variety. We identify three kinds of 
infrared continuum behaviour: 1 normal’ objects in which the infrared continuum peak 



occurs between 5 and 20 microns, and the flux falls off strongly at longer wavelengths; 
‘Strong Far Infrared’ objects with a more complex continuum in which the flux rises 
again after the mid infrared peak’ and ‘Weak Infrared’ objects where the spectrum is 
flat or falls off at near infrared wavelengths. 

In the (60 /xm/12 /im) versus (1.2 /xm/0.36 /zm) color-color classification 
diagram of Ward et al (1987), all of our objects are ‘Type A’, lying in the optically 
blue, far-IR weak quadrant corresponding to their ‘bare AGN’ classification, with the 
exception of the two far-IR bright objects, 3C 48 and Mkn 876, which lie in their 
‘empty’ quadrant (Fig XX.) 

Defining Ljr to be the sum of the far and near IR decades of Table 18, 
and LiRBa.se = (21n 10)Z/Sase which converts the monochromatic luminosity to an 
equivalent integrated luminosity over two decades, we construct the color-color 
diagram of Fig. 7, which replaces the optical/infrared color of the Ward et al diagram 
with a purely infrared criterion. 

Of those 26 of 33 objects in the sample with mid and far infrared detections or 
useful upper limits, 15 are classified as ‘normal’, with a single infrared peak between 
3 and 20 /mi. Five have abnormally weak infrared emission, and a further five have 
a more complex ‘strong far IR’ spectrum. 

• Normal Infrared Continuum 

It is apparent from Fig. 4 that in many of our objects (e.g. Mkn 586) the 
turnover in the far infrared occurs well shortward of 100/xm. The ‘3 micron 
bump’ seen in 3C 273 peaks at longer wavelengths in most objects, between 5 
and 20 microns, and is the principal infrared continuum peak. This range in 
turnovers suggests a peak defined by a temperature or a size rather than some 
atomic process. 

• Strong Far Infrared 

These objects have log(L(60/nn)/Z(12/xm) > -0.1 because of the extra peak 
in their spectrum. In only two of our objects (Mkn 876 and 3C 48) is the 
flux strongly rising out to 100/rni; this may be due to the presence of a large 
starburst in the host galaxy, at least in the former object (Yeeand Green 1987). 
In the milder examples PHL 909 and PKS 2135-147 the flux rises again slightly 
toward longer wavelengths; these may represent an intermediate case. These 
four objects (McDowell et al, 1989) are also notable for their weak ultraviolet 
emission as can be seen in the IR/UV color-color diagram (Fig. 6). 3C 273 is 
also in the strong far IR group because of its millimetre-loud blazar component 

• Weak Infrared 

The weak spectra have a low value of L(l-100/xm) relative to the infrared 
baseline flux of Table 20, with L, n < Li RBa se (Fig. 7). They are PG0844+349, 



PG 1 116+215, 3C249.1, and the extreme examples Kaz 102 and PG0026+129. 
In the latter two objects, a small 3pm peak falls away rapidly in the mid infrared. 


7.3 The Strength and Shape of the Ultraviolet Bump 

In Fig. 9 we use the IR, V and UV integral octave luminosities to construct a color- 
color diagram illustrating the variation of spectral shape of the bump in our sample. 
Most of the objects lie in a relatively well defined range of IR to UV color (overall 
bump strength), but show a wide range of shapes. The width of the distribution is 
comparable to the width induced by ultraviolet time variability. We note that the 
locus occupied by the objects begins near the power law line (dotted) and extends to 
highly concave spectra which are flat in the red but very steeply rising in the blue 
(i.e. log L{IR)/L{V) ~ 0 but log L(V)/L(UV) ~ -0.6). The solid line represents a 
pure face on Kerr accretion disk model, marked with logarithmic luminosity values 
appropriate for Eddington-limited accretion, and the dashed line repesents mixing an 
accretion disk with LboI = 10‘ ,< erg/s with a power law extending from the infiared 
of unit slope in L„. No objects lie outside the region defined by mixing these two 
components. Such combined models fit the colors well; however, we leave detailed 
modelling to a future paper (Kuhn et al 1991, in preparation). 

A smaller number of objects are significantly redder than the mean. These 
four ‘weak bump’ objects (McDowell et al 1989) have normal x-ray to infrared ratios, 
but their ultraviolet to infrared ratios are unusually low. Judging from the presence 
of XUV excess emission in some of the objects, it is unlikely that the weakness of the 
bump is due to internal extinction. There is no correlation of bump strength with 
luminosity in this sample. The ‘weak bump 1 objects are also moderately, but not 
exceptionally, weak in their ultraviolet to x-ray ratio. 


7.4 Luminosity correlations 

We reproduce the results of Kriss (1988). 


7.5 The Infrared to X-ray Connection 

Following Carleton et al (1987), we studied the connection of the infrared component 
to the x-ray. Carleton et al found a tighter correlation between the lowest infrared 
point (the “infrared baseline”) and the hard x-ray flux, for their hard x-ray selected 
sample. For our sample, the scatter in the correlation of the IR baseline with the 
x-ray was about the same as for the other correlations, in agreement with Sanders 
et al (1989) and in contrast to the Carleton et al results where a tighter correlation 
was found. We note that we used 1 keV fluxes rather than 6 keV, and that the mean 


luminosity of our sample is a factor of 30 higher. We also note that Brissenden (1989) 
does find an effect for low luminosity objects found by HEAO-1, suggesting that the 
effect may hold preferentially for x-ray selected objects. 


8 Conclusions 

We have collected a set of data useful for many investigations. In particular, these data 
will serve as a benchmark against which to compare the energy distributions of high 
redshift quasars. The data are available by anonymous ftp from cfa248.harvard.edu. 

We have studied the properties of a large sample of complete quasar energy 
distributions. We find that 

1. The bolometric luminosity of the quasars is a factor of 12.7 ± 4.2 times the 
monochromatic visual luminosity, using a simple interpolation in the EUV; the 
more accurately defined ‘UVOIR luminosity’ is a factor of 8.5 ± 2.9 times the 
monochromatic visual luminosity; 

2. The IR emission can be weak relative to the optical; 

3. “The submillimeter break” often occurs at wavelengths of 60^im or shorter; 

4. The strength of the ‘ultraviolet bump’ varies by a factor of 10, but typical objects 
have an ultraviolet luminosity three times their near infrared luminosity; 

5. The shape of the ultraviolet bump varies significantly from object to object, 
rising steeply in the optical for some objects, but remaining flat there for others; 

6. Objects with weak ultraviolet bumps are luminous in the far infrared; 

7. We confirm the near constancy in wavelength of, and ubiquity of, the “1 /zm 
inflection” first noticed by Neugebauer et al (1979) and studied by Elvis et al 
(1986) and Neugebauer et al (1987). 

8. The ultraviolet variability of these objects (chosen to exclude strongly optically 
variable quasars) is typically less than a factor of 2 on a timescale of years, and 
many vary by less than twenty per cent. 

9. While we confirm the correlations found by Kriss (1988), we suspect that they 
may be induced by host galaxy contamination at low luminosities. 

10. The ‘infrared to x-ray connection’ discussed by Carleton et al (1987) is not 
present in this sample. 
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Figure captions 


Fig. 1. Examples of radio-loud (3C 273, top) and radio-quiet (Mkn 586, bottom) 
quasar energy distributions, illustrating the main continuum features. The energy 
distributions show the logarithm of the energy per unit logarithmic frequency interval, 
in the rest frame. 


Fig. 2. Histograms of sample properties. Shaded bins correspond to radio-loud 
objects, R l > 1. The absolute magnitude is calculated using the corrections of Veron 
and Veron (1987), but with our value of 0 0 - The bolometric luminosity is derived 
from the observed energy distributions as described in the text. The monochromatic 
x-ray luminosity is vL(i / ) at 2 keV in the rest frame, estimated from the IPC spectral 
fits. R b , the ratio of radio core to extended luminosity (Orr and Browne 1982), is an 
indicator of the source orientation; it was only possible to estimate this for a subset 
of the sources. 

Fig 3. The observed energy distribution of PG 1407+265 from the near infrared to 
the far ultraviolet. The expected positions of prominent quasar emission lines are 
indicated for an assumed redshift 2 = 0.94. The short and long wavelength IUE 
observations are not contemporaneous and the low level of the short wavelength data 
may be due to variability. 

Fig 4 (a-ag). Rest frame, dereddened continuum energy distributions of the quasar 
sample. For each object, the left hand panel shows the overall radio to x-ray energy 
distribution, and the right hand panel shows the details of the UVOIR (100/zm to 
1000A) region. 

Fig 5. The maximum black body in the EUV region, compared with the power law 
interpolation, for the quasar Mkn 586. This shows that the total energy in the EUV 
is not well constrained by the observations. 

Fig 6. Infrared-ultraviolet color-color diagram, plotting the far IR color 
L(60fim) / L(\2fim) against the UV bump strength L( 1 — 2fim)/ L(0.l — 0.2/rm). The 
objects with weak ultraviolet bumps are also luminous in the far infared. 


Fig 7. Infrared color-color diagram, plotting the far IR color Z,(60/rm)/T(12/xm) 
against the IR bump strength 1(1 - 100/im)/I Baje . The dashed lines indicate the 
regions corresponding to the different spectral classifications. 



Fig 8. (a) Mean energy distributions for radio-loud and radio-quiet quasars. 

(b). Overall mean energy distribution in UVOIR region. Dashed lines indicate the 
extreme departures observed from the mean distribution. 


Fig 9. Ultraviolet-optical-near infrared color-color diagram, indicating the shape of 
the ultraviolet bump. 



EINSTEIN UV SAMPLE QUASARS 



References to NH: (1) Elvis, Lockman and Wilkes 1989; (2) Heiles and Cleary 1979; (3) Dickey, Salpeter, and Terzian 1978; (4) Stark et al, 1991 
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TABLE 5(a) (continued 
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...... 0.11 0.10 0.08 0.13 0.09 0.08 

1982 Dec LWP1739 SWP18741 log(i/ 0 P v „) ------ 1198 1172 1175 1183 1169 1183 

0.14 0.13 0.13 0.18 0.11 0.07 
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ULTRAVIOLET CONTINUUM FLUXES: UV sample quasars 
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Q0134+329 1 Oft (Vo) 15.12 15.10 15.09 15.08 15.04 15.02 15.00 14.98 14.96 14.94 14.93 14.91 14.87 14.86 14.84 

1982 Oct log(i/ o F„ 0 ) 11.79 11.87 11,81 11.71 11.72 11.78 108 IT08 - I I I I I ? 

0.23 0.25 0.26 0.23 0.19 0.18 0.28 0.39 














TABLE 5(a) (continued 
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TABLE 5(a) (continued 



Q1211 + 143 log(i'o) 15.22 15.20 15.19 15.18 15.14 15.12 15.10 15.08 15.06 15.04 15.03 15.01 14.97 14,96 14.94 

1985 Dec log(«/ o /v 0 ) 12.78 12.73 12.69 12.69 12.71 12.73 12.69 12.75 1277 12/79 12/76 12/76 1278 I 7 

0.04 0.03 0.05 0.08 0.12 0.11 0.06 0.05 0.04 0.03 0.03 0.03 0.06 

1987 Feb logfi/oFVJ 12.65 12.58 12.55 12.58 12.60 12.67 12.64 12.71 12.72 12.72 12.70 12.71 12.72 

0.06 0.03 0.04 0.09 0.12 0.10 0.08 0.05 0.04 0.03 0.04 0.03 0.06 



TABLE 5(a) (continued 



Q1307+085 log(i/ 0 ) 15.19 15.18 15.17 15.15 15.11 15.09 15.07 15.05 15.03 15.02 15.00 14.99 14.94 14.93 14.92 

1980 May log(o 0 Fi, 0 ) 12.36 

0.04 

1984 Feb log(l > 0 F Vo ) 12.25 12.21 12.19 12.17 12.10 12.14 12.12 12.13 12.13 12.16 12.16 12.16 

0.01 0.01 0.00 0.03 0.07 0.02 0.03 0.01 0.01 0-01 0-02 0.01 ; ; ; 







TABLE 5(a) (continued 



Ql 61 3+658 logt^o) 15.20 15.19 15.18 15.16 15.12 15.10 15.08 15.06 15.04 15.03 15.01 15.00 14.95 14.94 14.93 

1981 Sep log(i/ o F„ 0 ) 12.19 1 2.12 1 2.02 1 2.02 1 2709 TTH I I I I I I 1 1 7 











TABLE 5(a) (continued 




0.13 0.08 0.07 0.12 0.13 0.11 0.07 0.10 0.07 0.09 

1982 Dec log(i-of„ 0 ) 11.96 11.88 11.96 11.85 11.83 11.83 11.86 11.89 11.96 12.00 

n 17 0 15 0.18 0.15 0.13 0.10 0.07 0.05 0.04 0.04 
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fABLE 5(b) (continued 




ULTRAVIOLET CONTINUUM FLUXES: Other IPC quasars 



Q1635+1 19 log^o) 15.19 15.18 15.17 15.16 15.12 15.10 15.08 15.06 15.04 15.02 15.00 14.99 14.95 14.93 14.92 

1988 Apr log(i/ o F„ 0 ) 11.70 11.40 . "™" . 

0.36 0.41 - 















TABLE 5(b) (continued 
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OPTICAL CCD PHOTOMETRY: UV sample Quasars 
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OPTICAL CCD PHOTOMETRY: Other IPC Quasars 
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OPTICAL CONTINUUM FLUXES: UV sample quasars 
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NEAR INFRARED PHOTOMETRY: Other IPC quasars 
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N AND Q BAND INFRARED PHOTOMETRY: UV sample quasars 
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IRAS FLUXES: UV sample quasars 
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IRAS FLUXES: Other IPC quasars 
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RADIO CORE FLUXES 
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Abstract 


We have measured integrated Galactic 21cm column densities toward 200 
quasars and active galactic nuclei using the NRAO 140 ft telescope at Green Bank. 
These data have been corrected for stray radiation with the technique of Lockman 
et a/.(1986). The 21 arcmin beam size of the 140 ft is small enough to minimize 
the uncertainty in N H due to angular variations in the HI of the Galaxy at high 
latitudes. The resulting column densities are accurate to 1 x 10 19 atoms cm 2 or 
^ 5%, whichever error is larger. Opacity uncertainties dominate the errors above 
N// ~ 4 X 10 20 atoms cm 2 . 



I. The Need for Accurate N// Values 

The extension of sensitive extragalactic x-ray measurements into the band 
below the 0.28 keV Carbon edge ( the ‘C-band’, McCammon et a/.1983) has made 
it important to determine accurately the opacity of the interstellar medium of our 
own Galaxy at these energies. Interpretation of spectral results from the Einstein 
IPC (Imaging Proportional Counter, Gorenstein et a/.1981) and the EXOSAT CMA 
(Channel Multiplier Array, Taylor et a/.1981) depend critically on a knowledge of the 
Galactic N //. The inferred properties of the recently discovered ‘ultra-soft excesses’ 
in the x-ray spectra of quasars and active galactic nuclei (AGN) are particularly 
sensitive to this value (Wilkes and Elvis 1987, Branduardi et al. 1987, Giommi and 
Tagliaferri 1987). 

Most x-ray astronomers currently use the unpublished but widely circulated 
‘Bell Labs’ survey (Stark et al., 1984) of 21cm HI emission to determine the 
column densities of interstellar material toward extragalactic sources. This survey 
is superior to previous all-sky HI surveys because the telescope’s unblocked aperture 
makes it almost free of ‘stray radiation’ entering the signal from distant sidelobes. 
The small size of the antenna, however, has the disadvantage that its angular 
resolution is poor, typically ~ 2° x 3°. Small scale structure in the HI distribution 
on the sky is thus liable to introduce errors in the column toward any particular 
line of sight. Elvis et a/.(1986, Appendix B) estimated that the 90% errors 
introduced by small scale structure were ±1 x 10 20 atoms cm -2 . For a typical 
high Galactic latitude position with N//~ 3 x 10 20 atoms cm -2 this uncertainty of 
1 x 10 20 atoms cm -2 in N# leads to an uncertainty in the intrinsic x-ray flux density 
of a source at 0.2 keV of a factor 2.4. More accurate column densities than the Bell 
Labs survey provides are thus clearly needed. 

A method for achieving high angular resolution HI measurements free 
from stray radiation has recently been developed for the Green Bank 140 ft 
telescope (Lockman, Jahoda and McCammon 1986, Appendix A). We have used 
this technique to produce accurate HI column densities toward a large number of 
quasars and AGN that are bright and well-observed x-ray sources. We present these 
column densities here. We also analyze the results to provide an assessment of the 
remaining errors and to determine more carefully the uncertainties involved in using 
the Bell survey values when no more accurate measurement is available. 


II. Observations 

The sample of objects includes a wide selection of x-ray observed AGN and 
quasars. It includes, so far as they were observable from Green Bank (t.e. declination 
> -40°): all the quasars with Einstein IPC spectra in Wilkes and Elvis (1987); all 
the BL Lac objects with Einstein IPC spectra in Madejski (1985) ; a substantial 


number of other Einstein observed quasars for which IPC x-ray ‘colors’ are available 
(Brunner et al., in preparation); all of the members of the complete hard (2-10 keV) 
x-ray selected ‘Piccinotti’ sample of AGN (Piccinotti et a/.1982); and a selection 
of the AGN detected by the Medium Energy instrument on EXOSAT ( Sternberg 
et a/.1986). In addition we mapped the accessible Einstein Deep Survey regions 
(Giacconi et a/.1979), the ‘Braccesi’ survey region covered by Marshall et a/.(1984), 
and SA57 and the North Ecliptic Pole regions that will be observed in depth by 
ROSAT. The maps for these extended regions will be presented separately. 

The observing technique was as detailed by Lockman et a/.(l986, Appendix 
A). Briefly, the 140ft telescope is used to map out the beam covered by the Bell 
Labs observation that was centered closest to the position of each AGN. The 140ft 
data are then convolved with the beam profile of the Bell antenna to give the 21cm 
spectrum that the 140ft would observe if it had the same beam shape as the Bell 
antenna, i.e.including the stray radiation seen by the 140ft at this position. This 
spectrum is dependent on the altitude and azimuth of the telescope as well as on 
the celestial position. The difference between the Bell spectrum and the convolved 
140ft spectrum is then the stray radiation seen by the 140ft at that position. This 
is then be subtracted from the 140ft spectrum taken at the precise position of the 
active galactic nucleus. 

In practice two independent spectra were taken at each AGN position. Total 
integration times were 2-6 minutes. Because of the low system temperature (18 
K at zenith), these integration times are sufficient to make the uncertainty in N H 
due to noise negligible. An observation of particularly low signal- to- noise is shown 
in figure 1. Each AGN position was observed both immediately before and after 
the Bell beam was mapped and the results averaged, to minimize time dependent 
effects. Remaining errors in the observation are primarily due to uncertainties 
in the baseline (Lockman et a/.1986) and amount to ~ 2 x 10 18 atoms cm -2 . 
The amount of the stray radiation removed was equivalent to HI columns of a 
fewxlO 19 atoms cm -2 , or ~ 15% for typical values of N# . Figure 2 shows an 
example of a particularly large stray radiation spectrum compared with the observed 
140ft spectrum. There is uncertainty in the stray radiation correction at a level of a 
few times 10 18 atoms cm -2 and the form of the uncertainty is not well determined. 

To derive an integral N h measurement from the observed spectrum requires 
a correction for the opacity, r, of the gas. Although at high galactic latitudes this 
correction is generally s ill it is >metimes dominant. The correction is intrinsically 
uncertain. The observer spectrum is made up of contributions from different clouds 
in the interstellar medium which generally have a wide range of temperatures and 
optical depths and many possible geometric arrangements (see discussions in Dickey 
and Benson 1982 and Lockman and Dickey 1989). The distribution of these clouds 
cannot be derived from a 21cm spectrum. The size of the uncertainty can be 



estimated using the peak observed brightness temperature in the spectrum, T pea k , 
since 

Tpeak = T,pi „(1 - e- T ) 

if we ass um e that the emission arises from gas of a fixed kinetic temperature (called 
the spin temperature, T, p i„, for historical reasons). If T pea k is small then r is always 
small for reasonable values of T spin , but for larger T peo * t is more likely to be 
significant and one must make an opacity correction. Figure 3 shows the uncertainty 
in N h due to opacity corrections, AiV//(T jptn ), plotted against T pea k- A./V//(T,p, n ) 
was calculated using T, pin values of 10 4 K and 125 K as extreme cases. A best fit 
to figure 3 gives 


A N H (T apin ) = (0.073 ± 0.005)T peo * - (0.063 ± 0.007) (10 20 atoms cm 2 ) 

This expression, or figure 3, can be used to estimate the opacity correction 
uncertainty in our N // measurements. 

The best estimate column densities are given in Table 1. Values of N// 
are given for an assumed effective spin temperature of 250 K which is a good 
approximation to that of the interstellar HI (Dickey 1988, private communication, 
see also Dickey 1988). When the HI signal is brighter than several tens of Kelvins 
the Bell Labs survey data, and hence our techniques for removing stray radiation, 
become unreliable. But the brightness of the HI means also that stray radiation 
will be a relatively unimportant component of the total N# and the “uncorrected" 
140 ft spectra are thus fairly accurate. Sources so affected are flagged in table 1; 
the uncertainty in their N w is dominated by the uncertain opacity correction. We 
also tabulate the peak brightness temperature (T pea fc). The dominant error in the 
tabulated column densities may be due to either the baseline noise referred to 
above or to the uncertainty due to the opacity correction. Figure 4 shows how 
the uncertainty in N H due to T, pin ,AiV//(T jpin ) dominates above the measurement 
errors of 1 x 10 19 atoms cm -2 (dashed line) for N// > 4 x 10 2 ° atoms cm -2 . If 
the area under Galactic HI profiles could always be measured to ~2%, then opacity 
correction uncertainties would dominate at all values of N//. 

Figures 5 a,b compare the optically thin 140ft and Bell values of N// as both 
a ratio and a difference as a function of N// (140ft). The ratio of 140ft/Bell stays 
constant with N h while the difference seems to grow slightly. This suggests that 
small scale structure is a constant fraction of the total. That is, at high latitudes 
the galactic N H in a 21’ field is equal to that in the 3° x 2° Bell fields with a 1 <t 
uncertainty of 11% on average. 


IV. Conclusions 


The values of N // tabulated here have typical uncertainties of 
~ 10 19 atoms cm -2 and so axe about ten times more accurate than column densities 
derived from the Bell Labs survey alone. It is worth noting here that for the purposes 
of investigating small amounts of x-ray absorption it is the column density due to 
hydrogen and helium alone that contributes to the photoelectric cross-section in 
the ‘C-band’ below 0.28 keV (Morrison and McCammon 1983). Hydrogen in any 
form other than atomic appears to be uncommon at high Galactic latitudes (Blitz, 
Magnani and Mundy 1984) so that the 21cm data give an unusually direct and clean 
means of determining the column that must absorb x-rays. 

The extreme x-ray-ultraviolet telescopes on ROSAT (Pye 1984) and EUVE 
(Bowyer 1987) can only detect extragalactic sources if they lie in directions of 
unusually small Galactic column density. There are two AGN in table 1 with 
N// below 0.9 x 10 2 ° atoms cm -2 and seven with N h below 1.1 x 10 2 ° atoms cm -2 . 
These low N# AGN will make good targets for the extreme x-ray-ultraviolet 
telescopes. Our nearly 200 AGN form an essentially random sample of the high 
latitude sky. We can estimate then that roughly 1% of the sky has column densities 
of 0.9 x 10 2 ° atoms cm -2 or less. 

The need for more accurate column densities will rise rapidly with the many 
soft x-ray sources that are expected to be discovered with the soft x-ray sky survey 
satellite ROSAT (Triimper 1984, due for launch in February 1990). The effects of 
uncertain Galactic N h can be important to logN-logS studies and hence to questions 
of quasar evolution and the origin of the x-ray background (Zamorani et al. 1988). 
An all-sky survey of Galactic N// with a similar beam size to that used here and 
free of stray radiation would substantially reduce these problems and would be of 
great utility for x-ray astronomy. 
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Table 1: Galactic Column Densities Toward AGN 
From Green Bank 140ft Data 


Object 

R.A. 

decl. 

m. 

Tpcak 

MCS 18 

00 02 46.2 

05 07 30 

3.87 

9.6 

PHL658 

00 03 25.0 

15 53 08 

3.94 

14.7 

IIIZW2 

00 07 56.7 

10 41 47 

6.09 

24.5 

3C 9 

00 17 50.0 

15 24 17 

4.02 

18.7 

PG0026+129 

00 26 38.1 

12 59 29 

4.93 

18.6 

0038-020 

00 38 23.7 

-02 02 51 

2.57 

5.9 

UM 275 

00 43 39.5 

00 48 04 

2.03 

5.6 

PKS0044+030 

00 44 31.3 

03 03 34 

3.00 

9.0 

PG0049+171 

00 49 16.5 

17 09 40 

4.26 

16.5 

I ZW I 

00 50 57.7 

12 25 20 

5.07 

20.3 

0052+251 

00 52 11.1 

25 09 24 

4.50 

13.9 

PHL909 

00 54 31.9 

14 29 59 

4.20 

14.0 

0106+013 

01 06 04.3 

01 19 01 

2.63 

6.9 

MKN1152 

01 11 21.7 

-15 06 36 

1.67 

4.3 

0112-017 

01 12 43.4 

-01 42 59 

5.37 

15.3 

PG0117+213 

01 17 34.6 

21 18 05 

4.85 

18.6 

PG0119+229 

01 19 56.9 

22 54 36 

5.65 

20.1 

NGC 526 A 

01 21 37.2 

-34 19 33 

2.33 

6.3 

3CR48 

01 34 49.8 

32 54 20 

4.35 

9.5 

0205+024 

02 05 14.5 

02 28 41 

2.99 

8.3 

MKN590 

02 12 00.4 

-00 59 58 

3.07 

8.7 

NGC931 

02 25 14.4 

31 05 22 

7.07 

21.7 

0226-038 

02 26 22.2 

-03 50 55 

2.37 

6.5 

PHL 1377 

02 32 36.5 

-04 15 09 

2.29 

6.0 

0235+164 

02 35 53.0 

16 24 05 

7.60 

20.6 

PICS0237-233 

02 37 52.5 

-23 22 06 

2.23 

6.6 

0241^622 

02 41 01.2 

62 15 29 

73.30 

73.8 

H0323+022 

03 23 38.0 

02 14 47 

8.68 

24.8 

NRAO 140 

03 32 22.0 

32 08 3S 

14.22 

59.0 

H0414+009 

04 14 18.0 

00 58 03 

9.15 

26.9 

3C 111 

04 15 00.6 

37 54 19 

32.61 

63.0 

0420-388 

04 20 29.6 

— 3S 51 42 

1.91 

5.8 

0420-014 

04 20 43.3 

-01 27 24 

9.43 

39.2 

3C120 

04 30 31.6 

05 14 58 

12.32 

56.0 

3C 138 

05 18 16.9 

16 35 27 

23.28 

60.5 



PKS0521-365 

05 21 12.8 

-36 30 15 

3.37 

13.4 

0537-286 

05 37 56.7 

-28 41 26 

1.95 

3.6 

3C 147 

05 38 43.0 

49 49 45 

21.64 

45.1 

H0548-322 

05 48 48.9 

-32 16 60 

2.49 

3.9 

NGC2110 

05 49 46.3 

-07 28 01 

18.60 

53.4 

MCG8-11-11 

05 51 09.7 

46 25 51 

20.27 

39.7 

H0557-385(A) 

05 56 21.0 

-38 20 15 

3.35 

6.6 

0642+449 

06 42 52.9 

44 54 31 

11.61 

26.2 

3C 175 

07 10 15.9 

11 51 26 

11.51 

26.7 

MKN376 

07 10 35.7 

45 47 09 

9.32 

16.3 

PKS0735+17 

07 35 14.0 

17 49 09 

4.35 

9.2 

0736+017 

07 36 42.3 

01 44 00 

8.91 

26.5 

MI<N 79 

07 38 46.8 

49 55 46 

5.89 

22.3 

3C 186 

07 40 56.9 

38 00 32 

4.87 

9.7 

3C 191 

08 02 03.9 

10 23 56 

2.54 

4.7 

PG0804+76 

08 04 35.4 

76 11 31 

3.12 

12.2 

3C 196 

08 09 58.9 

48 22 09 

4.93 

14.3 

0830+112 

08 30 35.5 

11 15 30 

3.43 

7.7 

3C 204 

08 33 17.9 

65 24 04 

4.85 

18.3 

3C 205 

08 35 09.9 

58 04 53 

4.34 

20.8 

3C 206 

08 37 27.9 

-12 03 53 

5.85 

11.9 

3C 207 

08 38 01.9 

13 23 05 

5.40 

14.1 

PG0844+34 

08 44 33.9 

34 56 08 

3.39 

9.7 

3C 208 

08 50 23.0 

14 03 58 

3.60 

10.2 

OJ 287 

08 51 57.2 

20 17 58 

2.75 

5.4 

3CR215 

09 03 44.1 

16 5S 16 

3.75 

7.8 

0906+015 

09 06 34.9 

01 33 46 

3.38 

8.9 

MKN 704 

09 15 39.4 

16 30 59 

3.15 

6.3 

H091 7-074 

09 17 03.0 

-07 22 57 

3.30 

7.8 

0923+201 

09 23 05.6 

20 07 06 

4.16 

11.2 

PG0923+129 

09 23 20.1 

12 57 07 

4.03 

12.0 

4C39.25 

09 23 55.3 

39 15 23 

1.69 

1.8 

PG0934+013 

09 34 26.4 

01 19 12 

4.69 

12.4 

NGC 2992 

09 43 17.4 

-14 05 44 

5.56 

16.2 


PG1001+054 

10 01 43.3 

1011+25 

10 11 05.6 

1012+008 

10 12 20.7 

1020-103 

10 20 04.0 

B2 1028+313 

10 28 09.8 

3C245 

10 40 05.9 

3CR249.1 

11 00 27.3 

MKN 421 

11 01 40.5 

3C254 

11 11 53.3 

1115+080 

11 15 41.4 

PG1116+215 

11 16 30.0 

MKN734/PG 

11 19 10.9 

1121+422 

11 21 52.0 

MKN 180 

11 33 29.9 

3C263 

11 37 09.3 

PG 11 38+040 

11 38 42.4 

1146-037 

11 46 23.8 

NGC 4051 

12 00 36.3 

GQ COMAE 

12 02 08.9 

PG1211+143 

12 11 44.8 

ON235 

12 15 21.1 

1217+023 

12 17 38.3 

1219+305 

12 18 51.6 

MKN 205 

12 19 33.7 

1225+317 

12 25 55.8 

PG1229+204 

12 29 32.9 

NGC 4593 

12 37 04.6 

PG1241+176 

12 41 40.8 

PG 1244+026 

12 44 02.0 

3C277.1 

12 50 14.9 

1252+119 

12 52 07.5 

3C279 

12 53 35.7 

MKN 231 

12 54 04.9 

PG1307+085 

13 07 16.1 

NGC 5033 

13 11 09.1 


05 27 35 

1.88 

6.2 

25 04 09 

3.10 

16.3 

00 48 33 

3.22 

8.9 

-10 22 32 

4.89 

12.5 

31 18 20 

1.98 

2.7 

12 19 16 

2.70 

8.1 

77 15 08 

2.92 

6.9 

38 28 42 

1.45 

1.5 

40 53 41 

1.75 

1.7 

08 02 25 

3.61 

10.4 

21 35 42 

1.44 

2.2 

12 00 48 

2.74 

5.1 

42 16 52 

2.33 

5.2 

70 24 60 

1.27 

1.5 

66 04 26 

.82 

1.3 

04 03 39 

1.87 

3.5 

-03 47 29 

2.77 

4.9 

44 48 34 

1.31 

1.6 

28 10 54 

1.72 

2.6 

14 19 53 

2.83 

5.9 

30 23 41 

1.60 

2.3 

02 20 20 

1.97 

4.5 

30 27 15 

1.78 

2.9 

75 35 17 

2.74 

5.0 

31 4 r 12 

1.23 

1.1 

20 26 03 

2.58 

6.7 

-05 04 15 

1.97 

4.4 

17 37 26 

1.93 

7.5 

02 38 30 

1.93 

4.6 

56 50 38 

1.03 

1.1 

11 57 19 

2.63 

9.3 

-05 31 06 

2.22 

5.8 

57 08 36 

1.03 

1.3 

08 35 48 

2.20 

6.3 

36 51 29 

1.00 

1.9 



3C287 

13 

28 

15.9 

3C286 

13 

28 

49.9 

1331+170 

13 

31 

09.9 

MCG-6-30-15 

13 

33 

01.8 

PG1333+176 

13 

33 

36.5 

3C288.1 

13 

40 

29.8 

IC4329A 

13 

46 

27.8 

MKN 279 

13 

51 

51.8 

PG1352+183 

13 

52 

12.5 

PG1352+011 

13 

52 

25.7 

PG1402+266 

14 

02 

58.7 

PG1404+226 

14 

04 

02.7 

PG1407+265 

14 

07 

07.6 

NGC 5506 

14 

10 

39.1 

PG1416-129 

14 

16 

21.3 

3C298 

14 

16 

3S.9 

PG1425+267 

14 

25 

21.7 

PG1426+015 

14 

26 

33.7 

3C309.1 

14 

58 

58.0 

PG1501+106 

15 

01 

36.3 

PG1519+226 

15 

19 

02.0 

MKN 290 

15 

34 

45.3 

3CR323.1 

15 

45 

31.0 

1546+027 

15 

46 

5S.2 

PG1552+085 

15 

52 

19.0 

1611+343 

16 

11 

47.8 

TON256 

16 

12 

08.5 

MKN 876 

16 

13 

36.2 

MKN 877 

16 

17 

56.7 

3C334 

16 

18 

06.9 

3C336 

16 

22 

31.8 

PG1630+377 

16 

30 

15.1 

1634+706 

16 

34 

51.8 

1635+119 

16 

35 

26.0 

3C345 

16 

41 

17.6 


25 24 


30 

45 

58 

1.14 

2.2 

17 

04 

23 

1.76 

4.6 

-34 

02 

26 

4.06 

9.4 

17 

40 

31 

1.76 

4.7 

60 

36 

49 

2.09 

2.8 

-30 

03 

41 

4.55 

10.4 

69 

33 

13 

1.64 

1.4 

18 

20 

01 

1.84 

4.6 

01 

06 

51 

2.10 

6.1 

26 

10 

01 

1.42 

4.0 

22 

37 

59 

2.00 

6.4 

26 

32 

29 

1.38 

3.3 

-02 

58 

24 

4.22 

16.2 

-12 

56 

59 

7.20 

21.1 

06 

42 

21 

2.06 

7.3 

26 

45 

3S 

1.54 

5.1 

01 

30 

27 

2.64 

8.1 

71 

52 

13 

2.41 

1.8 

10 

37 

57 

2.23 

6.3 

22 

38 

22 

3.88 

19.1 

58 

04 

00 

2.32 

3.4 

21 

01 

2S 

4.04 

11.5 

02 

46 

07 

6.78 

28.9 

08 

31 

08 

3.47 

11.2 

34 

20 

21 

1.44 

3.4 

26 

11 

47 

3.77 

12.6 

65 

50 

37 

2.66 

3.3 

17 

31 

34 

4.35 

14.0 

17 

43 

31 

4.22 

13.4 

23 

52 

03 

4.52 

17.4 

37 

44 

11 

.90 

1.8 

70 

37 

37 

5.74 

11.7 

11 

55 

40 

4.29 

11.7 

39 

54 

10 

.74 

1.8 


MKN501 

3C351 

1721+343 

1725+044 

I ZW 186 

PKS1739+17 

1803+676 

3C380 

NGC6814 

PKS2121+053 

3C433 

2126-158 

2128-123 

IIZW136 

PKS2 134+004 

PHL1657 

NGC 7172 

BL LAC 

2201+315 

PG2209+184 

MKN304 

2216-038 

3C446 

CTA102 

NGC7314 

MR2251-179 

3C354.3 

PKS2251+113 

PKS2254+023 

NGC7469 

2304+042 

NGC7582 

PKS2344+092 

2345+184 

MIvN 541 


16 52 11.6 

17 04 03.4 
17 21 31.9 
17 25 56.0 
17 27 04.2 

17 39 26.7 

18 03 37.3 
18 28 12.8 

19 39 55.5 
21 21 14.7 
21 21 30.5 
21 26 26.7 
21 28 52.5 
21 30 01.1 
21 34 05.2 
21 35 01.1 

21 59 07.2 

22 00 39.2 
22 01 01.2 
22 09 30.1 
22 14 45.2 
22 16 15.9 
22 23 11.0 
22 30 07.7 
22 33 00.1 
22 51 25.7 
22 51 29.4 
22 51 40.5 

22 54 44.4 

23 00 44.2 
23 04 30.0 
23 15 3S.3 
23 44 03.4 
23 45 56.7 
23 53 28.0 


39 50 24 
60 48 32 
34 20 41 

04 29 29 
50 15 31 

17 21 59 
67 37 53 
48 42 40 

-10 26 33 

05 22 28 
24 51 18 

-15 51 48 
-12 20 20 
09 54 59 
00 28 26 
-14 46 26 
-32 06 36 
42 02 08 
31 31 06 

18 27 01 
13 59 26 

-03 50 33 
-05 12 16 
11 2S 24 
-26 18 31 
-17 50 54 
15 52 55 
11 20 40 
02 27 15 

08 36 15 
04 16 40 

-42 38 31 

09 l-_ 06 
18 27 29 
07 14 36 


1.73 

3.8 

2.26 

3.2 

3.06 

7.8 

7.03 

29.2 

2.58 

10.2 

5.44 

20.8 

5.00 

6.7 

6.60 

14.0 

9.80 

33.9 

6.27 

33.5 

8.75 

29.9 

4.85 

22.0 

4.83 

10.2 

4.20 

8.5 

5.12 

19.1 

4.45 

11.9 

1.65 

5.0 

20.15 

58.0 

9.91 

20.6 

4.82 

10.8 

5.23 

14.7 

6.18 

24.8 

5.26 

18.9 

5.05 

13.1 

1.45 

5.3 

2.84 

6.0 

7.13 

14.4 

5.53 

21.6 

5.88 

22.2 

4.82 

15.4 

5.52 

16.1 

1.48 

4.3 

5.05 

19.9 

4.25 

15.8 

5.02 

17.7 



a. For an assumed T, pin of 250 K (see text). _ 

1. Bright HI spectrum. Was not corrected for stray radiation (see text). 
Opacity correction dominates over stray radiation correction in most cases. 

2. Bell Labs survey was extrapolated to correct this spectrum. Error may 
be larger than normal. 


figure Is A low signal to noise 21cm spectrum from the NRAO 140ft. The integral cleaned 
N H of this spectrum corrected for stray radiation is 1.23 x 10 20 atoms cm -2 (1225+317). 


figure 2: A sample of the stray radiation spectrum (lower line) removed from the observed 
data (upper line). 


figure 3: Percentage error in N// due to opacity uncertainties (T, pin ) as a function of 
N w (140ft). 


figure 4: Uncertainty due to opacity uncertainties (N//(10000 K)- N//(125 K) ) as a 
function of peak brightness temperature (T peak ) in the 140ft spectrum. The horizontal 
line is the uncertainty due to other effects. Opacity uncertainties dominate above about 
4 x lO 20 atoms cm' 2 . 


figure 5: (a) Ratio, and (b) difference, between N w (140ft) and N w (Bell) as a function of 
Nj/ (140ft) in units of 10 20 atoms cm' 2 . 
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ABSTRACT 


The recent emphasis on ‘Big Bumps’ dominating the ultraviolet continuum of 
quasars has obscured the fact that bump properties vary widely, and that there 
are objects in which no such component is evident. As part of a survey of quasar 
continuum spectra, we have identified a class of quasars in which the optical- 
ultraviolet continuum ‘Big Bump’ feature appears to be weak or absent, relative 
to both infrared and x-ray. These ‘Weak Bump’ quasars are otherwise normal 
objects, and constitute a few percent of the quasar population. 

Subject headings: Quasars - spectrophometry - ultraviolet: spectra 



x. INTRODUCTION 


The ultraviolet ‘Big Bump ’ is the most striking feature in the continuum energy 
distributions of quasars (Malkan 1983). This bump is part .of* a larger *-*» < hatex> 
tends through the optical and, probably, soft x-ray bands (Amand et al. 1985, Czerny 
and Elvis 1987). The feature is thought to be a signature of thermal radiation rom 
the accretion flow (which is often modelled as a disk) at a few tens of Schwarzsch d 
radii from the central compact object (Malkan 1983, Bechto Id et at. 1987 and 
others). During a study of Big Bump properties in a sample of quasars with fu y 
observed energy distributions (100pm-4 keV), we found that in a minority of quasars 
(5 of the 31) the bump appears unusually weak or possibly absent. 


II. WEAK BUMP QUASARS 

The range of bump strengths in our sample is illustrated in Fig. 1 by three quasar 
energy distributions. The Einstein x-ray observations are plotted as i a bo - 
illustrating the best fit power law slope and 90% errors (see Wilkes an 
hereafter QED1, for details). Spectrophotometric and photmetnc observations are 
plotted as solid lines and individual points respectively. Optical s P ectr ° P ^° t °” e * 
is from Neugebauer et at. (1979,1987). IBAS AO photometry, near IR and °P t, “' 
photometry, and IUE spectrophotometry are from our QED Atlas, Elvis et at. ( 

preparation). 

The data have been corrected for Galactic reddening and converted to the 

. . r • it ^flkm s” 1 Mdc" 1 and Qo = 0.5; Emission lines have 

emitted frame assuming H 0 - 5Ukm s lvipc w . ’ , , r i ti 

been removed and a reddening correction has been made using a standard Gala 

extinction law (Savage and Mathis 1979) normalized using .small be ™^» 

obtained in the direction of the quasars (Elvis, Lockman and Wilkes 1989). In cas 

where substantial variability is evident in the data, an attempt has been made to 

keep only observations of similar dates. In most of our objects only mild variability 

is evident. 

3C 263 (Fig. la) has one of the most prominent bumps in our sample, 
while PG1116+215 (Fig. lb) is a more typical object; its ultraviolet bump is ess 
pronounced (It does have a large 3000A bump caused by Balmer continuum an 
unresolved Fe II line emission, Wills, Netzer and Wills 1985). In contract, Mkn 
876 (Fig lc) clearly has a much flatter spectrum than the other two objects; it 
fs a ‘weak bump’ quasar. Our sample contains four more ‘weak bump’ quasar^ 
shown in Fig. 2 that have no significant rise in their continuum flux toward the far 
ultraviolet other than the 3000Afeature. 


While the objects noted as lacking strong ultraviolet bumps by Edelson and 
Malkan (1986) and Ward et al.( 1987) were clearly reddened, we will argue that the 
five objects discussed here are intrinsically weak in the ultraviolet. Properties of the 
five weak bump quasars axe listed in Table 1(a). For comparison we list the other 
quasars from Fig. 1, and seven quasars whose energy distributions were published 
in Elvis et al. (1986) and Bechtold et at. (1986), in Table 1(b). The luminosities 
of the weak bump quasars are typical of the others in the sample, a few times 10 45 
erg s -1 in the near infrared, bright enough that contamination due to host galaxy 
starlight can be ignored. (Mkn 205 is less luminous, at a few times 10 44 .) Most of 
our sample objects, and our five weak bump objects in particular, vary by at most 
0.1-0. 2 dex (i.e. A log L v ~ 0.1 - 0.2) based on our observations and those from the 
literature (Table 1). Variability closer to 0.5 dex would be needed to simulate the 
weak bump energy distributions. 


To define ultraviolet bump strength objectively we use a far-ultraviolet to 
near-infrared color, 


Cuv/ir = log( 


L(0.1 — 0.2/im) 
L(1 — 2/xm) 


)• 


Fig. 3 shows the histogram of bump strengths for all quasars in our sample analysed 
so far. Despite the near-constant location of the 1/tm inflection (Elvis et al. 1986), 
there is a large range in the bump amplitudes. Most of our QED sample objects have 
strong ultraviolet bumps, with colors bluer than Cuv/ir = 0.4 (Table lb). However 
there is a tail in the distribution to redder Cuv/lRi ^ those objects picked out 
visually as having weak bumps are exactly those with C < 0.15 (Table la). Could 
it be that the infrared is enhanced rather than that the ultraviolet is weakened? 
It seems not. Three-way comparisons between the infrared (1-2/im), ultraviolet 
bump, and x-rays show that the weak bump quasars differ from other quasars by 
being weak in the ultraviolet rather than abnormally strong in the infrared. Also 
the (far-IR to far-UV) luminosities of the quasars range from 10 45 to 10 47 erg s" 1 
and there is no evidence within this sample for a correlation of bump strength with 


luminosity. 


III. LIMITS 01* INTERNAL REDDENING 

An obvious possibility is that the range in bump strengths is due to differing 
amounts of internal reddening in the sources. To test this, we construct color- 
color diagrams using ratios between the estimated luminosities in various frequency 
octaves. Fig. 4 is an (IR/visible, visible/UV) color-color plot. As noted earlier, 
the bump occurs at A < 1/im, so that the ratio of visual (4000 — 8000 A) and far 
ultraviolet (1000 - 2000A) luminosities largely measures its shape, while the near 



infrared (1 - 2,/m) to visual luminosity ratio measures their ‘power law continuum’ 
relative to the bump component. (We avoid the 2000 - 4000A region where the 
3000A bump is important, Wills, Netzer and Wills 1985). The reddening lines 
illustrate that the colors of the weak bump objects could be explained by reddening 
the ‘normal’ objects by a plausible E(B-V) of 0.1 to 0.2 magnitudes. 

The x-ray data make it unlikely, however, that such strong internal reddening 
is present. The IPC x-ray spectra for our objects are no different from those of strong 
bump quasars at the same level of radio loudness (QED1), and are inconsistent with 
the large hydrogen columns (> 10 21 cm- 2 ) that would be associated with the internal 
reddening (0.2 mag) discussed above. For all the weak bump objects no single 
power-law x-ray fit with this column is acceptable at the 90% level (see Fig. 2 in 
QED1). The presence of E(B-V) = 0.2 of reddening would imply that the intrinsic 
spectrum must be curved and extremely steep in the soft x-rays. For instance, for 
PG1613+658 the extra intrinsic column density would imply an energy index in 
the 0.1-4 keV range of a E > 2.5, and neither a single power law nor a two power 
law model gives a good fit to the data if both are absorbed by 10 21 cm 2 . The 
x-ray column density is unlikely to be smaller than the ultraviolet column since 
the x-rays are believed to come from a smaller region on energetic grounds. Indeed 
there is some evidence that the column to the x-ray emitting plasma may be greater 
than toward that emitting at other frequencies (Reichert et al 1985). We therefore 
conclude that reddening is not likely to be the cause of the bump weakness, unless 
an extra unabsorbed x-ray component is present (e.g. NGC 4151, Elvis, Bnel and 
Henry 1983, Pounds et al. 1986). 


IV. DISCUSSION 


If reddening is not the dominant effect, there are a number of alternate 
explanations for the diversity in bump strengths. There does not seem to be any 
correlation of bump strength with the other main continuum feature, the level of 
radio loudness, since two of the five weak bump quasars are radio-loud objects. It 
is possible that the bump may be highly variable in a given object; in this case all 
quasars may spend some time in the weak bump state, rather than there being a 
separate class of objects. This would provide strong constraints on parameters for 
accretion disk models of the bump . While some quasars do vary violently m the 
ultraviolet, (e.g. GQ Comae, Sitko 1986) in the five cases where we have multiple 
IUE observations, the variation is small compared to the range in bump strengths 
in the sample. Nevertheless this is a possibility and would be checked by repeated 
simultaneous optical and ultraviolet observation of the objects in Table 1. 


The r emainin g possibilities require intrinsic differences between the bumps in 
these objects and those in most quasars. For example, the bump may be present but 
peak at an abnormally high frequency, in which case its strength would anticorrelate 
with excesses in the soft x-rays. At present no quantitative measure of soft excess 
strength is available, and qualitative comparisons are inconclusive; several weak 
bump quasars seem to have no soft x-ray excess, but PHL909, the weakest of all, 
does show an excess (Masnou et al. 1989). Alternatively, the bump may peak at 
the same frequency but be intrinsically weaker relative to the power law continuum 
in these objects, compared to the objects in Table 1(b). This could be due either to 
a true luminosity change or to an inclination effect (c.f. Netzer 1985). The range of 
bump strengths relative to the total luminosity is a factor of 20-50, reasonable for 
an inclination effect. We note that the expected x-ray behaviour with inclination is 
uncertain, as the central region may become occulted at large inclinations. Finally, 
the slope of the bump component may be flatter in the weak bump objects. This 
could occur in accretion disk models if the disks in weak bump quasars were 
accreting using a different mechanism (with a different geometry?) from those 
with strong bumps and had a wider range of contributing temperatures. 

If the bump is really absent in some of our objects, what is causing the optical 
and ultraviolet flux that we do see (especially in the well observed 0.1 — 0.5/im 
range)? It has been argued (e.g. Sanders et a/.1989) that the infrared continuum 
is entirely due to dust, but dust could not contribute in the optical since the most 
refractory grains evaporate at around 2000K. A fairly weak bump and a thermal 
infrared continuum could conspire to give a continuous infrared to ultraviolet 
power law, but such a coincidence seems contrived. We suggest that the most 
probable explanation is still that the radiation is nonthermal in origin. The weak 
bump objects offer new possibilities for studying quasars and may allow us to 
examine the ‘bare’ non-thermal component in the optical and ultraviolet. 
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Table 1 

Object 

Coord name 

z 

Bump 

Strength 

Cuv/ir 

Observation date 
Near IR Optical IUE 

Optical 

Var. 

(dex) 

(a) Weak Bump Quasars 

PHL 909 

Q0054+144 

0.17 

-0.20 

1985 

1978 

1987 


3C 48 

Q0134+329 

0.37 

-0.01 

1975 

1978,1985 

1982 

0.1 

Mkn 205 

Q1219+755 

0.05 

0.15 

1986 

1972,1986 

1983 

0.1 

Mkn 876 

Q1613+658 

0.13 

0.14 

1985 

1980,1985 

1981/3 

0.1 

PHL 1657 

Q2135-147 

0.20 

0.09 

1985 

1969,1985 

1982 

0.2 


(b) Other Quasars 


PG 11 16+215 

Q1116+215 

0.18 

0.50 

1986 

1980,1986 

1982 

0.2 

3C 263 

Q1137+660 

0.65 

0.66 

1986 

1986 

1982 


PG1211 + 143 

Q1211 + 143 

0.09 

0.56 

1986 

1980,1985 

1982 

0.2 

3C 273 

Q1226+023 

0.16 

0.64 

Many 

1980, 1986 

1982 

0.3 

PG 1307+085 

Q1307+085 

0.16 

0.62 

1988 

1980,1986 

1980 

0.1 

PG 1416-129 

Q1416-129 

0.13 

0.59 

1988 

1980,1985 

1980 

0.2 

Mkn 1383 

Q1426+015 

0.09 

0.82 

1983 

1980,1986 

1983 

0.2 

Mkn 841 

Q1501+106 

0.04 

0.21 

1986 

1980,1986 

1983 

0.2 

3C 323.1 

Q1545+210 

0.27 

0.46 

1985 

1980,1986 

1983 

0.1 

References: Elvis et al (1986,1989), Bechtold et al (1987), Neugebauer et al (1979,1987) 



FIGURE CAPTIONS 


Figure 1: Sample rest frame energy distributions for three quasars, illustrating 
the range of ultraviolet bump strengths. The energy distributions are plotted 
as log {uL v ) against log„; (a) 3C 263 has one of the strongest bumps of any of 
our sample quasars, while (b) PG1116+215 is typical, (c) Markanan 205 is 
wpak bumo Quasar. 


FiRure 2: Rest frame energy distributions of weak bump quasars, (a) PHL 
909; (b) 3C 48; (c) Mkn 205; and (d) PHL 1657. Note the presence of the 
3000 A bump at \ogv = 14.5 - 15; it is especially prominent in (a) and (d). 


FiRure 3- Histogram of bump strength colors C IR/UV = logL(°.l 
0.2/xm)/L(l - 2 fxm) for a sample of 30 quasars, including those described in 

the text. Weak bump quasars are shaded. 


FiRure 4: Near infrared/optical/ultraviolet color-color diagram comparing the 
locus of sample objects with the direction of the reldemng line. The ' lines 
shows the effect of reddening PG1211+143 (center left) and 3C 273 (lower 
richtl by E(B-V)=0.1 and 0.2. Quasars are identified by their right ascensions, 
weak bump quasars are denoted by crosses, while quasars from Table 1(b) are 
represented by filled squares. 
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SUMMARY f 1 1 Rorri v rav selected active galaxies at 800 and 1 100 

throu^a spherical dust cloud, we fmd that to explain both the cont.nuum shape over 
S Zw 1 mm and the lack of silicate absotption, we requtre a region wrth r„. 10. 

models are the unexplained fine-tumng required, and a possible energy 
problem. 


1 INTRODUCTION 

Our picture of the spectral distribution of the continuous 
radiation from active galactic nuclei (AGN) has improved 
dramatically in the last few years, especially following the 
availability of far-infrared (far-lR) data from the IRAS satel- 
lite (Miley, Neugebauer & Soifer 1985; Neugebauer et al. 
1986 Slson & Malkan 1986; Ward « al. 1987). The 
authors of the present paper, along with other collaborators, 
have been involved in collecting multi-wavelength data for 
AGN seen in hard X-ray all-sky surveys, and ^particular for 
a complete flux limited sample taken from the HE AO- A. 
survey (Piccnotti e, al. 1982). Broad-band data from 00 
(IRAS) to 0.3 nm (l/-band) have been discussed by Ward 
et al. (1987) and Carleton et al. ( 1 987). VLA radio maps at 6 
and 20 cm have been discussed by Unger et al. ( 1 987 ). X-ray- 


spectra from 0.1 to 20 keV. measured with EXOSAT are 
discussed by Turner & Pbunds (1989). Ultraviolet (UV 
spectra, measured with IVE. will be discussed in Boisson et 
al. (in preparation). 

Comparison of far-lR and radio data shows that the milli- 
metre-wavelength region is of key importance, in the most 
general sense for delimiting the overall continuum shape^and 
specifically for critical tests of emission mechanisms. 1 his is 
especially true for radio-quiet AGN. where there is a drop of 
several decades in vF v between available far-lR and radio 
data (see the discussion by Lawrence 1987). Millimetre faci- 
lities have recently improved to the extent that observations 
of extragalactic radio-quiet objects are now feasible: at wave- 
lengths of 1 mm and shorter on the James c,erk ™ 

Telescope (e.g. this paper and Edelson et “ 

wavelengths of 1.3 mm and longer on the 1RAM facility 
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(Chini, Kreysa & Salter 1987; Chini et al 1988; Chini, 
Kreysa & Biermann 1989). In this paper we present 1.1- and 
0.8-mm measurements of three radio-loud and eight radio- 
quiet hard X-ray selected AGN (which together constitute 
roughly one-third of the complete Picinotti sample), and four 
other radio-loud AGN. 

We then proceed to discuss the overall radio to optical 
energy distributions, in the light of the new measurements. 
The most contentious issue remains whether the continua of 
radio-quiet objects are dominated by thermal or non-thermal 
emission mechanisms. It seems likely that varying combina- 
tions of both are required (Carleton et al 1987), but for the 
purposes of this paper we examine critically the conse- 
quences of extreme assumptions in either direction. 

It should be noted that throughout this paper, we present 
plots of log vF, versus log v. There are two reasons why vF v 
is preferred to F v for diagrammatic purposes: (i) vF v is energy' 
per unit log,, v, and so is easier for the eye to mentally inte- 
grate when plotted against log v; and (ii) over a large range in 
frequency. AGN energy distributions show vF,. constant to 
within an order of magnitude. Plotting vF x then enables one 
to reduce the vertical plot range, and so emphasize details in 
the shape. We shall, however, preserve the traditional nota- 
tion for spectral index in the positive form, i.e a is defined by 
F ( .= v a . Where applicable we assume a Hubble constant of 
/y„ = 50 km s ~ 1 Mpc “ 1 . 

2 OBSERVATIONS 

The observations were made during 1988 January 30 to 
February 1, on the UK-Netherlands-Canada submillimetre 
facility, the James Clerk Maxwell Telescope (JCMT), on the 
island of Hawaii. The JCMT is a 15-m segmented parabolic 
dish with chopping secondary. The secondary' mirror dis- 
tance is continuously adjusted as a function of telescope alti- 
tude and azimuth to keep focus during the homologous 
deformations of the dish. The photometer used was the 
common-user millimetre/submillimetre bolometer system, 
UKT14 (Duncan et al 1990). To reduce background noise 
and sky fluctuations, standard techniques of chopping and 
nodding were employed, with a beam separation of 60 arcsec 
in azimuth at a chopping frequency of about 7 Hz. 

The UKT14 system is equipped with filters that are rela- 
tively narrow and well matched to the atmospheric bands. 
This both reduces sky-noise, and minimizes colour effects in 
calibration. The effective frequencies of the observations are 
given by the expression 

f* vF(v) D{ v)dv 
v = , 

J i. F(v)D\v) dv 

where F[v) is the source spectrum, and D(v) is the responsi- 
vity curve, including telescope, detector, filters and atmos- 
phere. For reasonably good conditions, with Precipitable 
Water Vapour (PWV ) = 1 mm, and assuming a source spec- 
trum F,. = constant x v a with a - + 2, the effective frequen- 
cies (wavelengths) of the T 100’, ‘800’ and ‘450* broad-band 
filters are 277 GHz (1082 ^m), 385 GHz (779 ^m) and 672 
GHz (446 ^m), respectively. Changing the assumed source 
spectral index from a=+2toa=-l changes these values 
by less than 10 per cent in all cases. Most of our calibrators 
and targets have spectra rising steeply towards short wave- 


lengths, so we choose the effective frequencies appropriate 
to PWV= 1 and a= +2 to be the frequencies at which we 
estimate a monochromatic calibration. The uncertainty 
caused by colour effects will almost certainly be smaller than 
our main uncertainties due to (i) determination of 
atmospheric optical depth, and (ii) pointing errors. 

The primary' calibrators used were the planets Jupiter and 
Mars. The quasar 3C 273, bright enough on the JCMT to be 
visible in real time on the strip chart, was also used to check 
the estimation of optical depth, on the assumption that it did 
not vary' significantly in the course of a night. Planet fluxes 
and angular sizes for the dates in question were obtained 
from the JCMT fluxes program, as discussed by Duncan et 
al (1990). To perform the calibration at our chosen frequen- 
cies, we assume that both planets have spectral index 
a - + 2. We also assumed that all of our targets were point 
sources. (This may not be quite true, as for some of the radio- 
quiet AGN, dust emission from the parent galaxies may be 
comparable with the nuclear emission). The fluxes from 
extended planetary’ sources therefore require beam correc- 
tion. To perform this correction, we assumed that the planets 
were discs of uniform surface brightness, and that the tele- 
scope beam was of Gaussian shape, with Half Power Beam 
Width (HPBW)= 19 arcsec at both 1100 and 800 pm. As the 
observed ratio of signals from Mars and Jupiter agreed w ith 
prediction within a few' per cent, we have reasonable confi- 
dence in both the beam correction and the overall calibra- 
tion. 

At the time of our observations, the system zero-point was 
not very well known, so we determined both zero-point and 
atmospheric optical depth from observations of Jupiter and 
Mars at a range of airmasses. We assumed a plane parallel 
atmosphere, and so determined calibrations between instru- 
mental units and F v of the form log^i Jy/m V ) — A + B x 
airmass. The optical depths B on the two good nights of the 
run were found to be ~ 0.05 and 0.15 for 1100 jum and 
- 0.3 and 0.65 for 800 /<m. From the scatter in these calibra- 
tions, we believe that changes in optical depth may lead to 
errors in the derived fluxes of typically 10 per cent up to 20 
per cent. 

The other major uncertainty arises from pointing errors. 
The brighter objects allowed ‘peaking-up’ before a photo- 
metric measurement. (In subsequent data reduction we dis- 
carded a few objects measured in this way. where we felt that 
the peaking-up measurements were not of high enough sig- 
nal-to-noise ratio.) The procedure for blank fields was to 
peak-up on the nearest mm-bright object and then to offset to 
the target. Planets were peaked-up by hand. Other bright 
objects were measured at five positions - nominal pointing 
position, and north, south, east and west bv half a beam- 
width. Local pointing offsets were then de~:*>ed from the 
fitted profile. To ensure greater accuracy, this procedure was 
repeated until a good centred fit was obtained. However, 
there remain systematic errors in the telescope pointing 
it lei. and secular drifts over periods of time. Since our 
o rving run in 1 988 January/February, both kinds of prob- 
lem have been substantially reduced. Given the beamsize of 
~ 19 arcsec HPBW, it is necessary to maintain pointing to 
within - 3-4 arcsec to keep flux calibration errors within 10 
per cent. Experimentation on the nights in question showed 
that this required re-pointing roughly every fifteen minutes, 
which was done. Even so, we cannot, of course, completely 



dismiss occasional larger errors. The upper limits quoted in 
this paper are derived from many observations over 
extended periods at a variety of telescope altitudes and 
azimuths, so that completely missing a detectable source is 
very unlikely. Nonetheless, occasional pointing problems for 
blank-field sources may have led to substantial net errors. 

The primary targets were AGN selected by hard X-ray 
emission. Eleven of these were measured. We also examined 
four other radio-loud AGN. The results are summarized in 
Table 1(a) and (b). The signal-to-noise and formal error are 
based only on the statistical noise, estimated from the scatter 
within individual integrations; it should be remembered that 
the real uncertainties are of the order ~ 30 per cent for 
bright objects, and possibly substantially larger for blank 
fields. Where the formal significance of the measurement is 
below' 5 a, we have also quoted an upper limit. Because of 
the possible pointing problems, we quote a conserv ative 3- a 
limit, rather than the more usual 1-or 2 - o limit. 


2.1 Comparison with other measurements 

Observations of the radio-quiet members of our X-ray 
selected sample represent the first ~ 1 -mm wavelength 
measurements of these objects deep enough to allow detailed 
discussions of their emission mechanisms. [Note, however, 
that Edelson et al. (1988) have measured NGC 4151 at 438 
fi m using the same system.] Three radio-loud objects in the 
sample have been measured previously by other workers. 
The quasar 3C 273 has been measured many times and is 
clearly variable at mm wavelengths (Elias et al. 1978. Ennis. 
Neugebauer & Werner 1982; Clegg eial. 1983. Robson etal. 
1983; Landau ei al. 1986; Robson ei al. 1986; Courvoisier el 
al. 1988; Robson et al.. in preparation). We found it at a 
fairly typical flux level. It seems likely that the 1-mm 
measurement by Landau et al. (1986), at 69.9. Jy, is in error, 
especially as their 1-mm measurement is higher than both 
their 3-mm and their 350-//m measurements. Our 
measurement of 11IZW2 is lower than, but statistically con- 
sistent with, the 4-o measurement of Ennis et al. (1982). Our 
measurement of MKN 501 is also lower than that of the 
latter authors, but in good agreement with Joyce & Simon 
(1976). assuming a spectral index of -0, as is typical for 
blazars. 

For our other targets, PKS1413 and OJ287 are in good 
agreement. Our measurement of 3C 446 is somew hat lower 
than the UKIRT measurement made by Gear et al. (1985) 
but again the difference is statistically marginal, as the older 
measurement is of low significance. This object is known to 
be variable, as shown by the results of the monitoring cam- 
paign by the Preston group (Gear et al. 1985; Brown et al. 
1986; Brown et al. 1 989a.b ). 

3 RADIO TO OPTICAL CONTINL A 

We now combine our new millimetre wavelength 
measurements with the radio core data (from Unger et al. 
1987, and references therein), with IRAS data, with Kuiper 
Airborne Observatory (KAO) data, and with ground-based 
near-mid-lR and optical broad-band data (all from Ward et 
al. 1987, and references therein). We also add some data 
from pointed IRAS observations (Elvis et al.. in preparation s. 
For NGC 4151 we also include far-IR and submillimetre 
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Table 1 . (a) Hard X-ray selected AGN. (b) Other AGN. 


(a) 


RESULTS FROM THIS PAPER 


Previous measurements 


0) 

(2) 

(3) 

(4) 

(5) 

NAME 


date 

F*| m Jy ) 

3<? limit 

N2992 

779 

30/1/88 

80±24 

<152 


1082 

30/1/88 

47±34 

<149 

N3227 

1082 

31/1/88 

22±42 

<148 

N4151 

779 

30/1/88 

5±17 

<55 

N 4-593 

1082 

31/1/88 

31 ±22 

<97 

N5506 

1082 

31/1/88 

11±1? 

<62 

N7314 

1082 

31/1/88 

-20±29 

<87 

IC4329A 

1082 

31/1/88 

3±14 

<45 

3C 273 

1082 

30/1/68 

15150 



1082 

31/1/88 

15996 



779 

30/1/88 

11666 


MKN 50 1 

1082 

30/1/88 

411 ±34 



779 

30/1/86 

375±21 


MKN59Q 

1082 

31/1/88 

-34 ±21 

<63 

11IZW2 

1US2 

31/1/68 

207 ±22 


lb) 

RESULTS FROM THIS PAPER 


(6) 

(7) 

(8) 

(9) 

OtH-mn, 

+2.63 

A(^m) 

F„(mJy) 

ref 

> +2.03 

> +2.47 

438 

<200 

16) 

> +1.76 

> +2.15 

> +2.13 

> +1.56 

3000 

<270 

(7) 

-0.6 

1000 

see text 


-0.6 

1000 

800 ±200 

m 

> +134 

3200 

370±60 

fb.i 

-0.16 

1000 

800 ±200 

(U 


Previous measurements 


(1) 

(2) 

NAME 


3C446 

1082 


P KS 1413 

1082 

OJ2S7 

1062 


(3) (4) 

dale F„(mJy) 

31/1/68 27S3±6S 

31/1/88 1734 ±14 

31/1/88 3331 ±63 


(5) lb) 
3 <t limit 

800 

1100 

2000 

1100 

1100 

1000 

1000 

1000 

3300 


17/ (8) 

F t (mJy| ref. 

3900± 1 000 (4) 

6$00±1700 (4) 

5800 ±4 00 (4) 

1700±4UU 14) 

4000± ) 000 (4) 

4900 ±500 ( 1 ; 

8090 (5) 

3500 - 6000 |2 1 

26403:290 (3 1 


NR. A 0530 1082 30/1/88 1563 ±41 

30/1/88 1251 ±31 


Notes to columns: 1 2 j effective wavelength as explained in text: i3; 
local calendar date at start of night: 4 error quoted is statistical 
error only. Systematic uncertainty in fluxes is of the order 30 per 
cent: (5) limn is formal estimate or 0 i whichever is greater + act. i 6 
spectral index between 100 um and 779 or 10S2 /im as appropri- 
ate. For MKN 501. have used 60-,um flux, as 100-^m measurement 
is an upper limit. ( 8 1 References: 1 1 ! Ennis et al. 1. 1 982 ): ( 2 ^Roellig et 
ai ( 1 986); (3) O’Dell et al. , 1 978 : (4 Gear et al. 1 1 985 Landau 
et al. (1986); (6: Edelson et ai <1988;: i'7i Blitz. Mathieu & Ball) 
( 1986 j; (8) Joyce & Simon 1976/. 


data from Engargiola et al. (1989), and Edelson et al. (1988 ). 
Except where noted in the text. optical-near-lR data are 
nuclear point-source measurements, having been corrected 
for extended parent galaxv starlight by the methods 
described in Ward et al (1987). Whether this correction 
actually removes all the starlight will be discussed in Section 
4.3. The broad-band measurements will be contaminated to 
some extent by emission lines, but as we are here considering 
the gross shape of the continuum over many decades in 
frequency, this problem should not be too important. It 
should also be remembered that the measurements at various 
wavelengths are not in general simultaneous, so there may be 
some distortion of the continuum shape. This is much more 
likely to be a problem for the radio-loud objects in the 
sample, which all have ‘blazarMike cores. 

The continuum energy distributions are shown in Fig. 
Ka)-(c) corresponding to groups classified as explained in 
the subsections below. In this paper we are concentrating on 
the radio to optical continua but. for reference, we also indi- 
cate the vF r level at 6 keV. 
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Figure 1. Ratlin through optical enerpv distributions for the I I target galaxies. Data points at l .l and 0.8 mm are from this paper, other data as discussed in text. The vertical axis has been 
arbitrarily scaled in each case to accommodate several objects on one graph I riangles indicate upper limits I he dotted lines are present only to connect points referring to the same object, and 
should not be taken as an indication of the intervening continuum. Horizontal dotted lines indicate the vi\ level of 6-keV X-rays. The targets are divided into three groups as discussed in the 
text. 



3.1 Radio-loud objects 

Three of the 11 objects show energy distributions which are 
continuous throughout the optical-radio region (see Fig. la). 
These, of course, are all objects which have traditionally 
been classified as ‘radio-loud’. The overall shape could be 
described either by a continuously curving spectrum (Lan- 
dau et al. 1986), or by a series of power laws with self- 
absorbed synchrotron components in the mm to cm regime 
(Brown et al. 1989a). To second order, there are clear com- 
plications to the simple picture. The quasars 3C 273 and 
11IZW2, show UV turn-ups, as is typical for quasars in 
general. The BL Lac object MKN501 shows a peak in the 
near-IR, presumably due to residual starlight - note that in 
this case the near-IR data are from uncorrected 5-arcsec 
aperture measurements. 3C 273 shows another peak at 4 
which could well be due to hot dust in or near the broad 
emission-line region (Courvoisier etal. 1988). 

So we see that, although the overall smoothness of the 
energy distribution over many decades of frequency is con- 
sistent with the idea that the continua of these objects are 
dominated by non-thermal radiation of some kind, in detail 
the explanation of these energy' distributions may not be alto- 
gether simple. However, for these objects there is strong 
independent evidence for relativistic effects and generally 
non-thermal nature: superluminal motions, high polarization 
and large-scale variability. What we will stress below is that 
such supporting evidence is noticeably absent in radio-quiet 
objects. 

3,2 Radio-quiet objects 

The energy distributions of the remaining eight objects in the 
sample all show a dramatic fall from 100 ^m to 1 mm (see 
Fig. lb and c; Table la). Where the spectral index is greater 
than 2.5, this does not strictly rule out synchrotron emission 
as is commonly assumed, as steeper slopes are possible (de 
Kool. Begelman & Sikora 1989; de Kool & Begelman 1989). 
A break of some kind between the far-IR and radio data was, 
of course, previously required, but our new data refine the 
position of this break by two decades in frequency. Regard- 
less of the nature of the 1R continuum, these objects may all 
be called radio-quiet in the objective sense that the vF v level 
at 6 cm is typically 2-3 decades below that at 6 keV, com- 
pared to the objects in Fig. 1(a), where the difference is typi- 
cally 1 decade. We caution, however, that radio-loudness 
may be a continuous variable, rather than a bi-modal quan- 
tity (Unger et al. 1987; Lawrence 1987). These objects, in 
contrast to the radio-loud objects, are not known to show high 
polarization, dramatic variability, or superluminal motions, 
so any postulated non-thermal nature of the optical-IR con- 
tinuum is on rather weaker ground. In particular, of course, 
given the steep fall from 100 gm to 1 mm. it is likely (though 
not certain) that the far-IR emission is dominated by thermal 
emission from dust. Similar conclusions have been arrived at 
recently for optically selected radio-quiet quasars, following 
measurements at 1.3 mm (Chini et al. 1989). 

For the objects in Fig. 1(b), the case for thermal emission 
in the far-IR is fairlv strong. In all these cases, vF v is rising to 
long wavelengths through the far-IR, with a slope typical of 
normal spiral galaxy disc ‘cirrus’ emission (Rowan-Robinson 
& Crawford 1989). These are all relatively nearby, low- 
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luminosity Seyfert galaxies, so the large IRAS apertures are 
dominated by normal galaxy emission rather than the rela- 
tively weak nuclear emission, as discussed by Ward et al. 
(1987). The level of any underlying non-thermal continuum 
emission may be indicated either by the X-ray vF,. level 
(assuming the spectral index to be roughly - 1 ), or by the low 
points in the 1-100-^m continuum, as argued by Carleton et 
al. (1987). Such a power-law component would still need to 
become self-absorbed by - 1-200 ^m, if it is not to exceed 
our millimetre upper limits. We also note that there is a 
tendency for these objects to show peaks around 1-2 ^<m. 
despite the starlight subtraction’ performed by Ward el al. 

(1987). . 

The objects in Fig. 1(c), unlike those in Fig. 1(b). are not 
obviously dominated by parent-galaxy interstellar cirrus 
emission. The steep fall from the far-IR to our mm upper 
limits requires that a power-law component dominating in 
this region must certainly cut-off by 1 00 /<m. (More stringent 
limits, of a few tens of ftm, will be discussed in Section 4.1 ) 
However, it is not at all clear that such a power-law 
component is dominant. Both N4151 and IC4329A show 
clear peaks around 25 ,um, and curve slow ly down thereafter, 
strongly suggestive of a thermal component. Both N5506 
and IC4329A show downturns in the optical. This is almost 
certainly due to reddening of a few magnitudes, consistent 
with the observation of large Ha/H/J ratios, and small Ha/ 
X-rav ratios (Ward et al. 1987). 

For N4151, N5506, and 1C4329A. the mid-lR vF, level is 
a factor of several above the 6-keV level, whereas for radio- 
loud objects it is level or below, apparently indicating an 1R 
excess' in radio-quiet objects. Alternatively, this could be 
interpreted as an ‘X-ray excess' in radio-loud objects. Indeed, 
it seems that radio-loud quasars are more X-ray luminous for 
a given optical luminosity, than are radio quiet quasars 
(Zamorani et al. 1981; Kembhavi. Feigelson & Singh 1986). 
possibly due to an extra flat X-ray component in radio-loud 
objects (Wilkes & Elvis 1987). At first glance. MKN 590 
seems to be a flat power law, but in detail it has a turn-up 
from 60 to 100 ft m (suggesting a galaxy disc component) and 
also shows a weak maximum at 25 ^m. 


4 DISCUSSION 

The importance of dust is clear in objects with reddened 
optical continua or steep IR continua. However, even in 
radio-quiet objects with flat energy distributions (MKN 590. 
NGC 4151), the presence of local peaks and curvature rules 
out pure power-law models. Contributions due to dust 
heated by the UV excess source, and by interstellar cirrus 
and star-formation regions in the parent galaxy, must at least 
be taken into account. However, non-thermal emission could 
still be the dominating component in many wavelength 
regions. In Section 4.1 we fit maximal power-law models to 
MKN 590 and NGC 4151, and discuss briefly whether the 
physical characteristics of such synchrotron sources are 
reasonable. Given the lack of direct evidence (polarization, 
variability) for non-thermal emission, we might conserva- 
tively assume that all emission is thermal. The immediate 
problem for such a position is that around 1 nm neither dust 
emission nor the UV excess will contribute significantly. In 
Section 4.2 we discuss the possibility that the required extra 
component is starlight from a nuclear star cluster. In Section 
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4.3 we then proceed to compare the energy distributions of 
NGC 4151 and MKN 590 with multi-component thermal 
models. A reasonable match requires rather fine tuning of 
model parameters. 

4.1 Maximal power-law model fits 

Examples of maximal power-law fits are shown in Fig. 2(a) 
and (b). There is, of course, no unique solution unless models 
for the obvious excess components are also included, but the 
range of models shown gives us an indication of constraints 
on the self-absorption frequency, which in either case must 
be A a <75 /4m, and probably much less if we wish to mini- 


mize the mid-IR excesses. We can then derive physical 
characteristics of the synchrotron source; diameter d , bright- 
ness temperature 7~ h , equipartition magnetic field B. and 
ratio R of photon to magnetic field energy densities. In the 
discussion below' we assume that the synchrotron source has 
a filling factor of unity, and that the ratio of proton to 
electron energy content is 100. 

For NGC 4151, / A = 7 5 /4 m implies d=l 0 light-hours, 
i? = 1 48 G, T h = 4.8 x 1 0 10 and R = 1 0. If we do not assume 
equipartition but instead that 7 h = 10 i: , as usually seems to 
be the case for compact radio sources (Kellerman & Pauliny- 
Toth 1981 ), then d= 2 light-hours, B= 543 G and R = 10\ 
In either case, most energy should emerge via Compton scat- 


(a) mm IRAS LKHJ opt-UV 
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Figure 2. Self-absorbed homogeneous synchrotron models compared to the observed energy distributions. The turnover region is modelled as 
F„ * constant ( v/r' A ) 5 -(l - e~ '*), where the optical depth to self-absorption follows r v * ( v/v A ) a ' 5 2 . (a) NGC 4151. Models shown have 
spectral index and self-absorption wavelengths of a = - 1.2, A a =75 /4m and -1.4, A a =54 //m. (b) MKN 590. Models shown have 
a = -1.2, A a = 24 um and a~ -1.0. A a = 75 /4m. 
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tering rather than synchrotron radiation. The observed 
X-ray luminosity of NGC 4151 is similar to the observed 
synchrotron luminosity, but as the spectral index of the X-ray 
emission is ~ 0.55 (Yaqoob, Warwick & Pounds 1989), it is 
patently not the corresponding Inverse Compton compon- 
ent. For MKN 590, A a =75 ptm implies d = 2 light-hours, 
B= 1.2 x 10 3 G, T„=6.2x 10" and R= 100. 

With Hq = 100, T b would be a factor of 4 lower, and d a 
factor of 2 lower. The large uncertainties in proton/electron 
ratio and filling factor result in uncertainties in d and B of 
factors of several, but not orders of magnitude. The very 
small sizes derived here therefore constitute a problem for 
synchrotron models, as an extensive study shows no vari- 
ability in the IRAS bands for radio-quiet AGN (including 
NGC 4151 ). to a level of 15 per cent over one year (Edelson 
& Malkan 1987). In contrast, the same study demonstrates 
strong variability in several radio-loud objects. 


represent a - 50 per cent dilution of young stars, i.e. red 

supergiants (Terlevich el al. 1990). 

The continuum at - 1 ftm cannot be entirely due to stars, 
as there is some variability, but it is consistently smaller than 
the blue-light variability, even for quasars (Cutri et al ■ - i 

For example, in the case of NGC 4151, Lebofsky & Rieke 
(1980) describe observations in an 8.6-arcsec aperture where 
the C-band varied by a factor of 2.5. but the R-. /-. 7-bands 
bv only 25 per cent. (In this example, variability increased 
again at longer wavelengths.) If the wavelength dependence 
of variability amplitude seen in NGC 4 1 5 1 is due to dilution, 
we can conclude that the variable component contributes 
only 14 per cent of the light in the low-state, or 35 per cent in 
the high-state. 


4.2 A nuclear starlight component? 

The data for NGC 4151 and MKN 590 that we use have in 
principle been corrected for parent galaxy starlight (Ward l et 
al. 1987). This involved modelling CCD photometry with a 
point source of arbitrary colours, and an exponential disc 
and r 1 ' 4 bulge with fixed colours representative of a late-type 
stellar population. However, if a nuclear star cluster is 
present, this will have been included in the point-source 
component. Even in quasars, an inflection point in the con- 
tinuum consistently occurs around 1 and noticeably, the 
minimum is deepest for the most luminous quasars (Knss 
1988) It seems possible therefore that a point-source ste ar 
component is common in AGN. Note that even at the rela- 
tively small redshift of NGC 4151 (z- 0.0033), requiring 
that such a starlight component be unresolved ( ~ 1 arcsec) 
onlv requires it to be - 100 pc across. [Ebstein, Carleton & 
Papaliolis (1989) find from speckle observations that NGL 
4151 has an unresolved component of size 0.3 arcsec, but 
this is at 5000 A.] The luminosity in such a stellar compon- 
ent for NGC 4151 is roughly equivalent to 10 K0 giants. 
Such a nuclear star cluster may provide the material for the 
broad-line clouds, and the fuel for accretion (Hills 1975, 

Penston 1988). . . . 

The obvious stellar feature to look for is the Ca n triplet at 
- 8542 A. Unfortunately Ca n often appears in emission in 
broad-line objects (Persson 1988). It certainly appears in 
absorption in NGC 4 1 5 1 . but there is some controversy over 
the strength. Malkan & Fillipenko (1983) give a value of 
1.1 A for the equivalent width of one component of the trip- 
let, Ca n 8542 A. in the unresolved nuclear seeing disc, com- 
pared to an average value for spiral galaxies of 3.2 A. Recent 
higher resolution work by Terlevich. Diaz & Terlevich 
(1990), however, indicates an equivalent width oi 
These authors quote an equivalent width for the whole triplet 
of 6 A - not very' different from the average value they find 
in normal early-type spirals (7.7 A). Similar results are found 
by Terlevich ei al. for two other Type 1 Seyferi galaxies, 
NGC 3227 and 3516. This suggests a more or less undiluted 
normal old stellar population of red giants, in contrast to the 
substantial dilutions of the G-band and Mgb features seen at 
shorter wavelengths (Malkan & Fillipenko 1983; Terlevich el 

al 1990) Alternatively the large Ca u equivalent width could 


4.3 Purely thermal model fits 

Here we consider all components that may be present to 
some degree - the UV excess, interstellar cirrus, star-forma- 
tion regions, starlight, and dust heated by the UV excess. T e 
latter is of course, the most controversial and interesting 
part, so we develop a physical model for this region, whereas 
for the other components we use templates of fixed shape. U) 
The UV excess is modelled as a single blackbody at 
30000 K. (ii) Following Rowan-Robinson & Crawford 
(1989) we model the interstellar cirrus as the sum of two 
grevbodies. at 2 10 and 30 K. which provides a good fit to the 
spectrum of the interstellar medium in our own Galaxy 
(Boulanger, Baud & van Albada 1985; Hauser et al. 1984). 
We should note, however, that the 12-,um excess may be due 
to verv small grains, so lhat in a galaxy with a powerful 
Sevfen nucleus, the hotter component may be absent. (m» 
For the star-formation component, we use the model derive 
by Rowan-Robinson & Crawford (1989), which fits the IR- 
mm spectra of both hot-centred clouds in our own Galaxy, 
and the starburst ring in NGC 1068. (iv) The nuclear starlet 
component is modelled by the sum of blackbod.es at 3000 
and 5000 K, which we find adequately fits the bulge starhgh 

colours quoted by Ward etal.{ 1987). 

Finally we consider a more detailed model of the dust 
intimately connected with the AGN itself. We assume that a 
dust shell surrounds and is heated by the UV excess source, 
with a density that falls off with radius as a power law. We 
then numerically solve the radiative transfer through this 
region (using the same software as Rowan-Robinson & 
Crawford 1989) in order lo predict the detailed emergent 
IR-mm spectrum. The grain properties assumed are as in 
Rowan-Robinson (1982). The variable parameters are then 
the grain condensation temperature, the density-law slope, 
the UV optical depth, and the ratio of inner to outer radii. To 
first approximation, the overall spectral slope is fixed by the 
density law and optical depth. The melting temperature fixes 
the high-frequency cut-off, the outer radius fixes the low- 
frequency cut-off, and the optical depth also fixes the silicate 
feature strength. As well as matching the overall energy dis- 
tribution. we aim to minimize the strength of the 9. n 
silicate feature, as this is normally weak m Seyfert galaxies 
(Aitken & Roche 1985). 

A fit to NGC 4151 is shown in Fig. 3(a). We assume tha 
the far-lR is dominated by a strong starburst. and that most 
of the light in the 3-band is due to starlight. What remains is 
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Figure 3. Thermal models compared to the observed energy distributions, (b) NGC 4151. The components used (discussed in more detail in 
the text) are the 30000-K blackbody, starlight, a star-formation region, and a dust shell with density law following r ~ 1 and UV optical depth 
~ 12. (b) MKN 590. The components used are the 30000-K blackbody, starlight, interstellar cirrus, and a dust shell with r~ 1 density law and 
UV optical depth ~ 6. 


reasonably well-fitted by a dust-shell model with melting 
temperature 940 K, density following r” \ Ri/R 2 * 0.001, and 
r uv * 12. The optical depth needed to explain the continuum 
shape is close to that needed to minimize the size of the 9.7- 
jum silicate feature, but it is difficult to get these both perfect 
simultaneously. A fit to MKN 590 is shown in Fig. 3(b). Here 
we assume that the far-lR is dominated by cirrus, and that 
starlight dominates at J. What remains is well fitted by a dust- 
shell model that in some resj ;ts is very similar to that 
required for NGC 4151. Once again the density must follow 
r~\ and the optical dep u is of a similar order (6). However, 
the melting temperature . equired is slightly higher, 1 380 K, 
and RJR 2 ^ 0.01 differs from NGC 4151 by a factor of 
10 . 


With present broad-band data quality then, we can fit 
purely thermal models, but the physical properties of the 
dust shell need fine tuning. With future higher resolution 
data, the confrontation between model and data should 
become crir il, as is obvious from the small-scale curvature 
of the mode, spectra. In particular, it is impossible to get the 
silicate features in both the 10- and 20-^m windows simul- 
taneously minimized. From our model fits we expect to see 
the 20-/im feature weakly in emission. We therefore strongly 
encourage 20-//m spectroscopy of Seyfert galaxies. The most 
obvious contradiction in the modelling is that a spherical 
dust shell of such high optical depth would completely 
extinguish the UV source. One solution would be a patchy 
dust shell, and another would be a thick disc of dust, as 


postulated to exist in the Seyfert 2 galaxy NGC 1068 to 
explain the optical spectropolarimetry (Antonucci & Miller 
1985). Two of us (AE and MRR) have recently developed a 
two-dimensional numerical code for solving the problem of 
radiative transfer through dust in the axisymmetric situation 
(Efstathiou & Rowan-Robinson 1990). This will be applied 
to the Seyfert galaxy problem in future work. 

The deduced luminosities in the heated dust component, 
and the observed UV component (up to 1500 A) in NGC 
4151 are -5x10” and 1 x 10 4 ’ erg s' 1 , respectively. The 
XUV continuation of the UV component is extremely un- 
certain, and could result in several times more energy being 
available for heating, but a large covering factor would still 
be required to achieve consistency. However, if an optically 
thick disc model, such as envisaged for NGC 1068 turns out 
to be feasible, the energy available for heating could be very 
much larger, as hard X-rays will also be absorbed (via 
Compton scattering). The (2-10 keV) luminosity of NGC 
4151 is variable, but of the order -5x10 - erg s 
However, the y-ray continuation suggests a total luminosity 
possibly as large as - 10 4 - (Lawrence 1980). 
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ie respectively by requiring a particular optical depth, and 
by requiring a fairly constant covering factor. In considering 
thermal versus non-thermal phenomenological models, we 
are then left then with the choice of (i) arbitrarily introducing 
an extra component not suggested by independent evidence, 
i e the synchrotron component, or (ii) keeping only to com- 
ponents definitely present, but having to fine-tune he 
properties of those components, with no good physical 
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5 CONCLUSIONS 

New millimetre measurements have enabled us to examine 
the shape of the continuum from optical to radio wave- 
lengths for hard X-ray selected AGN. Objects with strong 
radio cores are almost certainly dominated by non-thermal 
emission throughout the IR. mm and radio regions, but also 
show evidence for several secondary- components and effects 
- the UV excess, reddening, residual starlight and some ther- 
mal dust emission. Objects with weak radio cores show a 
steep fall from the far-lR to mm wavelengths, as has also 
been recently found for optically selected quasars (Chini et 
at. 1989). Part of the emission from these objects is very 
likely to be thermal, but it is empirically possible to model 
the continuum shapes with and without a synchrotron 
power-law component. Such a synchrotron power-law com- 
ponent must, however, become self-absorbed by a few tens 
of urn This implies a size of a few light-hours tor the syn- 
chrotron source, in conflict with the lack of far-IR variability. 
Furthermore, such a synchrotron source should be heavily 
Compton-dominated. Confronted with these problems, and 
the fact that there is no supporting evidence for non-thermal 
phenomena, one might take the conservative stance that for 
such radio-quiet objects all the emission is thermal. The con- 
tinuum shape longward of 2 Jim can be duplicated by a com- 
bination of interstellar cirrus, star-formation regions and a 
dust region heated by the observed UV source, if the dust 
shell has a density law following r'.An extra component is 
then needed, peaking around 1 /zm. We suggest the existence 
of a nuclear star cluster, and argue that this idea is supported 
by the presence of Ca n absorption at 8542 A. and by the 
fact that variability is a minimum at 1 ^m. 

The strongest arguments in favour of a power-law com- 
ponent underlying the continua of radio quiet objects are (i) 
the smoothness of the 8-13-/zm spectrum, and (n) the rela- 
tive constancy of the overall continuum shape from object to 
object (reflected in the good correlation between near-lR 
and X-ravs. and in the consistency of spectral indices defined 
by ‘baseline’ techniques). Both of these facts can be explained 
in purely thermal empirical models, but only by fine tuning , 
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Abstract 


Two Einstein and three EXOSAT observations of PG1211 + 143 over a six 
year baseline show that strong steep-spectrum low energy x-ray emission is a 
persistent feature of this quasar. EXOSAT observations of PG1211 + 143 detected 
a factor 2.3 increase in its steep soft x-ray (C-band) flux in 18 days. The bulk of 
the soft x-ray emission of PG 121 1 + 143 therefore comes from a region < 5 x 10 16 cm 
across. In another (quasar rest frame) interval of 193 days the soft x-rays decreased 
by a factor 3.7. Either an intrinsic emission variation or an increase of N// from 3.4 
to / ,6x lO 20 atoms cm -2 could produce the same change. In the same time interval 
the hard (2 - 10 keV) x-rays decreased by a factor 1.6+0.05. This suggests a 
connection between the two energy regimes, and argues against variable absorption 
causing the soft X-ray variations. In order to decrease in luminosity so rapidly a 
thermal source in PG1211 + 143 would have to be optically thick to both electron 
scattering and free-bound absorption. The same argument applies to 6 of the 7 other 
rapidly variable soft excesses reported in the literature, suggesting that optically 
thin thermal models are ruled out in general. 



1. Introduction 

‘Soft X-rav excesses’ are a new component of the quasar continuum. At 
energies below the Carbon edge at 0.28 keV (>44A, the ‘C-band’) half of all studied 
quasars and AGN are dominated by large excess flux over an extrapolation of eir 
higher energy power-laws (Arnaud et al. 1985, Wilkes and 

Pounds 1989, Turner et al. 1990, Masnou et al. 1991, Comastn et al. 1 ) • 

has been suggested that these soft x-ray/extreme-ultraviolet excesses are emitted 
from the inner edges of accretion disks (Arnaud el al. 1985, Bechtold et al 1 . 

Czerny and Elvis 1987). In this case their radiation comes to us from the innermost 
identifiable region of a quasar. The spectral shapes of the soft X-ray excesses might 
then give us information on the peak temperatures in disks and hence determine 
the allowed values of the central masses and accretion rates. 

Whatever their origin this new component deserves careful investigation. 
Existing constraints on their spectra are poor because the X-ray measurements 
cover a small energy range (-0.1 - 0.3 keV) and have almost no energy resolution 
Also very few objects were observed more than once so that variability informati 
scarce. Clearly more studies of these ultra-soft excesses are needed Variability data 
in particular will be most important to distinguish between extended and compact 
emission regions, and to limit the properties of accretion disk models (S.emiginowska 

and Czerny 19S9). 

PG1 9 11 + 143 is an extreme case of a quasar with an soft x-ray e-vcess. It 
has an excess extending to higher energies than in other cases (-2 keV Bechtold 
et al 19S7), and is the only one of the 33 quasars with soft x-ray spectra in h 
Einstein IPC survey to be dominated by an extremely steep, soft x-ray spectrum 
(Wilkes and Elvis 1987). We therefore undertook a program of x-ray o serva ions 
of PG 121 1 + 143 with the EXOSAT Observatory. This paper reports this series ol 

observations. 

The value of the absorbing Galactic column density is crucial to studies of 
soft X-ray excesses (see . Elvis, Wilkes mrd McDowell 1990). We use i the accurate 
Galactic column density of 2.8 ±0.1 xl0*> atoms cm" from the ' study by El , 
I orltman and Wilkes (1989). Values of the Hubbl. constant H 0 =50 km s Mpc 
altZ = 0 are used throughout. The redshift ot PG1211+143 is 0.085 (Schmidt 

and Green 19S3). 


precise ‘Soft X-ray excess’ which could apply to any energy from ~0.2 to 


2. Observations 

PG1211 + 143 was observed on three separate occasions by EXOSAT with 
intervals of 209 and 19 days between the observations. (For details of the EXOSAT 
instrumentation see de Korte et al. 1981 and Turner, Smith and Zimmerman 1981). 
Table 1 gives details of the observations, including the derived count rates. 

The quasar was seen in all of the three energy bands defined by the ‘thin 
lexan’ (3LX), ‘aluminum/parylene’ (ALP) and ‘Boron’ (BOR) filters, which each 
cover the ~0.1 to 2 lceV band with somewhat differently weighted responses, in front 
of the LE (Low Energy Channel Multiplier Array) during the first observation. In 
the later observations the source was fainter so the filters giving smaller effective 
areas (BOR, ALP) could not always be used. There were no serendipitous EXOSAT 
observations of PG1211 + 143. 

The position of the source in all the LE detections axe consistent with one 
another and with the optical (Schmidt and Green 1983) position of PG1211+143 
to within the 6 arcsec uncertainties of the EXOSAT aspect (Osborne and Angelini 
1986). The LE source appears point-like in ail the exposures. The mean LE position 
for the source is: a = I2 h ll m 44.9 s ±6”; 6 = 14° 19' 53"±0. s 4 (1950.0). 

The LE count rates (Table 1) were estimated using a square box centered 
on the source centroid and of a size which maximizes the signal-to-noise. The 
background was measured from a nearby source-free area of the image and rescaled 
to the source position to take into account spatial non-uniformities in the LE 
(Giommi 1985). The count rates were corrected for the fraction of counts outside 
the box (Davelaar and Giommi 1985), for vignetting effects, and for telemetry and 
instrumental dead times (Osborne 19S5). 

On the first two occasions the source was clearly detected in the ME (Medium 
Energy) Argon detectors. Poor background subtraction in the third observation 
makes detection by the ME uncertain. The data from the ‘inner’ detectors of the 
ME were much noisier than the ‘outer’ detectors since they axe less well shielded 
against paxticle background (Arnaud et al. 1985). We therefore used only the 
‘outer’ detectors in our analysis. Background subtraction for the ME was performed 
according to the method of Smith (1984). The background in each case was taken 
from the same detectors during times when they were pointed away from the source 
during an ‘array swap’. Periods of high count rate due to solar flares were excluded 
by visually inspecting the light curves of each detector ‘half’. The source is detected 
over the PHA range 7-40 which corresponds to 2 - 10 keV. 

A serendipitous observation of PG1211 + 143 (seq no. 6982) was made by 
the Einstein IPC (Imaging Proportional Counter, ~0.2-3.5 keV, Giacconi et al. 
1979, Gorenstein, Harnden and Fabricant 1981). This IPC observation was made 



in 1979, a year earlier than the targeted observation. This significantly extends our 
temporal baseline. In this observation the quasar was 21.8 arcminutes off-axis, 
such positions the gain (‘BAL\ which measures the conversion from recorded pulse 
height to incoming photon energy) must be estimated using the spatial gam map 
(‘DGNI\ Harnden et al. 1984). The ‘BAL’ was 13.6 for this position (on a scale from 
_ ^19). The IPC count rate has been corrected for vignetting and the loss o 
counts due to the degraded off-axis point response function. The targeted Einstein 
IPC observation of PG1211 + 143, reported by Elvis, Wilkes and Tananbaum (1 ), 

Bechtold et al. (1987) and Wilkes and Elvis (1987), had a 20% higher count rate. 

3. Persistence: A 6-year soft excess 

The unusual feature of PG1211+143 is its steep soft X-ray spectrum It 
is important to understanding this emission to know whether it is persistent or 
transient. The five EXOSAT LE and Einstein observations allow us to test this. 

We have fitted a simple power-law plus absorption model 2 to each of the 
data sets. Table 2 gives the results of these fits. (The spectrum from the first 
observation is reported by Comastri et al. 1991). The errors on each fit are quite 
large but the error regions overlap in both spectral parameters. It is notable that the 
well-constrained fits (i.e. excluding the 79/339 and 86/008 data) give ste ^P e ^^ s 
for the power law index than are typical of either low energy quasar spectra (W likes 
and Elvis 1987), or of higher energy AGN spectra (Mushotzky , urner an 
Pounds 1989). When the Galactic column density is imposed as a mminuim v ue 
then slopes of around a E ~2.0 are required, even for the 79/339 (off-axis IPC) a a^ 
However the x 2 values even for the best fit values (Table 2) are not good suggesting 
that the spectrum has curvature, steepening toward lower energies. 

The LE data for 1985/164, having three data points, allow a fit to be made^ 
This fit suggests an even steeper slope, ~3.0. A black body fit to the same 1985/164 
LE data gives a best fit temperature of 75±5 eV . 

The unusually steep spectrum of PG1211+143 is thus not a transitory feature 
but is seen in observations spanning six years. 


2 a v~ aE e~ SH ° , cross-sections (c), were 


taken from Morrison and McCammon (1983). 


4. Change: Optically thick thermal emission 

On timescales of days to months significant variability is seen between each 
of the EXOSAT observations. In particular, in the quasar’s frame, the LE 3LX data 
show a factor 3.7 ±0.3 decrease in 193 days and a factor 2.3 ±0.1 increase in 18 
days (table 1). In the 193 day interval the hard (2 — 10 keV) x-rays also decreased, 
though only by a factor 1.6±0.05. This may indicate some physical link between 
the emission in the two energy regimes. 

The variations may be due either to intrinsic emission variability or to 
changing intervening absorption. At these low energies quite small column densities 
can produce such changes. To produce the decrease by a factor 3.7 in the 3LX filter 
seen between the first (Jun 19S5) and second (Jan 1986) EXOSAT observations 
needs an increase in N h from 3.4 x 10 20 atoms cm -2 (the best fit value for the first 
observation) to 7.6 x 10 20 atoms cm -2 . This is a minimum value. Material with 
a larger column density could obscure up to half the source to produce the same 
effect. 


The contemporaneous decrease in the higher energy ME count rate by a 
factor 1.6 argues against absorption as the only variable since 2-10 keV X-rays are 
negligibly affected by column densities less than 10 21 atoms cm -2 . 

If the XUY excess were due to an optically thin thermal plasma then its 
spectrum would be rich in diagnostic emission lines for the inner regions of quasars 
(Marshall 1990). However the variability information rules out an optically thin 
plasma quite strongly, assuming an intrinsic change in the emission, since it requires 
large optical depths to both bound-free and electron scattering in the emitting 
plasma, as shown below. 

If we take the decrease in 3LX count rate between the first and second 
EXOSAT observations to define a minimum cooling time, we can derive an optical 
depth for the emission region, taking into account electron scattering within the 
source, which smears out variations and increases the observed timescales by a 
factor (l + r es ). Of course more rapid variations may have taken place between 
our observations, but the shorter implied cooling times only increase the required 
optical depth derived below. For simplicity we consider a uniform density sphere of 
constant temperature plasma. 

A spherical optically thin source of radius r with electron density n = T es /roT 
and temperature T has luminosity in the 0. 1-0.3 keV range of 

I = -xr 3 F 0 fn 2 T' /2 B(T) 

u 



r 


= ^r^Fo/T^exp 


O.lkeV 

kT 


where 


the bandpass factor 


F 0 = 2.4 x 10 -27 erg cm 3 K 1/2 s \ 


B(T) = exp 


O.lkeV 


kT 


which is sensitive to the uncertain effective lower frequency hunt of the soft excess 
observation, and f is a factor giving the additional ^ntr.but.on_of hne «™ss.on to 
bremsstrahlung cooling, 1 < / < 10 (Raymond and Smr.h 19. ., B*ymond 1990). 
The radius in this equation can be related to the minimum cooling time bj th 

variability, t var by 

Cl V ar 

r£ (TT77 

(Optically thick sources vary slowly for their size, so the inferred size i mailer for 
high depth) and so, combining the above equations 



1 + T„ 


> 9 


where, 


9 = 


3 Lot 


4tz ct vaT FofT l / 2 B{T) 


r \-i/2 


exp 


- 1 - 2 (nc) 


-n 


= 7.1 x 10 2 


10 44 erg s 


-) r ' fe)’‘ (*ot) 1/1 exp \~ 0M (ftwk 


If o < 1 it is the square of the electron scattering depth. If 9 1’ 

assumption of optical thinness leads to a contradiction and we may deduce that the 
“urce must be optically thick; , is then equal to the electron scattering depth as 

long as the formuia for L remains valid. 

Even if the source is optically thick to electron scattering, it 

interactions are dominant in modifying the spectrum toward a black body The 


bound-free optical depth r b j is just r e (where the a are the corresponding cross- 
sections). The bound-free cross-section at these sub-kilovolt energies depends only 
weakly on temperature in the range above 1 million degrees; cr b j/a e3 at 0. 1-0.2 keV 
is in the range 1-5 with T between 1 and 10 million degrees (Krolik and Kallman 
1984, Kallman private communication) and so the spectrum will indeed be modified 
towards thermal form if the source is optically thick. The upper limit to the optical 
depth obtainable from comparing the blackbody radius to the variability timescale 
is not interesting. 

The decrease in 3LX count rate between the first and second EXOSAT 
observations of PG1211-I-143 gives a value g of ~14, implying that it is optically 
thick. If instead the flux decrease is due to changing absorption in front of the source 
then electron scattering delays do not apply. The variability gives the physical size 
directly if we assume an absorber local to the quasar which moves to occult at least 
the fraction of the source whose luminosity disappears. For a transverse velocity of 
the absorber. r, the optical depth to electron scattering is then ~ 180(-)“i, which 
is considerable even for a relativistically moving intervening absorber. 

Variability of soft x-ray excesses is now known to be common in quasars 
and AGN (see summary in Elvis, Wilkes and McDowell 1990). In some cases the 
variations are far more rapid than in PG1211-H43 ( e.g . Mkn 335, Turner and 
Pounds 198S, Lee and Balick, 19S8). Table 3 lists the seven cases in the literature 
for which luminosities and timescales for luminosity decreases are given. (Only 
decreases in luminosity were used since only these define cooling times.) These 
variations imply the values of g listed in Table 3 and are shown in Figure 1 against 
lines showing the y=0.5 (r=l) and g = 10 loci for electron scattering. 

PG1211 + 143 is seen from Figure 1 (and Table 3) to be one of the least 
extreme of these soft X-ray excesses. Most are factors of 10 to 100 times more 
optically thick, to the extent that large values of g imply a large optical depth. 
(E1615+061 is the only exception but here the excess may not be real, Buckley et al. 
19S7). We conclude that optically thin emission rarely dominates the production 
of quasar soft X-rav excesses, so that strong line emission from this component is 
unlikely. Optically thick emission would tend to support an origin in the central 
regions of an accretion disk, although we note that the physics of the inner edge 
are not well understood (Shakura and Sunyaev 1973), and that other processes 
{e.g. Comptonization, Czerny and Elvis 1987) could well modify the emitted disk 
spectrum. 

The changes in the LE/3LX flux of PG1211+143 between the EXOSAT 
observations trivially rule out a source for the soft excess x-rays in extended regions 
on the scale of the optical/ultraviolet narrow emission line cloud zones (~100 pc- 
1 kpc), such as suggested by Halpern and Filippenko (1984). They would also 



have ruled out emission from the broad emission line regions (Krolik, McKee and 
Tarter 1981 ) under earlier assumptions about their size, however recent observations 
(Peterson ct ai 1990) suggest that sizes of 5-30 light days axe qmte p aus.ble for 
these regions, comparable to the light travel time size limits derived for the soft 

excess emission regions. 


5. Conclusions 


Over five observations covering six years the steep low energy x-ray spectrum 
of PG1211 + 143 was always present. In this quasar the XUV excess is a persis en 

feature. 

EXOSAT observations of the quasar PG1211+143 have detected large am- 
plitude (factor >2) short timescale changes in its low energy X-ray emission, an 
increase in IS days and a decrease in 193 days. The small size imp ie or 
bulk of the soft x-ray emission from PG1211 + 143 implies large optical depths m 
the emitting plasma and so makes the presence of strong emission lines unli 'e y. 
This conclusion is true of most known variable XUV excesses m AGN and quasars. 
Optically thick emission is instead likely from the inner regions of an accretion dis 
although modifications (due to e.y. Comptonization) can occur, and the detailed 

phvsics is not well understood. 

Observations of variability and spectra of the XUV excess component of 
AGN continuum emission are proving valuable. ROSAT observations are likely to 
extend this work substantial!} . 
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figure 1: Luminosity in the ‘soft X-ray excess’ component us. observed timescale for a 

decrease in the soft X-ray component for the AGN from table 3. The lines show values of g 
= 1 (marginally optically thick) and g — 10 for thermal emission with electron scattering. 
Larger optical depths lie to the bottom right, implying that most soft excesses are optically 
thick. 



Table 1: Einstein and EX OS AT Observations of 

PG1211 + 143 


year/day instrument/filter exposure (s) counts 


1979/339 

IPC 

1980/347 

IPC 

1985/164 

ME a 

LE/3LX 

LE/ALP 

LE/BOR 

1986/008 

ME 0 

LE/3LX 

1986/028 

LE/3LX 

LE/ALP 


b. PHA 7-40, 2-10 keV 


6306 0.89± 0.02 

1795 1.28± 0.03 

22370 0.84±.02 

1257 0.206±.012 

3539 0.0867±.006 

16560 0.0087±.0012 

17030 0.52±.02 

15753 0.055±.0025 

10065 0.13±.004 

15753 0.05±.009 


Table 2: Power-law fits to X-ray data for PG1211+143 



1979/339 IPC'-® 1*318:1 

1980/346 IPO 7 2.0lo6 

1985/164 ME h 

1985/164 ME h +LE 

1985/164 LE 2.9to.9 

1986/008 UE h 0.9±H 

1986/ 008 ME^ + LE 0.9to:4 


0.9^‘j 

4 If 

12, 5 

0.04 

2.3 toi 

5.5+°;* 

12, 5 

0.04 

0 lo 5 


43, 31 

0.08 

9 0 + 1-5 

4-8tU 

56, 34 

0.01 

3.4 Ifl 

2.4tll 

0, 0 

— 

o.o tli 

-i ■» +5.5 
1-0.6 

47, 31 

0.04 

0.3lg1 

1 1+ 1 ' 3 
A ' A -0.4 

47, 32 

0.04 


a. 90% confidence ranges for 3 interesting parameters (xLn + 6 ‘ 25 ’ Avni 1976 )’ 

b. 10 20 atom cm -2 . 

c. /xJy at 1 keV. 

d. degrees of freedom. 

e. probability of obtaining the observed x • 

f. PHA 2-9. 

g. errors include IPC gain uncertainties at this off-axis position. 

h. PHA 7-40 


Table 3: Optical Depths for Variable Soft Excesses in AGN 


Name log L° s log At variability g h reference 

{erg s~ l ) s factor 


Mkn 335 

43.6 

4.3 

2.0 

FairaU 9 

44.9 

5.7 

1.9 

MSI 

40.5 

2.7 

2.0 

3C273 

45.2 

7.2 

1.7 

NGC 4051 rf 

41.5 

3.0 

2.0 

PG 1211 + 143 

44.7 

7.2 

3.7 

NGC 5548 

44.6 

4.4 

1.3 

Mkn 841 

44.4 

4.6 

1.5 

El615+061 e 

44.0 

8.0 

1.5 


1000 Turner, Pounds 1988, Lee et ai 1988 
450 Morini et ai 1986 
30 Barr et ai 1985 c 

65 Turner et ai 1990, Masnou et ai 1991 
300 Lawrence et ai 1985 
14 this paper 

8100 Kaastra and Barr 1989 
3200 Arnaud et ai 1985 
0.5 Pravdo et ai 1981 


а. luminosity in the soft excess, ^0.1 - 0.3 keV. 

б. for T = 5 x 10 6 K. Values are ~5 times smaller for T = 1 x 10 6 K, while r^/ is ~85 times larger. 

c. quoted in Barr and Mushotzky 1986. 

d. note the warnings of Lawrence et al. (1987) about picking timescales from 1/f variations. 

e. but see Buckley et al. (1987) who note a nearby x-ray bright RSCVn system. This may be the 
cause of the apparent soft excess in this AGN. 




log t (s) 



log L (erg s _1 ) 
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Abstract. 


We study the low energy excess above a high energy power-law in the X-ray spectra 
for a signal-to-noise limited, sub-sample of 14 quasars from the Wulkes &: Elvis 
(1987) sample. Detailed analysis of the Einstein Imaging Proportional Counter 
(IPC) data, combined with Monitor Proportional Counter (MPC) data where 
possible, and the use of accurate Galactic N H values allow us to estimate the 
flux of any detected excess and improve constraints on the spectral parameters 
at higher X-ray (> 0.6 keV) energies. We find a significant excess in 8 of the 
14 objects, confined in all but one case to below 0.6 keV. The excess is typically 
strong, dominating the <0.3 keV X-ray flux and being 1—6 times as strong as the 
high energy component at 0.2 keV. In 3C273, multiple observations show that the 
excess is variable. The lack of spectral information, uncertainty in the intrinsic 
line-of-sight absorption and the possibility of variation prevents the detection of 
any relation with the optical/ultra-violet blue bump. Comparison with an existing 
accretion disk model for the blue bump suggests that, for this model to explain 
both the blue bump and the soft X-ray excess, a range of disk parameters axe 
required to explain the whole sample. 


Key words : Quasars — X-Rays:spectroscopy — X-Rays:general 
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1 Introduction. 


Recent soft X-ray observations have shown that Active Galactic Nuclei (AGN) 
and quasars commonly show strong excess emission above the extrapolation of 
their higher energy slopes (Singh, Garmire & Nousek 1985, Arnaud et al. 1985, 
Branduardi-Raymont et al. 1985, Wilkes & Elvis 1987, Turner & Pounds 1989, 
Elvis, Wilkes &: McDowell 1989). These excesses occur primarily in the C— band, 
i.t. below the 0.28 keV Carbon edge, and rise steeply to lower energies. Informa- 
tion is generally limited to a single discrepant point at low energies when compared 
with the slope at higher energies (1 — 10 keV) ( t.g . Turner & Pounds 1989). In this 
paper we use the highest signal-to-noise Imaging Proportional Counter (IPC) ob- 
servations from the Wilkes &z Elvis (1987) quasar sample, combined with higher 
energy (2—10 keV) Monitor Proportional Counter (MPC) data where it is use- 
ful, to give data extending from 0.1 — 10 keV. These data combined with accurate 
Galactic Nh values (Elvis, Lockman Wilkes 1989), allow two-component spectral 
fits and thus a well-constrained estimate of the observed flux in the soft excess. 

Although there are viable alternatives, such as free-free emission (Barvainis 
1990, Ferland, Korista and Peterson 1990), the optical/ultra- violet (UV) blue 
bump is generally modeled as thermal emission from an accretion disk possessing 
a range of temperatures (Shields 1978, Malkan & Sargent 1982), an interpretation 
which is enticing both observationally and theoretically. An obvious identification 
of the soft X-ray excess is then the high-energy tail of the thermal emission (Arnaud 
et al. 1985, Bechtold et al. 1986, Pounds et al. 1986, Czerny & Elvis 1987). In 
this case it would originate in the hottest material, at the innermost edge of the 
accretion disk at ~ 10 Schwarzschild radii. Should this be confirmed, the excess 
will provide far stronger constraints on the physical conditions of this part of the 
disk them have been possible so far. Quantitative measurement of the soft excess 
is thus of considerable importance. 


2 Analysis. 

2.1 The Data. 

W T e have selected the highest signal-to-noise observations in the Wilkes Elvis 
(1987) sample of quasars observed with the Einstein IPC, those with >1200 net 
counts. This yielded a subset of 14 quasars, 5 of which are radio-loud (radio- 
loudness, Rl > 1.0, Wilkes Elvis 1987). Details of the X-ray observations along 
with redshift and Galactic column density of neutral hydrogen are given in Table 1. 
PG 1211 + 143, which is a member of the sample, is included for completeness in 
the discussion although it is not re-analyzed here. We include two additional 



observations of 3C273 not reported by Wilkes ft Elvis (1987). For completeness 
and comparison with the earlier results, the results of smgle power-law fits wrth 
free Nh to the IPC data alone for these two observations, using Brown & Gou 
1970) absorption coefficients, are given in Table 2a and Figure 1. Throughout 
the remainder of our analysis, absorption coefficients are taken from Morrison k 

The IPC X-ray data were analyzed following the standard procedure de- 
scribed in Elvis el al. (1986) and Wilkes & Elvis (1987) Briefly counts ,n the 
energy range -0.15-3.5 keV were obtained from a 3' arde centered on the posi- 
tion of the source. A correction was applied for counts falling outside this arc e 
and background counts were estimated from a 5'-6' annulus centered upon the 
source position. The first bin of the pulse height distribution was excluded since 

the errors are not well determined. _ _ t , . , e nn 

In this new analysis we have also used data from the MPC for which > 500 

net counts were detected. Apart from early studies by Halpern (1981) and Elvis 
e( al. (1986), the MPC has been little used for quasars until now due to systemati 
uncertainties in the background subtraction process for these faint sources^ Tese 
uncertainties have now been reduced and quantified for all observations prior to 
Oct 1980 (Arnaud 1991 and see below) allowing simultaneous fitting of IPC ana 

MPC data. „ • , * 

The MPC data have the advantage of covering a higher energy range wi 

the counts di : buted into 8 pulse height channels extending from 1.2 to 20 keV. 
Channels 7 ana 8 (> 10 keV) were not used in this analysis as they provide no 
additional constraint on the spectral fits. The energy of our 

ally covers the two decades 0.1-10 keV. In most cases the MPC and IPC d 
were obtained simultaneously. The exceptions axe PG 0026+129, II Zw 136 and 

PKS 0637— 752 where no significant variation was detected between t e wo a es 

and non- simultaneous data were used to maximize the signaJ-to-noise ratio. 

2.2 The Monitor Proportional Counter (MPC). 

The MPC was a collimated proportional counter filled with Argon and sealed with 
a 38 nm window of beryllium foil. Its effective area was 667 cm with a triangular 
field of view having full width at half maximum of 43 and a flattened i top > ~ 
wide Counts were collected in 8 pulse height (P :hannels from 1.2-0 keV 
with an energy resolution of ~ 20 7c at 6 keV. The instrument was mounted on the 
outside of the Einstein Observatory telescope and aligned within 1 of the op ica 
axis. The MPC is described in detail by Gaillardetz et al (1978) and Gnndlav 

ti al (^°J sed the latesti Revls time-ordered MPC processing with data extracted 
following the method of Arnaud (1991). Briefly, all data taken during points in the 
orbit with Mcllwain l> 1.2, B< 0.25, which includes South Atlantic Anomaly pas- 
sage were discarded. The background consists of two components. A backgroun 



model for the cosmic background was constructed using source-free observations. 
The background is estimated using correlations with the pulse risetime discrimi- 
nator rates and Mcllwain l,B, parameters modified by long-term time variability 
which includes the decay of the in-flight calibration source (Cd 109 ), modeled as an 
exponential with a half-life of 450.4 days). In addition the response matrix was 
modified following the discovery that the earlier, RevO, processing system assumed 
an energy resolution ~ 10 times better than reality. The resulting effective area 
curve is shown by Vrtilek et al. (1990) and Arnaud (1991). 

2.3 IPC/MPC Relative Normalization. 

The IPC and MPC instruments were calibrated independently and the estimated 
accuracy of these calibrations are about 10% and 20% respectively. Before combin- 
ing data from the IPC and MPC, the relative accuracy of their calibrations must 
be determined. Since the instruments have very different sensitivities they did not 
observe the same calibration sources. Instead we made a comparison utilizing the 
quasar data under the assumption (consistent with Ginga results, Pounds 1990) 
that quasar spectra do not change drastically (A a > ±0.5) in the 0.6-10 keV 
energy range. 

A single power law model, with the line-of-sight absorbing column density 
fixed at the Galactic value, was fitted to the high-energy band of the IPC (0.6-3. 5 
keV) to minimize the effects of any soft excess. A power law model with the 
IPC best fit slope and Galactic N H was then fitted to the MPC data but with 
the normalization (flux density at 1 keV) free. These fits should give comparable 
normalization in both instruments. In all but the highest S/N observations, those 
of 3C273, the ratio of MPC to IPC normalization (Figure 2) is consistent with 
unity. However a weighted mean gives a value of IPC/MPC = 1.14±0.04, similar 
to those for 3C273 and implying that the MPC calibration is too low relative to the 
IPC. This conclusion is reinforced by a comparison of the best fit power law slopes 
determined by fitting the IPC and MPC data individually and in combination 
for the two 3C273 observations where this is possible (Table 2b and Wilkes &: 
Elvis 1987). The IPC/MPC combination gives a steeper slope (a£~ 0.6) than 
either instrument individually, which each gave o ;£~ 0.4. The latter is in better 
agreement with measurements of 3C273 by other experiments (Worrall et al. 1979, 
Turner et al. 1990a, b). Since the ratio of IPC to MPC normalization is not sensitive 
to the spectral slope assumed, our result does not critically depend on the accuracy 
of the spectral fit. 

We conclude that the relative normalization of the MPC is systematically 
low relative to the IPC by a factor of 1.14. This is within the absolute accuracy 
of the calibration of each instrument. The MPC data were corrected accordingly 
in our subsequent analysis and all numbers presented here include this correction. 
We note that the decision of which instrument to correct is arbitrary. 



2.4 Spectral Analysis. 

The results of conventional, single power law models, with Nh ted Galactic 

value, are presented in Table 3. They axe poor (P(x>> <10%) for 9 ob^rvabons 
of 7 of the 13 quasars analyzed here. In 11 observations the residuals clearly d. 
the lowest 1 or 2 PHA (pulse height analysis) channels to be sigmfic y 
the fit, suggesting the presence of a soft excess. 

Conventional two-component fits give poorly constrarned parameters (slop 

and normalization at 1.0 keV, f(1.0 keV)) for the low energy 
el al. 1985, Pounds e( ai. 1986). However, the excess flux can be much better de 
termined than its spectrum since a small change in flux at 0 2 keV results in a large 
change in slope and normalization (f(1.0 keV)). We have therefore applied broken 
power-law fits to the data, minimizing with respect to the 0.2 ke ux in eac 
component and to the 1.0 keV flux of the high energy component ( f*s(0.- . k ), 
f„ E (0.2 keV), f(1.0 keV) respectively ) instead of the more conventional spectrd 
slope. The results are given in Figure 3 and Table 3, with errors quoted at la 
for one interesting parameter, A* J = 1.0 (Avni 1976). N„ was aginn ted atjhe 
Galactic value. The energy resolution and S/N of the data are ins 
mnt he JddLn of a fourth free parameter and thus it is not possible to determine 
"ak energy between the two power laws. The IPC is best suited to detecting 
soft excesses in its lowest independent energy band i.e. with break energy in 
range 0 4-0.6 keV. The break energy was thus fixed at 0.6 keV, yielding a conser- 
vative estimate of the observed flux in any detected excess. Decreasing this v Jue 
to 0.4 typically results in a small but insignificant « 2a) increase m the deduced 

flux of ^ s “ f ‘ h e ““® r Uw plus black . bo dy components were also applied leaving 

free the 0.2 keV flux in each component and the normalization of the high energy 
component The energy resolution is not sufficient to fit the black-body temp 
ature which was therefore fixed at 0.05 keV. Both the high energy ““P°« nt 
andthe soft excess flux axe consistent with corresponding result from the broken 
power law fits as shown in Figure 4. Our discussion wrll be confined to the results 

of the br ° k “ “Stance of an additional parameter shows that a bro- 

ken power law model represents a significant improvement over a single power 
law fit for 8 of the 16 observations (5 of the 10 with MFC data) at 5% (MS) 
confidence or better (Table 3). Including PC 1211+143, this lead b , U > ® * 
excesses detected in observations of 8 quasars from our sample of 14. Upper hm 
at 3a are given for the remaining 8 observations for which no significant excess 

was detected. 


2.5 The Effect of Additional Absorption. 

The presence of any additional absorbing material (molecular and partially ionized 
Galactic or material intrinsic to the quasar) along the line of sight to a quasar would 
increase the requirement for a soft excess. We have quantified this effect for the case 
of partially ionized hydrogen in the Galaxy. Reynolds (1989) estimated that 37% of 
the gas along a high latitude line of sight through the Galaxy is ionized. Increasing 
the metal contribution to our X-ray opacities by 37% results in the detection of 
4 more excesses for 3 quasars: (3C273 (9310), 1426+015 (both observations) and 
1613+658). For those quasars with a detected excess, the 0.2 keV flux of the excess 
increased by 20-43%. 

2.6 Notes on Individual Objects. 

PG 1211+143. 

PG 1211 + 143 is already known to have a steep low energy spectrum (Elvis, Wilkes 
&z Tananbaum 1985, Bechtold et al. 1986, Czerny Elvis 1987). This soft excess 
is unusual since it is reported to display a break energy around 2 keV, much 
higher than the 0.3— 0.6 keV generally seen. The MPC data for this object axe 
unreliable, as are all MPC observations made after October 1980, due to a change 
in instrument configuration on that date 1 . The lack of higher energy data combined 
with the 2 keV break energy between the excess and the high energy component 
make it impossible to constrain the spectral parameters and so it is not analyzed 
here. Since it is part of the sample, we include it in our discussion. A detailed 
discussion of its spectrum both from Einstein and EXOSAT data is given by Elvis 
et al. (1991a) and from Ginga data by Turner et a/.(1990a). 

3C273. 

There are three observations of 3C273 over 17 months. The first observation has a 
weak excess (Table 3). In the second, 6 months later, the excess has almost twice 
the flux level of the first and in the third, 11 months later than the second, the 
upper limit is consistent with both the earlier measurements. The high energy flux 
density also varies but is uncorrelated with a 20% decrease over 6 months followed 
by a 93% increase over 11 months. Similar, uncorrelated low and high variations 
were found in Ginga and EXOSAT observations of 3C273 (Turner et al. 1990b). 

PHL 1657 (2135-148). 

The X-ray spectrum is not fitted well by a power law or a black body form. A 
good fit was obtained with a partial covering model (see e.g. Turner and Pounds 

*This change in state of the MPC was not known at the time of the earlier IPC/MPC study by 
Bechtold ei al. (1986) 


1989): a£=2.1±0.5, NH(additional)=7.1±1.6 1 0 20 cm 2 (fraction 92±7%), f(1.0 

keV)=4.1±1.6 /iJy (x 2 =15, n=16). . 

Optical images of this source show two nearby companions and surrounding 

nebulosity (Hutchings et al. 1984) extending to > 6" so that additional X-ray 
absorption is likely to be present. X-ray spatial structure at this level canno 
be confirmed at the resolution of the IPC and indeed the IPC image shows no 
evidence for extended emission. There is a possibility that the closest companion, 
a fainter quasar at the same redshift only 1".9 away, may contribute to the X-ray 
flux leading to the complex spectrum. Two, higher energy, EXOSAT observations 
of this quasar show no evidence for a complex spectrum. A single power law fit 
gives an energy spectral index ~ 0.7 (Singh, Rao and Vahia 1990). However, these 
authors also report a variation in the (0.2-1.0 keV) soft X-ray flux during this 
time which might indicate complexity at lower energies. 

2.7 Excess above an IR— X-ray Baseline. 

Under the hypothesis that the soft X-ray excess is the high energy- tail of the blue 
bump, we would expect its flux to correlate with the strength of the blue bump. 
However, the soft X-ray excess is measured with respect to the higher-energy X-ray 
power law while the blue bump is generally estimated relative to the IR emission. 
While it has been suggested that the IR and X-ray emission is related (Malkan 
1985, Carleton et al 1987), it is by no means proven and the lack of any relation 
could invalidate a search for a blue bump/soft excess correlation due to the large 
continuum range involved. Alternatively, both components could be measured 
above the common baseline of an effective IE/X-ray power-law We have therefore 
estimated the soft X-ray excess above this power law as a check on our earlier 

reSU ’ tS An estimate of the spectral slope from IR to X-ray was made from luminosi- 
ties integrated between 1-2 pm and 0.6 keV using data from Elv.s e! «I.( 199 b, 
and Figure 8), with the exception of 1253-055 (3C279) for which we have lit- 
tle multi-wavelength data. The monochromatic luminosity at 0.2 keV in the rest 
frame, for a soft X-ray excess above the IR-X-ray power law was computed. Such 
an excess is present at > U for 7 of the 8 quasars which show an excess above 
the high energy X-ray power law. Or- quasar, 3C273, no longer shows an excess. 
The fitted results for all -bjects an- upper limits for the remainder are given in 

Table 4. 


3 Results. 

3.1 Soft Excess above the High Energy X-ray Power Law. 

The flux density in the soft X-ray excess at 0.2 keV ranges from ~ 1-6 times 
that in the high energy component at the same energy (Table 4). It is clear y 


very strong, dominating the energy output in the soft X-ray region ( < 0.3keV) in 
quasars. 

The soft excess was detected, at >90% confidence, in 7 of the 13 quasars 
analyzed here and also in PG 1211-1-143. Five of these were reported by Wilkes 
& Elvis (1987), the remaining two, NAB 0205+024 and B2 1028+313, are new 
detections in this analysis. Including PG 1211 + 143, this yields 8 of 14 quasars 
with an excess. Studies of lower luminosity Seyfert 1 galaxies (Turner Sz Pounds 
1989, Kruper, Urry & Canizares 1990) show a similar 50% detection rate. This 
number is likely to increase if, as discussed above, a significant fraction of the 
hydrogen in our galaxy is ionized (Reynolds 1989). Our results confirm that soft 
excess components are common in many types of active galaxies and over a wide 
range in luminosity. 

3.2 Soft Excess above an IR-X-ray Baseline. 

When measured relative to an IR-X-ray power law, 7 of the 8 quasars earlier 
found to have an excess still have one at > 1 a (Table 4). 3C273 is the exception, 
leaving 1 radio-loud and 6 radio-quiet quasars with soft X-ray excesses. Figure 5 
shows the comparison of the excess determined by the two methods with the line 
of equality shown. With the exception of 3C273, the soft excesses determined 
according to the two different methods agree. Although insignificant, the sense of 
any difference is that the excess strength decreases when measured with respect 
to the IR-X-ray power. This implies that the IR-X-ray slope is generally steeper 
than the high energy X-ray power law (see Tables 3, 4). 

A decrease in the apparent strength of the soft X-ray excess when it is 
measured relative to an IR-X-ray power law is expected for radio-loud objects such 
as 3C2/3. They are stronger X-ray sources than their radio-quiet counterparts, 
this is thought to be due to presence of an additional, radio-linked, high energy, 
X-ray component (Zamorani et al 1981, Wilkes & Elvis 1987). This picture also 
predicts flatter IR-X-ray power law slopes for radio-loud quasars, though for our 
small sample the difference is only 2ct: ai RX = 1.18+0.02 (radio-loud), = 1.22+0.02 
(radio quiet). 

3.3 Relation to other Properties. 

The magnitude and errors of our detected flux densities at 0.2 keV suggest that 
a redshift of 0.5 would be sufficient to make it undetectable in most cases. The 
lack of a detected excess in the two highest redshift objects in our sample (z>0.5), 
PKS 0637—752 and 3C279, is consistent with this prediction. 

In the current sample, 2 of the 5 radio-loud and 6 of the 9 radio-quiet 
quasars show excesses. The probability of seeing such a distribution if the paxent 
populations of these two samples Eire the same is 28% (Fisher test), i.e. there is 
no evidence from this small sample that radio-loud and radio-quiet objects are 



different in either the presence or the strength of their excess Pf°P ert ' es - * S 
conclusion does not change if we instead determine the soft excess above the 1R A- 
ray baseline. It is strengthened if we omit the two higher redshift quasars which 

axe both radio-loud. . , 

No relation was found between the soft X-ray excess, measure in ei er 

way, and the high energy X-ray luminosity density (at 1.0 keV) (16% probability 
of a chance occurrence). 


3.4 The High Energy Component. 


Figure 6 shows a comparison of the high energy slopes determined from the broken 
power law fits with those reported by Wilkes and Elvis (1987) for a single power 
law fit with free N H - The results are consistent in all cases, confirming i es an 
Elvis (1987) conclusion that allowing free N H in a single power law fit with the 1PL 
successfully models the soft X-ray excess. We confirm and extend their primary 
result that a range of slopes is present in the soft X-ray region, 0.2-10 keV, with 
the slopes for radio-loud and radio-quiet quasars being a E ~ 0.5, 1.0 respectively. 

A synchrotron self-Compton origin for the a E ~ 0.5 X-ray component in 
radio-loud 'quasars is generally accepted (Zamorani et al. 1981, Wilkes and Elvis 
1987) For radio-quiet quasars and Seyfert 1 galaxies, the current y popu ar re 
flection model, in which the X-rays are reflected off 

surrounding the central continuum source (Lightman and White 1988 Gmlbert 
and Rees 1988), nicely explains all the available X-ray spectra. The mo e pre ic s 
a E ~ 1.0 at low energies, <5 keV, as seen for radio-quiet quasars with Einstein (Ta- 
ble 3), and a broad hump of emission extending from ~ 10 keV to high energy 
(~ 300 keV). It was the detection of this hump by Ginga , for a sample of Seyfer 
1 Galaxies (see Pounds 1990 for a review), which first brought this model to t e 
fore. The intermediate, a E ~ 0.7, 2 - 10 keV slope, frequently reported for Seyfert 
1 galaxies by HEAQ1-A2 (Mushotzky 1985) and EXOSAT (Turner and Pounds 
1989), is the natural result of fitting a single power law to the complex spectrum 
predicted by the reflection model. Data for individual Seyfert Is observed with 
both HEA01-A2 and the Einstein IPC (0.2-3.5 keV) confirm this conclusion in 
that the HEA01-A2 (~ 0.7) slopes are flatter than the IPC (~ 1.0, Kruper, Urry 
and Canizares 1990). Although we see no evidence for a reflection component in 
the current quasar sample, the S/N of the MPC data is so low that it is unlikely 

to be detected. 

3.5 Relation to the Blue Bump. 

The far-infrared to X-ray energy distributions for those quasars with a soft excess 
are shown in Figure 8. These distributions are corrected for Galactic reddening 
and emission lines have been subtracted. The data are being presented in detail 
by Elvis et al. (1991b) and so will be discussed no further here. 


It is clear from the figures ( e.g . NAB 0205+024) that the blue bump gener- 
ally dominates the energy output of these quasars and that the soft X-ray excess 
may well be its high energy tail. However, we found no evidence for a corre- 
lation between the strength of the blue bump (Figure 7), as measured by the 
ratio of blue bump to IR-X-ray power law flux at 0.1 - 0.2 /im (f B Ei), 2 and that 
of the excess, measured by its ratio to the high energy component at 0.2 keV, 
( fxs(0.2 keV)/fnE(0.2 keV), Table 4 ). This is true even when the two higher 
redshift >0.5 objects are removed and when radio-loud and quiet objects are con- 
sidered separately. While this result is disappointing in view of the popular inter- 
pretation of the soft excess as the high energy tail of the blue bump, it is clear 
from the figure that the uncertainties in the soft excess measurement are suffi- 
ciently large to mask any existing correlation. Further uncertainty is added by the 
possibility of variability and the possible presence of intrinsic, absorbing material. 

No firm conclusions on the identity of the soft excess can be made until 
more data, in particular simultaneous UV and soft X-ray data with improved 
X-ray spectral resolution, have been obtained. This will be possible with ROSAT. 

3.6 Comparison with the Czerny— Elvis Model. 

Accretion disk models have generally not attempted to predict the soft X-ray flux 
in addition to the UV/optical regime, primarily due to the large uncertainties in 
the parameters of the inner regions of an accretion disk ■where this emission origi- 
nates. Czerny & Elvis (1987) do include the soft X-ray and successfully reproduce 
the strong soft excess in PG 1211 + 143 using an electron scattering atmosphere 
surrounding the disk itself. In comparison with this result, column 8 (“CE87 ra- 
tio”) in Table 4 gives the ratio (multiplied by 10) of the product energy-luminosity 
of the luminosity in excess of the IR-X-ray power law at 0.2 keV to that at 10 15 
Hz. This is similar to the parameter defined by Czerny & Elvis (1987, their equa- 
tion 28). They predict that, assuming the accretion disk parameters for every 
quasar are the same and their inclinations are random, ~ 15% of quasars will have 
observed values for this parameter > 1.0. Six of the quasars with a soft X-ray 
excess, t.e. 40% of our sample, exceed this value implying that, for this model to 
apply, either the quasars in this sample are predominantly viewed face-on, or a 
range of model parameters is present. 

4 Conclusions. 

Two component spectral fits to IPC/MPC quasar data show that the soft X-ray 
excess dominates the X-ray flux <0.3 keV in half our sample. The fluxes in the 
excess at 0.2 keV cover the range 1-6 times that in the high energy component at 

2 This ratio is strongly correlated with the bump measure, Cuv/Cir, McDowell et al. (1989) 
but is used in its place to be consistent with the X-ray measurement. 


the same energy. The lack of spectral information combined with uncertainty m 
the lhie-of-sighf absorption makes it impossible to draw any definite conclusions 
concerning its identity. The high energy tail of the blue bump remains the m 

attractive candidate and the common presence of both components rega, : ess 

the class of AGN, might be considered as circumstantial support for this pic . 
However, we cannot rule out such extreme alternatives as a —e y new com- 
ponent or a single X-ray component which steepens below ~ 0.3 keV (Band an 

Grindlay 1986, Schwartz, Qian & Tucker 1990). 

The high energy slopes determined for the broken power-law fits presented 

here are consistent whh those reported earlier by Wilkes and Elvis (198 ). Our 
results are fully consistent with the reflection model popularized by recent G.nga 
results for radio-quiet quasars and Seyfert 1 galaxies. For radio-loud quasars 
the flatter (o E ~ 0.5) slopes is consistent with the presence of a syn^ron ^ ' 
Compton (SSC) radio-linked component (Zamoram el of. 1981, Wilkes L E 
1987) which probablv dominates any reflection component in these quasars. 

Measurements of the soft excess relative to an IR-X-ray power law r agree 
well with those relative to the higher energy X-ray in all objects excep -< • 
This confirms that our soft excess measurements are valid for comparison with 
the blue bump flux which is measured with respect to the IR. Some weaken. S 
of the excess flux for radio-loud objects, such as is consis en wi 

additional SSC component and stronger X-ray flux. Measurement with respect to 
the high energy X-ray is the most sensitive way to detect the soft X-ray excess in 
^diotud quasars. It is once again clear that radio-loud and radio-quiet quasars 

should be treated as distinct in X-ray studies. „ 

Comparison of the observed ratio of the luminosity m excess of the IR X- 
ra y power law at 0.2 keV to that at 10“ Hz with that predicted by an accretion disk 
model imply that a range of parameters is necessary to explain the whole sample. 
However, large uncertainties in line-of-sight absorption (Galactic or intrinsic to 
the quasar) and the lack of simultaneous UV data preclude dialled modeling. 
More higher resolution, soft X-ray observations ( e.g . with ROSAT) are c eax v 
necessary! preferably simultaneously with UV data, in order to delineate the blue 

bump unambiguously. 
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6 Figure Captions. 


Figure 1: 68%, 90% and 99% confidence contours for u* and > N« .for smg 
power law (Wilkes and Elvis (1987)-style) fits to observations of 3C273 ' 
numbers: a) 5692 and b) 9310. These are included for comparison with the earl 

results. 


Figure 2: A histogram of IPC relative to MPC normalization for single power 
law fits using the IPC best fitted slope. 


Figure 3: Results of broken power law fits to the IPC plus MPC data. 


Figure 4: Comparison of the 
as determined using a power 
soft component. 


0 2 keV flux density in the soft excess component 
law and a black body as the spectral form for the 


Figure 5: Comparison of the soft X-ray excess determined relative to the high 
energy X-ray power law and the IR-X-ray power law. 

Figure 6: Comparison of the IPC/MPC slopes determined here with the IPC 
only results of Wilkes and Elvis (1987). 


Figure 7: The strength of the blue bump vs. that of the soft X-ray excess relative 
to the higher energy X-ray power law. 


Figure 8: Far-infrared to X-ray energy distributions for those quasars m our sam- 
ple with a detected soft X-ray excess. The data are fully presented and ismssed 
by Elvis et al (1991b). The scale of log pF.vslogn gives the energy output as 

function of frequency. 
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Table 3: Parameters of Spectral Fits. 
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Table 4: Luminosity of the excess : 
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ABSTRACT 


We have measured the X-ray spectra in the 0.1-10 keV energy range for 17 
AGNs (mostly PG quasars) using the EXOSAT Low and Medium-Energy arrays. 
The sample has been selected on the basis of soft X-ray flux, and includes all the 
AGNs in Zamorani et al. (19S1) and Tananbaum et al. (1986) with an Einstein 
flux greater than 6 x 10~ 12 er</s cm~ 2 sec~ l . 

For each object we have determined the best fit spectral parameters in the 
medium energy range (~ 2 - 10 keV), i.e. the power law slope, the normalization 
and the absorbing column density. In this energy range all the spectra are well 
described by a single power law model, with a wide distribution of energy spectral 
indices in the range 0.4 < a < 1.3. The average spectral index is < a >= 0.89±0.06, 
but most of the objects are distributed around o ~ 1.0. The overall distribution 
of spectral indices is significantly displaced toward steeper spectra, when compared 
to that found in the same energy band by Turner and Pounds (19S9),for a large 
sample of hard X-ray selected Seyfert galaxies. This difference can be understood 
in terms of the different selection criteria of the two samples (soft versus hard X- 
ray selection) and it is concluded that the true distribution of hard X-ray spectral 
indices in AGNs is probably wider than previously suggested. 

Inclusion of lower energy data (~ 0.1 — 2.0 keV) shows the presence of signif- 
icant soft excess in G out of 17 cases. This soft excess emission can be modeled by 
either a steep power law with o ranging from 2.7 to 4.2 or the high energy tail of 
a blackbody with temperatures in the range 40-80 eV. The intersection of this soft 
excess emission with the hard X-ray power law is in the energy interval 0.35-0.75 
keV, in the source rest frame. 

A correlation analysis performed on our data set does not reveal any correla- 
tion between the medium energy spectral properties and other physical parameters 
of the sources. However, fits performed in the energy range 0.1-4. 0 keV (compa- 
rable with the Einstein IPC energy range) show a significant correlation between 
radio-loudness and X-ray spectral properties, with the radio-loud objects having, 
on average, flatter slopes than the radio-quiet ones, in agreement with Wilkes and 
Elvis (19S/ ) findings. The fact that the same correlation is not seen in our data at 
higher energy seems to be explained by a correlation between presence and strength 
of soft excess and ”radioquietness": the objects with significant soft excess emission 
tend to be the most radio-quiet quasars in our sample. 

Subject headings: quasars- Seyfert nuclei-X-rays:spectra 


I. INTRODUCTION 

Quasars and Seyfert 1 galaxies are well known to be powerful sources of X-rays 
which constitute an important, and sometimes even dominant, fraction of their total 
energy output. Although the nature of the X-ray emission is still uncertain, it i 
widely believed that it might be associated with the activity induced by a cen r 
compact object, perhaps a massive black hole. The study of the X-ray spectra and 
their time variability, together with an intercomparison of the spectral behaviour 
other portions of the emitted radiation spectrum, may therefore provide important 
clues to understanding the nature of the emission mechanisms and the phvsica 
conditions of the regions surrounding the central compact object. 

In addition, the study of the X-ray spectra of quasars and Seyfert 1 galaxies is 

important to set constraints to the integrated contribution of these objects to the 

extragalactic X-ray background (XRB), whose origin is still the subject of debate 
(e e Setti 1990 for a recent review). Since the average spectium in the - 
energy interval of a close- by. bright sample of Seyfert 1 galaxies is significantly 
steeper than that of the XRB in the same energy interval, it appears that this class 
of objects cannot contribute * dominant fraction of the XRB unless there is a 
positive cosmological evolution in the number density and/or luminosity of Seyfert 
1 nuclei accompanied by a flattening of the average spectrum (m the observer s 
frame) as a function of cosmic epoch (Setti and ttoltjer, 19S9; Fabian et ah. 1990, 

Morisawa et ah, 1990). __ 

Early surveys in the 2-20 keV range of X-ray bright AGNs (mostly low-lumi- 
nosity and low-redshift Seyfert 1 galaxies) suggested the existence of a ‘umversa 
energy spectrum described by a power law of the form 5(E) oc . w iei 

the photon energy and o cr 0.7. The small spread in the observed distribu ion of 
the spectral indices was consistent with the expected dispersion introduced bj the 
observational uncertainties only (Mushotzky 1984). Recently, Turner and Pounds 
( 19S9 hereafter TPS9 ) reported EXOSAT observations of a larger sample of 4S hard 
X-uav' selected AGXs (again mostly Seyfert 1 nuclei) and confirmed .ha. the energy 
spectra in the 2-10 keV range are described by simple power-laws, with a nariou 
distribution (dispersion „ = 0.16) of spectral indices around a mean energy index 
of about 0.7. There are however some significant exceptions (NGC41ol. 3L-.J, 
MKN335, AKN120). which suggest a real spread in spectral slopes o\ei t i< ia> 
band (TPS9 and reference's therein). 

Hard X-ray spectra of quasars have been recently measured with the Japanese 
satellite GINGA in the energy range 2-10 keV (Turner et ah 1989). For tins sn.al 
sample of S objects (5 radio-loud and 3 radio-quiet) they found a mean spectral 
slope of 0.09, but with a relatively high dispersion of about 0.29 In .fact, two 
objects show spectra significantly flatter than the average: a = 0.4S for 3C-.3 and 
o = 0.05 for PG141G-129. For this sample no correlation between radio lou ne.. 

and the 2-10 kcY spectral slope is found. 

The X-rav spectra of a sample of 20 Seyfert 1 galaxies observed over a broader 
energy range (~ 0.15 - 10 keV; IPC and MPC instruments of the Einstein Obsev- 
vatorv ) show the "canonical" average energy index of 0.7 at energies > - ke\ ■ (Ui} 
et ah* 19S9). ;.From the combined IPC/MPC data it is concluded that, although 


a single power law model is generally acceptable, in most cases a power law plus a 
soft excess emission modeled by a blackbody give a significantly better description 
of the data. 

The simplicity of the X-ray spectra of AGNs at hard energies (t'.e. single 
power laws) apparently breaks down at softer X-ray energies. At these energies 
there are at least two effects which may complicate the description of the X-ray 
spectra. The first, intrinsic absorption, has been shown to be important mainly in 
low luminosity AGNs (Reichert et al. 19S5, Lawrence and Elvis 1982). The second, 
soft excess, has been quite clearly seen in the study of the Einstein IPC data of a 
sample of 33 quasars (Wilkes and Elvis 1987, hereafter WE87) and then confirmed 
to be a common feature in EXOSAT data of Seyfert galaxies (TPS9). The same 
Einstein Observatory IPC data (0.3-3. 5 keV) showed (Elvis et al. 1986; WE87) 
that, although the soft X-ray spectra of AGNs (with both low and high intrinsic 
luminosity) have a wide range of power-law spectral indices (—0.2 < n < 1.8), they 
are strongly grouped around a ~ 0.5 for radio loud objects and a ~ 1.0 for radio 
quiet objects. This difference of the soft X-ray spectra as a function of the radio 
properties has been confirmed and extended by the results obtained by Canizares 
and White (1989) on a larger sample of quasars, covering a wide redshift range. 
They also find that the X-ray spectral index shows no evidence of dependence on 
redshift. 

A recent analysis of Einstein IPC observations of a large sample of mostly radio- 
quiet Seyfert 1 galaxies (To objects) performed by Kruper, Urry and Canizares 
(1990) shows a mean spectral index a = 0.81 ± 0.03, significantly steeper than 
the “universal" hard energy slope. Although the number of radio-loud objects 
in their sample is small, their results are consistent with the correlation seen in 
quasars (WEST. Canizares and White 1989) between soft spectral indices and radio 
loudness. A weak correlation between the X-ray power law index and the X-ray 
luminosity (An ~ 0.2 for AlogLx ~ 5 decades) has been found by these authors 
for the sub-sample of well determined spectra of radio-quiet Seyfert 1 and quasars, 
with the higher-luminosity objects having a steeper spectral index. 

Elvis et al. (19SC) discussed various possibilities which may explain the ob- 
served difference in the spectral index distributions between the IPC results and 
those obtained at harder energies. They concluded that, while the lower energy 
range of the IPC is relevant in a fraction of the objects, the most likely explana- 
tion for the observed difference is biases induced by the selection of the samples to 
be compared. To further investigate this effect we have used EXOSAT Low and 
Medium Energy data to determine the energy spectral index of a number of AGNs, 
selected on the basis of soft (i.e. Einstein) X-ray flux. In §11 we define the selection 
criteria of our sample, we give the observation log and describe in detail our spec- 
tral analysis. The results of our analysis are presented and discussed in §111. the 
conclusions are given in tjIY. 


II. OBSERVATIONS AND ANALYSIS 

We present new X-ray data for a sample of X-ray bright Seyfert 1 nuclei and 



quasars (hereafter AGNs) observed with the EXOSAT Observatory. The objects in 
the sample were selected on the basis of soft ( Einstein ) X-ray fl ux - n P ar ^ icu ® 

sample includes all the AGN in Zamorani et al. (1981) and Tananbaum et ■ al. ( ) 

with an Einstein flux greater than 6 x 10 ~ u ergs cm 2 sec , or equivalently 0.- IPC 
counts per sec , which appeared in the EXOSAT database as of May 1990. A total 
number of 17 AGNs satisfy this selection criterion; most of them are PG quasars 
(Schmidt and Green 19S3). Table 1 lists some relevant data for these objects, lhe 
successive columns give : coordinates and alternate names of each object, redshift, 
absolute B magnitude, Einstein flux, luminosity and spectral index, o 0 *, and radio- 
loudness index. The definitions and/or the references for these data are given in the 
footnotes to the table. H 0 = 50 km sec~ 1 Mpc' 1 and g 0 = 0.0 are used throughout 

this paper. 

For all the objects in our sample the EXOSAT data were obtained with both 
the low energy (LE) imaging telescopes (De Korte et al. 1981 )™ d th * med ^ m 
energy (ME) detector array ( Turner, Smith and Zimmerman 1981). Although le 
ME detector was in principle sensitive up to about 50 ke\ , the signal/noise ratio 
becomes very poor above - S keV for the relatively weak sources of our sample. 
Details on the instrumentation can be found in White and Peackock (\9Sb). 

Table 2 gives the observing log for the 17 objects of the sample. The various 
columns list the source name, the date of the observations, the detector used (ME 
or LE filter), the exposure time (in seconds) and, finally, the count rate (in counts 
per second for ME and 1U“ 2 counts per sec for LE observations). M hen more than 
one observation was available in the EXOSAT database for the same objects, we 
have used the data-set with the longer exposure time and/or higher quality. As seen 
in Table X most of the objects have been observed with at least two low energy 
filters (3000 Lexan and Al/p). the only exception being PKS 0037-75 for winch 
there is onlv one observation with the Lexan filter. For ten objects additional LE 
data with the Boron filter are also available, while for one object (PG1229+204) we 
have also one observation with the PPL filter. For all filters (with the exception of 
the boron) we have used vignetting corrected and background subtracted data fiom 
the EXOSAT database. The Boron case is special since the point spread function in 
this filter stronglv depends on the source spectral shape. For this reason a care u 
re-analysis of all the original image data has been performed using the interactive 
analysis system of the EXOSAT Observatory. 

The effective areas of the LE filters can be found in Figure 2 of White and 
Peackock (19SS). While the Boron filter has a narrow energy response peaked at 
around 1 keV, the Lexan and Al/p filters are still sensitive at much lower energies, 
well below the C-band (0.2S keV). The combination of the count rates in the low 
energy filters with the more detailed spectral information given by the medium 
energy detector at higher energies not only provides a better determination of the 
spectral slope, but also allows to constrain the degree of interstellar and/or mtimsic 
absorption and to detect, if present, additional soft components. 


a) ME Spectral Fitting 


For all our objects the ME analysis has been performed only in the energy 
range where a net signal above the background was detected. The highest PH A 
channel used ranges from channel 25, corresponding to about 7 keV, to channel 35 
(about 9-10 keV). For each object the pulse height spectrum was fitted with a simple 
model corresponding to a source power law spectrum corrected for absorption due 
to interstellar gas of solar composition (Morrison and McCammon, 1983). For the 
fitting procedure we used the standard x 2 minimization technique available in the 
EXOSAT spectral fitting package (XSPEC; Shafer, Haberl, Arnaud 1989). 

In a first fitting run both the slope of the power law spectrum and the amount 
of absorbing hydrogen were considered as free parameters and the 90 % confidence 
limits on each parameter were calculated (see Avni 1976; Lampton, Margon and 
Bowyer 1976). Because of the relative faintness of our objects, with typical fluxes 
in the range 4 to S x 10 _12 e/£/s cm~ 2 sec- 1, the allowed region in the Nh - a plane 
is almost always very large (Figure 1). With the only exception of PG0049+171 the 
values of Galac tic Nh arc- well inside the allowed regions and mostly very close to 
the best fit values, so that we have decided to keep Nh fixed at the Galactic values. 
The results of these fits are shown in Table 3 where for each object we give the 
channels used in the fit. the flux in the 2-S keV energy range, the power law spectral 
index with associated 90 ( Z errors (A\ 2 = 2.71 for one interesting parameter, Avni 
1976). the Galactic Nh from the accurate radio measurement performed by Elvis. 
Lockman and Wilkes (19S9). and, finally, the reduced \ 2 and degrees of freedom. 

As seen from the last column of Table 3, the simple power law model with 
absorption from our own Galaxy provides an acceptable fit for essentially all the 
objects (reduced \ 2 of the order of unity). The overall distribution of the reduced 
X 2 for the entire sample (with Nh = Nngal) is consistent with that expected by 
chance according to a Kolmogorov - Smirnov test (p ~ 0.2). Therefore our data 
do not indicate the presence of any substantial amount of intrinsic absorbing cold 
matter. PG0049+171. with a reduced \ 2 of 1.98 (17 degrees of freedom), appears 
to be the only object for which a substantial Nh in excess of Galactic Nh seems to 
be required. For this object we give in the footnote of Table 3 the parameters corre- 
sponding to the best fit with free Nh- Inspection of the data for this object shows 
that an important contribution to the large x 2 is due to a single channel (chan- 
nel 15 at about 4.2 keV) which is significantly higher than the adjacent channels. 
Excluding this "peculiar" channel from the fit, the reduced \ 2 with Nh = Nngal 
becomes 1.31. We also note that PG0052+251 is the only object for which the best 
fit column density is higher than the Galactic value at 68 % confidence interval, 
however this value is consistent with the Galactic Nh at 90 %. 

Thermal bremsstrahlung models were also tried, but they generally provided signif- 
icantly worse fits to the data. 


b) LE and ME Combined Spectral Fitting 

The same simple model (single power law + cold gas absorption) was then 
fitted to the combined LE and ME data. The results of these fits are shown in 
Table 4, where for each object we give two sets of results, the first with both the 
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For some of the objects in our soft X-ray selected sample there are already 
published EXOSAT spectral data with which we can compare our results. These 
objects are MKN335, IIIZW2, 3C273 (TP89), PG0026+129 (Treves et al. 198S), 
MKNS41 (Arnaud et al. 1985). In most cases our fits are in good agreement, within 
the euois, with the published results. This is true both for the slopes of the power 
law and for the presence or absence of soft excess emission. The only discrepancy is 
in the temperature of the fitted blackbody for MKN841; the value found by Arnaud 
et al. (1985) (~ 17 eV) is significantly lower than our best fit value (57 eV). It must 
be noted, however, that Arnaud et al. fitted the blackbody emission combining data 
in the ultraviolet and in the soft X-ray ranges. Additional comparisons with other 
non EXOSAT, results will be discussed in the next section. 


III. DISCUSSION 
a) ME Results 

The hard X-ray emission spectra of all the 17 AGNs in our sample can be well 
described by a single power law. The distribution of the ME spectral indices is 
shown in Fig. 2. where the dashed part of the histogram represents objects with 
a one sigma error on the slope smaller than 0.20. In the same figure we show also 
the gaussian which best fits the data in TPS9. The distribution of spectral indices 
of our objects shows a wide range of slopes centered around o ~ 1.0 and with a 
pronounced tail toward hat slopes. The average spectral index for our 17 objects 
is < a >= 0.S9 ± U.OG. with an observed dispersion a = 0.25. In this case the 
weighted mean is somewhat meaningless because it is dominated by the slopes of 
two objects (3C2<3 and IIIZYV2) which, having a significantly higher ME flux, also 
have 1(7 errors on the slope much smaller than the rest of the objects. Since both 
these objects happen to be among the flattest of the entire sample, the weighted 
mean would be artificially displaced toward flat slopes. 

Our distribution of observed ME spectral indices is significantly different (at 
9S% confidence level, according to the Kolmogorov-Smirnov test) from that obtained 
with the same instrument by TPS9 for hard X-ray selected AGNs. In their sample 
the percentage of Seyfert 1 galaxies with o steeper than 1.0 is 7 %, to be compared 
with 41 % in our sample. Figure 3 shows the spectral indices of the Seyfert 1 and 
quasars in the two samples us a function of the X-ray luminosity. Two differences 
between the samples are immediately clear from the figure: first, our objects are on 
average more X-ray luminous; second, in the range of overlapping luminosities our 
objects tend to have a steeper spectral index (even if a few flat spectrum objects 
are present, also in our sample). Part of the latter difference may be explained by 
the different selection criteria of the two samples, as Elvis et al. (1986) concluded. 
While ours is a flux limited sample selected in the soft X-ray domain, the TP39 
sample comprises objects which have been primarily selected in the hard X-ray 
domain. Each of the two selections has its own biases, ours against flat spectra, 
theirs against steep spectra. 

However, essentially the same difference is present between our data and the 
Einstein MPC data of 20 Seyfert 1 galaxies of Urrv et al. (1989). who claim that 
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b) ME + LE Results 

The addition of the low energy (LE) data has allowed the analysis of the X-ray 
Hie addition soanning about two orders of magnitude (0.1 - 8 
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emission has been also lecovered by Mosnou et al. (1991). Vice versa they find 
soft excess in two out of five of the remaining objects in common for which a single 
power law gives a good fit to our data. 

We note that the objects with soft excess in our sample tend to lie in regions 
of lower Galactic Njj than those which are consistent with a single power law. For 
example, soft excess is detected in 5 out of 8 objects with a Galactic Njj lower than 
3.5 x 10 20 , and in 1 out of 9 objects with Galactic Nh greater than this limit. The 
same effect is present at some level also in the TP89 data: the percentage of Seyfert 
1 galaxies with soft excess in their sample increases to 55 % when the objects with 
Galactic Nh < 3.5 x 10 20 are considered. Some correlation between detection of soft 
excess and Galactic Nh s qualitatively expected because the same amount of soft 
excess is more easily detected in the absence of any kind of intervening absorption. 
To obtain a more quantitative estimate of the possible importance of this effect we 
have performed the following test: first, we have created simulated data for the 
six objects with soft excess emission (Table 5) assuming that the true spectrum of 
each source is the one found from the best fit with a double power law, or with 
the blackbody plus a power law models. For each source we have assumed the 
actually observed background. Then we have applied our fitting procedure on these 
simulated data through the actually used filters searching for the minimum amount 
of Galactic A# above which there would be no statistically significant difference 
between the single and the double power law (or blackbody plus power law) fits: 
had these objects been in a region of Galactic Nh greater than this minimum value, 
we would have classified their spectra as consistent with a single power law. For 
five objects we find that this minimum value of Nh lies between 8 x 10 20 and 
1-^ ^ 10 ( itoms cm 2 . For the remaining object (PG1229-I-204) the minimum 

value of Nh is about 4 x IQ 20 . These minimum values of Nh have been computed 
choosing the appropriate A\ 2 with an associated probability of 95 % according to 
the F-test. (F cr i t ~ 4 for the relevant number of d.o.f.). Since these values (with 
the exception of PG 1229+204) are much higher than the Galactic Nh of all the 
objects in our sample, we can conclude that a strong soft excess, if present, would 
have been detected in the objects classified as single power law, all of which have a 
Galactic A h lower than G x 10 20 atom s cm ~ 2 . A more marginal excess, however, 
(see the case of PG1229+204) could have escaped detection in objects with Galactic 
Nh > 4 x 10 20 . As a consequence, the fact that a large fraction of soft excess objects 
in our sample is located in regions of low Galactic N H appears to be accidental (see 
also §III.c) even if the probability of observing this distribution by chance is low 
(~ 4.2% according to the Fisher test). 

To give a more quantitative estimate of the number of objects for which the 
detection of a weak excess could have been missed by our fitting procedure, we 
have computed for each object the quantity B tx (last Column of Table 4) defined 
as the ratio of the observed minus predicted counts in the low energy filters to the 
predicted counts where the predicted counts have been computed by extrapolating 
to low energy the ME best fit. For the six objects classified “soft excess” the value 
°f R*x is significantly larger than zero, thus confirming the previous classification. 
For six out of eleven objects for which a single power law gives an acceptable fit 



over the entire energy range, we find an excess of counts (he. J?,, , > 0) » 

significance between l<r and 2<x, while for the remaining 5 »e .find a deficit of counts 
(i e Bex < 0). We have then considered the six objects which h « . , 

whose spectra are well fitted by a single power law and have* i, ted then data wdh 
a double power law model with the high energy parameters fixed from the ME fi 
anda steep low energy power law. To be consistent with the results shown m Table 
5 we varied o,„„ in the range 2.5-4.0 leaving as a free parameter the low energy 
normXation or! equivalently, the break energy. The results show that the break 
energies above which the fits become significantly worse than the single power aw 
fits lie between 0.2 and 0.5 keV depending on the adopted steepness of the low 
energy power law and on the value of Jt„. These values have to be compared w,th 
the SeL energies derived for the objects with soft excess emission, which are in 
he energy range 0.35 - 0.75 keV. The partial overlap between the energy range o 
the breaf energies for objects with and without detected soft excess suggests that 
we cannot exclude that also some of the AGNs in our sample whose spec i a at 
consistent with a single power law do really have some kind of soft excess, although 
confined to softer X-rny onoigies* 

Also it should lie noted that the objects with detected soft excess emission tend 
to have better constrained ME spectra, i.e. smaller error on a. For this reason thej 
have on average, smaller errors on the expected extrapolated counts at soft energy 
L Thai "he soft excess, if present can be more easily detected. On the basis of rite 
results of the two previous tests and of this consideration we cone ude that the ob- 
served percentage (35 'XI of soft excess objects in our sample should be consideied as 
Ttowe. limit on the true frequency of such objects. Even more so . we consider that 
the presence of a possible absorption intrinsic to the objects which was assui 
be absent in our fitting procedure, would cause a systematic underestimate of the 
soft excess strength. We have quantified the effect of additional absorbing materia 
assuming an intrinsic column density at the source of 2 x 10 atom* cm ■ 
corresponding increase of the X-ray opacity results in the detection of soft excess 

in one more quasar, he. PG2130+099. 

An independent indication of the possible presence of some excess emission at 
low energies in most of the objects is seen in Fig. 4. which shows the distribution of 
AA’ a dTied as the difference between the best Nh value obtained through the 
oTa single power law to the ME and LE data and the Galactic one. If all the object 
had a negligible intrinsic absorbing column, this distribution would be symmetric 
about o with a spread consistent with the expected small-scale fluctuations of 
Nh in our Galaxv and the measurement errors. Even excluding the six objects 
whh de “ ted soft excess (marked with an S in the histogram), the distribution 
Ts Inificmrtt disphiced toward negative values of AN* only two objects have 
Nh^EXOSAT) > .Y„(21nti). Following WEST, this can be taken as an evidence 
of a gc al presence of some soft excess in the population, although not detectable 
in the ,lvsis of the spectra of the single objects. Moreover, this is consistent with 
the avei .ge steepening of the fits to the ME + LE data with respect to the ME 
ones- for the 11 objects fitted with a single power law we find an average LE+ME 
spectral index < o >= 1.01 ± 0.0S, while for the same objects the mean of the - IE 


spectral indices are 0.95 ± 0.05. 

For the objects with a significant soft excess the resulting parameters from the 
fits with two power laws or a blackbody plus a power law are qualitatively similar 
to each other. Because of the lack of detailed energy resolution of the LE detectors 
the spectral shape of the soft excess is not well constrained, and for most of these 
objects we cannot discriminate between the power law and the blackbody models. 
However, it is interesting to note that for five out of six objects the blackbody 
model gives a better fit than the power law model and the difference between the 
two models is statistically significant (see columns 4 and 8 of Table 5) for 3C273 
and MKN335, i.e. the two sources with the highest soft flux. In particular, for each 
of the objects we find that: 

a) For the two power law model the best fit slope at low energy falls in the range 
2. 7-4. 2, much steeper than the high energy slope. 

b) For the blackbody plus power law model the best fit blackbody temperatures 
are in the range of 45 to SO eY, with associated errors of the order of 30 %. This 
range of temperatures is at variance with the bimodal distribution of the preferred 
temperatures (AT < 10 e\ or AT ~ 200 — 300 eV) of the fitted blackbody suggested 
by Urry et al. (19SD). on the basis of an analysis of IPC data of Seyfert 1 galaxies. 
However the) - note that this result is affected by the discontinuous response of the 
IPC in the energy range 250-G00 eY. 

c) For both fitting models the nominal energy of intersection of the soft and hard 
component is in the range 0.35 - 0.75 keV. 


c) X-RAY SPECTRA VS. RADIO-LOUDNESS 

We have searched for correlations between the hard X-ray spectral index and 
other physical parameters of the sources, but none has been found. Specifically, 
no significant correlation has been found with the intrinsic optical luminosity nor 
with the redshift spanning from 0.025 to 0.25 (with one object at z=0.Gb) nor with 
the X-ray luminosity of the objects, which encompasses about three decades. As 
discussed in §III.a, the fact that the mean spectral slope of our sample, containing 
mostly high luminosity objects, is steeper than the one found by TPS9 for lower 
luminosity AGNs. is probably due to selection effects. 

In agreement with the findings of TP89 no obvious correlation is evident also 
between the ME spectral indices (a me) and the strength of the radio emission (/?£,), 
defined as Rl = log(fR/ J r) where Jr and fs are the flux densities at 5 GHz and 
2500 A (Fig. 5a). Such a correlation was found by WE87, although at lower energies 
(0.3-3. 5 keV), with the radio-loud objects grouped around a ~ 0.5 and the radio- 
quiet ones around a ~ 1.0. Following WES7 we define a quasar- radio-loud when 
i?L > 1. In order to understand the reason for this apparent discrepancy (note that 
11 out of 17 objects in our sample are in common with the WES7 sample) we have 
performed spectral fits to our data in the 0. 1-4.0 keV energy range (similar to the 
IPC range, 0.2-3. 5 keV). The results (Fig. 5b) show that in this softer energy range 
there is a clear trend with radio-quiet quasars having steeper slopes than radio-loud 
quasars, in agreement with the WES7 work. However, for the 11 objects in common 



wi th WE87 we find that our 

2olZ^n F^g. 5b where most of our data points lie above the 
the WEST correlation could be due to either a systematic effect between the 
detectors (IPC and LE/ME), or the slightly softer energy intend 
spectral fit results. Here we point out the existence of four objects with 
FXOSAT steep slopes (o > 1.6) in the 0.1-4.0 keV energy range (PG0003+1SW, 
PGini+143 PG1229-* 204, PG2130+099). One of them (PG1211+143) was the 
Llv Object With a very steep slope (a = 1.8) also in the IPC energy range and was 
classified as exceptional by WES7 . Our results indicate that these steep slopes are 
probably a diffuse feature in the quasars spectra. 

The other apparently significant correlation in our data is between the presence 
of Jt exct a“d radio Lies, Five out of six objects with detected soft excess 
emission are radio quiet and four of these are among the last five most radio- quie 
object in the sample (see column 8 in Table 1), The only radio-loud object with 
soft excess (3C273) is somewhat peculiar with respect to the other objects m b 
;le Lass in tha, its value of ft, is significantly lower than the mean value o h 
other five objects (7?„ - 3). so that the detection of the soft excess in _ 
mainly due to its high X-ray flux, about one order of magnitude larger than the 
average flux of the other objects in the sample. 

In the previous Section we have discussed the apparent correlation between 

the soft excess and the amount of Galactic A'h- Therefore t e S ° aut l 

. + T? \ annpar<; to be COITclflteQ Wltll both i' H 

eauivalentlv. the excess counts H tx ) appears to pe Luimaicu 

ft. In order to determine which of the two correlations is the fundamental one 

have computed the Kendall partial rank-order correlation coefficient T xy x (bieg 1 

and Castellan 19SS). which allows one to compute the probabilities of correlation for 

^ ^ of variables with the other one held constant. The mrly s,gm leant (anti- 

correlation resulting from this test is that between the excess of counts /?„. and the 

radio-loudness, Rl- This result is consistent with the > conclusion $ ' ’ 

that the apparent correlation between soft excess and A h is mosth accidental. 

The different spectral behaviour of the radio-loud and radio-quiet quasars 
also shown in Figure G where we plot the ME spectral slopes versus the 0.1-4. 0 
slopes Two facts are immediately clear from the figure: first, the radio-loud objec 
are spread around the line of equal spectral indices and do not show any chang 
of slope over the entire energy range; second, the radio-quiet o jects en o » 
systematically steeper slopes in the 0.1-4.0 keV energy range. This is true also for 

those objects without soft excess. 

The observed anti-correlation between ft, ft. and the data shown in Fig- 
ure 6 can be interpreted in two physically different ways by adopting simple two- 
comnonents models. Either the radio-quiet objects may have an enhanced soft 
emission on the top of the power law component, which m this case would lie the 
same in both radio-quiet and radio-loud quasars; or, alternatively, the radio-loud 
ones mav have a relatively stronger power law component with a cqstant s ope 10 
tow to nhdium energy (see triangles in Figure 0) The existence of tta enhanced 
power law component would make more difficult the detection of the 


their spectra and would be consistent with the fact that the radio-loud quasars are 
stronger X-ray emitters (about a factor 3) than the radio-quiet ones (Zamorani et al. 
1981). These alternative pictures are perfectly self-consistent within our EXOSAT 
data sample. However if the explanation of the systematically steeper 0.1-4. 0 keV 
spectral slopes found in our sample is not due to instrumental calibration differences 
between EXOSAT and the Einstein IPC, but rather only to the softer spectral ex- 
tension of EXOSAT, then the emission spectra of the sources become much more 
complex and the simple modelling described above need to be modified. 



IV. CONCLUSIONS 


We have analyzed simultaneous LE and ME EXOSAT data to Mta 

X-ray spectra of 17 objects selected by their soft X-ray flux (mostly PG Quasars 
and Seyfert Is nuclei). The most important results of our spectral survey can 

summa^l^Uotvs : ^ ^ ^ _ 2 keV) of the AGNs in our 

sample can be described by a single power law model. These measnredspectra 
indices show a wide variety of power law slopes with a spanning from 0.4 to 1 . 6 . 
The average slope is < a >= 0.89 ± 0.06, but the distribution ^centered around 
o = 1.0. The narrow observed range of slopes around a ~ 0.< in high energy 
selected samples is a result of their tendency to select the most -10 ^ bright 
AGN, which tend to be those with the flattest X-ray spectra over the whole 0.2-10 
lceV range. These results suggest that the true distribution of hard X-ra\ spec i 
indices is wider and probably steeper than previously suggested. . , 

2) Approximately 1/3 (G out of 17) of the sources in our sample show e'iden 
for soft excess emission in the low energy filters (below 1 keV). For these objects a 
power law plus a soft excess modeled by steep power law or by the high eneigj 

of a blackbody is a significantly better description of the data than a single powei 
law over the Energy range 0.1-10 keV. Although the poor energy resolution of the 
LE filters precludes an accurate description of the spectral shape of the soft excess 
emission the blackbody model seems to be preferred to the two power la " _ 

3) The best fit slopes of the steep power law fits at low energies are in the range 
a _ 9 7 _ 4 9 while the preferred temperatures of the fitted blackbodies are in the 
range 40-S0 eV. The "break" energy between the soft and hard component is m the 

range 0.35-0.75 keV at the source. . , 

4) There is no obvious dependence of ME energy index with the other physical 
parameters of the sources such as the optical and X-ray luminosity, the redshift. 

and the radio- loudness. . . 

5) Spectral fits in the 0.1-4.0 keV energy range confirm the correlation between 

the X-ray slope in this band and the radioloudness, with the radio-loud quasars 

having flatter slopes than the radio-qiuet ones. 

6) Five out of six objects with detected soft excess emission are radio-quiet. 

The detection of soft excess emission in one radio-loud object (3C_<3) is pioba 
due to its much higher X-ray flux. 
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Table 1 

THE SAMPLE 

Coordinate 

Name 

t 

Mb 

El 

mm 


151 

■39 

0003+199 

PQ, Mkn 335 

0.025 

-22.55 

76.6 

2.0 

— 

1.25 

0.55 

0007+106 

PG, 111 Zw 2 

0.089 

-22.93 

19.0 

IKQ 

mn 

1.06 

2.24 

0026 + 129 

PG 

0.142 

-25.13 

ED 

BZ1 

dsh 

1.43 

0.34 

0049+171 

PG, Mkn 1148 

0.064 

-22.45 

ED 

mn 

DSfl 

1.24 

-0.49 

0052+251 

PG 

0.155 

-24.95 

7.1 

Bra 

DSfl 

1.36 

-0.61 

0637 - 752 

PKS 

0.651 

-27.65 

6.9 

230.0 

mn 

1.25 

3.52 

0804+761 

PG 

0.100 

-24.30 

DO 

5.8 

— 

1.34 

-0.21 

0923+129 

PG, Mkn 705 

0.029 

-21.70 

EZ3 

Hi 

— 

1.41 

0.32 

1211+143 

PG 

0.085 

-24.25 

32.3 

11.8 

DSfl 

1.22 

-0.72 

1217+023 

PKS 

0.240 

-24.60 

8.6 

26.0 

DSfl 

1.14 

2.38 

1219+756 

Mkn 205 

0.070 

-22.95 

17.0 

2.9 

DSfl 

1.25 

-0.31 

1226+023 

PG, 3C273 

0.158 

-27.43 

83.4 

ESO 

DSfl 

1.32 

3.06 

1229+204 

PG, Ton 1542 

0.064 

-23.58 

Bll 

1.5 

— 

1.47 

-0.96 

1426+015 

PG, Mkn 1383 

0.086 

-23.90 

16.0 


DSfl 

1.26 

1 

O 

ON 

1501+106 

PG, Mkn 841 

0.036 

-21.93 

28.5 

■23 

DSfl 

1.16 

-0.67 

2130+099 

PG, II Zw 136 

0.061 

-23.75 

8.8 

1.8 


1.46 

0.37 

2135 - 147 

PKS, PHL 1657 

0.200 

-25.18 

7.1 

17.0 

DSfl 

1.30 

1.70 


Note*: 

* unit* of 10 ~ lt ergs *ee“ , em”* in the energy range 0. 3-3.5 keV. Zamorani et al. (1981), 
Tananbaum et al. (1986). 

* units of 10 44 erps sec” 1 in the energy range 0.5-4.5 keV (the same references as Sx)> 

* from Wilkes and Elvis (1987), except for PG0049+171 (from Elvis et al. 1986) and 
PG0052+251 (from Elvis et al. 1991). 

* defined in Zamorani et al. (1981). 

* from Wilkes and Elvis (1987), when available, or computed from Keller man et al. (1989) 

data. 
























































































































Table 2 


Observations Log 


Object 1 | Pate* | detector j Exposure Time* 


ME 

31x 

0003+199 84.341 


0007+106 


0026+129 64.198 



Boron 


ME 

85.354 Six 


0049+171 85.334 


0052+251 85.354 


0637*752 84.252 


0804+761 85.075 


0923+129 85.129 


1211+143 85.164 


Boron 


ME 


Al/p 


ME 


Six 


Al/p 


ME 


31x 


ME 


Six 


Al/p 


Boron 


MB 


Six 


1217+023 84.029 


Boron 


ME 


Six 


Al/ 


Boron 


ME 


Six 


Al/p 


2 


20100 


18888 



22170 


10841 


18030 


14340 


26090 


6563 


10735 


12220 


S110 


6568 


1S026 



17880 


8052 


Count Rate* 


1.25 ± 0.03 


20.00 ± 0.90 


9.39 ± 0.56 


0.93 ± 0.09 


2.60 ± 0.04 


8.50 ± 0.20 


2.65 ± 0.43 


.51 ± 0.04 


1.61 ± 0.22 


1.16 ± 0.18 


0.43 ± 0.15 


.71 ± 0.04 


1.43 + 0.12 


0.82 ± 0.21 


.30 ± 0.04 


1.65 + 0.16 


.92 ± 0.17 


0.50 ± 0.08 


0.81 ± 0.03 


S.30 + 0.30 


2.50 ± 0.30 


.35 ± 0.18 


1.02 ± 0.08 


2.70 ± 0.40 


1.80 + 0.20 


.43 + 0.15 


.88 + 0.04 


.2 


8.8 _ 0.62 


0.88 + 0.25 


0.44 + 0.09 


2.54 + 0.28 


1.17 + 0.16 































































































Tkble 2 (cost.) 

Observations Log 

Object* 

Dai It* 

Detector 

Ewpoevre Time* 

Count Rate t 



ME 

mo 

0.57 ±0.07 

1219+756 

64.027 

1 81x 

1679 

2.94 ± 0.51 



Al/p 

6293 

1.34 ± 0.21 

- 


ME 

22820 

6.27 ±0.04 


DE 4 «6 

Six 

6424 

26.00 ± 1.00 


50. loo 

Al/p 

6332 

13.30 ± 0.82 



Boron 

10546 

2.20 ± 0.19 



ME 

49730 

0.40 ± 0.02 



Al/p 

11194 

2.95 ± 0.21 

1229+204 

65.159 

Boron 

21197 

0.56 ± 0.17 



41x 

2624 

6.21 ± 0.57 



PPL 

3156 

6.97 ± 0.58 



ME 

21610 

0.83 ± 0.04 

1426+015 

85.187 

Six 

4689 

4.63 ± 0.39 

Al/p 

3599 

2.71 ±0.35 | 



Boron 

10056 

0.27 ±0.12 



ME 

19860 

1.14 ± 0.04 

1601+106 

86.053 

31x 

1998 

9.55 ± 0.86 



Al/P 

S713 

3.92 ± 0.41 



ME 

26130 

0.49 ±0.03 

2130+099 


Six 

4090 

5.64 ±047 

50*010 

i 

Al/P 

S325 

3.26 ± 0.39 



Boron 

15327 

0.54 ± 0.17 



ME 

28560 

0.71 ± 0.05 

2135-147 

64.133 

Six 

4863 

2.14 ± 0.27 

Al/p 

8033 

1.20 ± 0.16 



Boron 

11691 

0.34 ± 0.12 


Notwi 
Year and Day 
Secondi 

Count* per second (unit* of 10“* tot ike LE filter*). 




Table 3 


Single Power Law fit* to the ME Data 
ngautting Qalactie^yw 


Object* 


0003+199 

0007+106 

0026+129^ 


C hanne l* 

7-35 


0049+171^ 


0052+251 

0637-752 


1217+023 


1219+756 


1229+204 

1426+015 


1501+106 


7-29 


0923+129 7-24 


1211+143 7-25 


2130+099 7-30 


2135-147 7-30 


K..i, xt IJ±LL 


12.6 


1.U-0JI 


0. o5 w a. or 


4.1 


i.03lS;» 
0.931JU. 


4.00 1 1-01/27 
0.75/27 
0.07/17 
1.98/17 
1.26/18 


0804+761 7’95 83 


ffl 

O.VZ.iijs 


ft 7^®'^ 
0 * 75.1144 


~4 80 I 0.61/21 

0.93/17. 
0.98/16 


3.12 


4.03 


2.83 1 1-36/IT 


4.3 


1.97 


6.0 


2J4 0.84/20 


1226+023 7-35 88,8 _ 


ft 

O.53_ 0t 01_ 


1.80 


ft 

O.4o. 0 rg 


1.07/21 

0.81/27 


A 1A+0.H 

0.42_ 0 


2.23 


4.20 


O.B3i; M 


4.45 


Note.: 

UniU ° 1 10 ~ l3 ‘ r *' ** ^ Id Wilke. (1989), except for 0003+199 

■ unit, of 10 ,9 otom. cm *, from Elvi» i>ocam»B 

;TP89) f and 1226+023 (WE87) 

- <U t ree. ot Mom w % ^fidcac* b> the 

1 Pox POtKHB+171 the Celeettc ATjf I. ex u e 

. •. V- free to wry, the beet fit peremetere ere o _ 3 5 ood 

fi t to the dot. Loorm, *■ *" » ^ l# d .,. (ee, di.cu .,00 to the 

~ 8 x 10 ** at cm. cm"*, with a reduced r 01 

te«t). 



















































































EXjSglE 
EEZBUI 


Tkble 4 


Single Power Low fits to (he ME+LE Data 


xl/d.o.f} A i?fn.o.J e Clarification 


Soft Excess 




lMtl A A 



OMi!:** 


fli+O.M 


1-0712;*: 


mi 



1.021*;** 


o O 40 U 


mm 

nssn 


msrnt 


1.1015!! 


mfflfflc 
DBSttll 


«i. 


.0-0.007 


39.75/3 


167.05/3 


4.86/2 


.08/2 


9 1.04/3 


20 | 1.21/3 


13.66/2 


4.26 1.90/19 2.45/2 


2.1±f*» 1.22/19 3.65/2 


4.50 I 1.25/20 | 2.01/2 


.61/21 


.62/22 0.80/1 


1.18/19 6.72/3 


3.12 I 1.14/20 6.99/3 


.94/18 I 1.12/3 


4.03 0.91/19 1.57/3 


2.63/19 


2.83 I 8.12/20 | 38.99/3 


.96/18 | 0.12/2 


.92/19 0.32/2 


. 83/21 


.81/22 I 1.10/2 



Single P.L. 


Single P.L. 


Single P.L. 


Soft Excess 



1.05/24 I 2.64/4 


2.58 I 1.38/25 11.97/4 


.0-0.023 


2.64 


.0-0.051 


2.23 


Soft Excess 


.41/2 


1.97/3 


.90/25 | 3.93/3 



•vo« 0 «* 


3 


1.81/3 


Soft Excess 


Single P.L. 


Single P.L. 


2.7412:" 


0.0012” 


Single P.L. 0.4912$ 


Single P.L. -0.3912$ 


1 14 +1 " 
*•**-0.™ 


Single P.L. -0.4512$ 


0 42■ f0 •* 4 
0.3T 


Single P.L. 0.5112$ 


3.8712$ 


Single P.L. —0 2712$ 


Single P.L. -0.2812:52 


Soft Excess 1 .0012$ 


1 Rfi* 0 ** 1 


Soft Excess 3.2612$ 


% 64 +# m 


1 ID * 1 '* 3 











































































































































Notes: 


• units of 10*° atom* cm" 3 
t degrees of freedom 

• drop on x* from ME to ME+LE peer the »>“*«' <* “ “*"* 

n „tio rf the obeereed minu. predicted count. to the prodieted eouau In the Al/p •»<> 
3000Lex filter.. The predicted count, nre computed by entropoUtinj the ME beet fit. The 

errors are 1 c. 


Table 6 


DOUBLE POWER LAW AND BLACKBODY + POWER LAW FITS 


Object 

! A|(l(l 

Bizga 

Ezsn 

EP1 

ltT(tV) 

nia 

P(x’) 


0003+199 

liiMl 

bmshi 

ESI 

0.46 

46t*® 

1.01/30 

0.45 

0.49 

1211+143 

cm 

BE3SH 

oza 

0.67 

BSM 

1.18/20 

EH 

0.61 

1226+023 

bissi 


ESI 

0.95 

ESfl 

■EJEJB 

EH 

0.50 

1229+204 

iessi 


ESI 

0.46 

Krai 

WXEEM 

E3I 

0.44 

1426+015 


KZ3EB 

ESI 

0.66 

■isa 

msnnm 

EB 

0.76 

1601+106 

JiZsuL 

■EXES 

ESI 

0.64 


WS3EM 

ESI 

0.54 


Notes: 

* degrees of freedom 

* keV (at the source) 









































FIGURE CAPTIONS 


The x 2 contours for the power law spectral fit to the ME data. The X axis 
is the logarithm of the absorbing column density of neutral hydrogen (units of 
10 22 atoms cm' 2 ) and the Y axis is the power law photon index. The contour 
are 68 % 90 %, 99 %, joint confidence limits. The vertical line marks the column 
density expected due to absorption in our galaxy as measured by 21 cm emission 
(Elvis, Lockman and Wilkes 1989). 

F'* 2 

Distribution of the hard X-ray (2-10 keV) spectral energy indices of our sample. 
The dashed part of the hystogram represents objects with la error on the s^pe 
smaller than 0.20. The gaussian curve is taken from Turner and Pounds (1989) ai 
describes the best fit energy index distribution of their sample. 

ME spectral energy index values versus the logarithm of the X-ray luminosity 
(er<,s sec' 1 )- We have marked with filled circles the objects of out r sample , with 
open circles the Seyfert 1 galaxies and quasars in the Turner and Pounds (19S9) 
sample and with open triangles the three objects in common. 

Distribution of A.Y H (in units of 10 20 atoms cm' 2 ) defined as the difference 
between the N H value obtained through the fit of a single power law over the full 
enercv range (~ 0.1 - 10 keV) and the Galactic value obtained by radio measuie- 
ments (Elvis, Lockman and Wilkes 19S9). The objects with detected soft excess 

emission are marked with an S. 

F'* 5 

a ) ^ME spectral energy index values versus the radio-loudness index 

b) 0 1-4 0 keV energv index values versus the radio-loudness index Hi- 

We have marked with open circles the objects with soft excess emission and 
with filled circles the others. The line gives the WEST correlation. 

Pjg t 0 

ME spectral energy index values versus the 0.1-4.0 ke\ ones, 
triangles represent radio- quiet and radio- loud objects respectively, 
with soft excess are represented with open symbols. 


Circles and 
The objects 















FIGURE 2 






FIGURE 4 
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FIGURE 5 
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